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EXECUTIVE  SUMMARY 


This  report  describes  the  theory  support  of  DTRA’s  Plasma  Radiation  Source  (PRS)  program 
carried  out  by  NRL’s  Radiation  Hydrodynamics  Branch  (Code  6720)  in  FY  2003.  Included  is 
work  called  for  in  DTRA  MIPR  03-2045M  -  “Plasma  Radiation  Theory  Support”.  This  work  was 
also  partially  supported  by  Sandia  National  Laboratories  as  part  of  its  high  energy  density  physics 
campaign. 

This  year  Code  6720  assumed  partial  responsibility  for  coordinating/planning  DTRA’s 
Analysis  and  Modeling  Working  Group  meeting.  As  part  of  this  effort  a  meeting  was  organized  and 
held  on  May  22, 2003  in  Colorado  Springs.  The  purpose  of  the  meeting  was  to  review  and  discuss 
the  physics  and  engineering  issues  affecting  long  implosion,  12  cm  diameter  nozzle,  experiments 
on  DE  and  DQ.  In  addition,  recommendations  for  the  priority  issues  to  be  studied  in  the  series 
of  DE  and  DQ  experiments  to  be  performed  in  the  later  half  of  2003  were  made.  The  minutes 
and  recommendations  of  this  meeting  are  reported  in  Section  (I)  of  this  report.  The  remaining  six 
sections,  three  of  which  (III)  -  (V)  have  been  submitted  for  journal  publication,  are  summarized  in 
the  proceeding  paragraphs,  as  follows: 

(II)  To  help  achieve  the  near  term  goal  of 40-60  kJ  of  Ar  K-yield  on  Decade  Quad  the  2003  DQ 
experiments  are  analyzed  from  a  OD  and  ID  model  perspective,  beginning  with  a  discussion  of  the 
analyzes  that  was  provided  to  the  PRS  community  prior  to  the  2003  DQ  shots.  These  analyzes  are 
then  updated  to  account  for  the  reduced  current  that  was  delivered  to  the  experimental  loads.  There 
is  substantial  contrast  between  the  two  perspectives.  The  OD  model  supports  the  contention  that 
the  12  cm  diameter  nozzle’s  performance  on  DQ  was  significantly  degraded  from  its  performance 
on  DE.  On  the  other  hand,  the  ID  model  assessment  is  that  the  12  cm  diameter  nozzle  performed 
similarly  on  both  machines.  In  general  the  ID  yield  results  followed  the  same  trend  as  found  in 
the  experiments.  Both  models  show  that  in  excess  of  9  kJ/cm  of  K-shell  emission  is  potentially 
achievable  with  a  12  cm  diameter  DQ  nozzle  at  the  present  level  of  current  delivery  from  DQ. 

(III)  Recent  experiments  conducted  on  the  Saturn  pulsed-power  generator  at  Sandia  National 
Laboratories  have  produced  large  amounts  of  x-ray  output,  which  cannot  be  accounted  for  in 
conventional  MHD  calculations.  In  these  experiments,  the  Saturn  current  had  a  rise-time  of  180 
ns  in  contrast  to  a  rise-time  of  60  ns  in  Saturn’s  earlier  mode  of  operation.  In  both  aluminum  and 
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tungsten  wire-array  Z-pinch  implosions,  2  to  4  times  more  x-ray  output  was  generated  than  could 
be  supplied  according  to  one-dimensional  magneto-hydrodynamic  calculations  by  the  combined 
action  of  j  x  B  acceleration  forces  and  Ohmic  heating.  In  this  section,  we  reexamine  the  problem 
of  coupling  transmission  line  circuits  to  plasma  fluid  equations  and  derive  expressions  for  the 
Z-pinch  load  circuit  resistance  and  inductance  that  relate  these  quantities  in  a  1-D  analysis  to 
the  surface  resistivity  of  the  fluid,  and  to  the  magnetic  field  that  is  stored  in  the  vacuum  diode, 
respectively.  Enhanced  energy  coupling  in  this  analysis,  therefore  comes  from  enhancements  to  the 
surface  resistivity,  and  we  show  that  plasma  resistivities  approximately  three  orders  of  magnitude 
larger  than  classical  are  needed  in  order  to  achieve  energy  inputs  that  are  comparable  to  the  Saturn 
experimental  x-ray  outputs.The  challenge  to  the  DTRA  program  is  to  channel  this  energy  into 
enhanced  x-ray  outputs  that  are  of  long  term  interest  to  the  program. 

(TV)  Confinement  and  compression  of  magnetic  flux  by  plasma  shells  is  of  interest  for  a 
variety  of  applications  associated  with  keV  x-ray  production  in  pulsed-power  driven  Z-pinch 
plasma  radiation  sources  (PRS).  Confinement  of  a  field  with  the  field  of  the  pinch  current  can  help 
stabilize  the  implosion  of  a  shell  from  a  large  initial  radius.  Compression  of  the  azimuthal  magnetic 
flux  in  the  PRS  (secondary)  circuit  with  a  plasma  shell  driven  by  field  of  the  primary  circuit  may 
represent  a  new  opportunity  for  using  low-cost,  relatively  slow  pulsed  power  to  generate  large 
amounts  of  keV  x-rays.  The  magnetic  field  has  to  be  compressed  and/or  confined  by  low-mass 
plasma  shells  emitting  soft  x-ray  radiation  that  limits  the  temperature  and  conductivity  of  the  shell 
plasma.  The  thickness  of  a  plasma  shell  is  established  self-consistently  during  the  implosion,  and  it 
is  not  obvious  that  it  becomes  thick  enough  to  confine  or  compress  the  magnetic  flux.  The  results 
of  analytical  theory  and  numerical  simulations  demonstrate  that  the  flux-compressing  capability 
of  a  low-beta  plasma  shell  is  surprisingly  good  because  the  shell  is  shown  to  dynamically  adjust 
its  thickness  so  that  it  always  remains  of  the  order  of  its  skin  depth.  The  self-similar  profiles  of 
confined  magnetic  field  predicted  by  the  theory  are  consistently  reproduced  in  the  simulations. 

(V)  Imploding  a  few  high-Z  wires  at  large  radius  onto  a  massive  wire  core  is  another  idea  for 
producing  copious  amounts  of  >  5  keV  x  rays  being  pursued  to  meet  DTRA’s  long  term  goals.  In 
this  section  energetic  implosions  using  two  or  three  load  wires  to  create  a  focused  axial  stagnation 
of  dense  wire  cores  midst  the  assembled  precursor  plasma  are  examined  with  respect  to  the  trade- 
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olf  between  the  implosion  mass  lost  to  precursor  ablation  and  the  mass  or  kinetic  energy  available 
at  stagnation.  The  calculated  kinetic  energy  at  stagnation  serves  as  the  primary  source  for  the 
output  x-radiation. 

(VI)  The  12  cm  diameter  argon  experimental  shots  taken  during  the  October  2003  campaign 
on  DQ  all  exhibited  substantial  evidence  of  current  flashover  at  the  insulator  stack.  In  this  section, 
the  effect  of  this  current  flashover  on  K-shell  yield  production  is  analyzed  using  ID  and  2D 
radiation  MHD  models.  From  this  analysis  we  conclude  that,  in  terms  of  the  general  energetics 
of  the  flashover  event,  the  K-shell  yield  attained  on  DQ  was  not  significantly  reduced  by  the 
presence  of  flashovers.  This  is  not  to  say  that  azimuthal  perturbations  and  a  more  pronounced 
r-z  perturbation  than  were  modeled,  perhaps  brought  on  by  an  asymmetric  power  flow  from  the 
flashover  or  tube  misalignment,  could  not  account  for  the  observed  yield  performance. 

(VII)  This  section  is  available  as  a  separate  document,  upon  request  from  Code  6720.  It 
is  the  Final  Report  from  NumerEx  to  the  Naval  Research  Laboratory  (contract  No.  N00173- 
03-P-6201.  The  title  of  this  document  is  Gas  Puff  Nozzle  Design  for  Argon  Z-pinches  Using 
Mach2  Simulations.  The  report  discusses  the  improvements  NumerEx  has  made  towards  Mach2’s 
capability  of  simulating  the  implosion  dynamics  of  a  Z  pinch.  These  improvements  include:  1)  A 
grid  boundary  shape  created  by  interpolating  along  a  discrete  curve  read  from  a  file.  2)  Inclusion 
of  a  tabular  collisional  radiative  equilibrium  model  for  equation  of  state  and  radiation  transport. 
3)  Corrections  to  the  energy-based  voltage  model  for  coupling  the  external  circuit  to  the  plasma 
and  diode  region  for  Lagrangian  problems,  and  4)  Determination  of  the  proper  plasma-to- vacuum 
coupling  scheme  for  enhanced  circuit-to-simulation  coupling  fidelity  for  Eulerian  simulations. 
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I.  Key  Issues  and  Future  Directions  in  PRS  Research 


On  May  22,  2003,  a  DTRA-sponsored  meeting  was  held  in  Colorado  Springs  to 
discuss  outstanding  issues  in  PRS  physics  and  possible  means  of  resolving  and 
clarifying  those  issues  via  collaboration  between  experimental  and  modeling  efforts. 
The  meeting  was  chaired  by  J.  Apruzese  of  NRL  Code  6720.  The  following  is  a 
summary  of  the  meeting  and  the  group ’s  findings. 


The  emphasis  of  the  PRS  Working  Group  meeting  was  to  share  common  experiences 
being  gained  in  physics  and  technology  of  long  implosions  from  large  radii  (within  the  technical 
constraints  of  Double  Eagle  (DE),  which  limits  K-shell  producing  implosions  to  about  260  ns 
and  no  more  than  12  cm  outer  diameter)  and  to  formulate  recommendations  for  the  future  work. 
Its  purpose  was 

1)  Review  and  discuss: 

-  the  engineering  and  physics  issues  to  be  addressed  in  long  implosion 
experiments  on  DE  and  Decade  Quad  (DQ)  with  12-cm  outer  diameter 
nozzles  and  in  the  supporting  modeling  efforts, 

-  the  design  aspects  of  the  15  cm  diameter  nozzle  and  its  testing  on  DQ; 

2)  Develop  recommendations  for  the  priority  issues  to  be  studied  in  a  new 

campaign  on  DE . 

The  review  part  included  informal  status  reports  and  presentations  of 

-  Henry  Sze  on  the  TPSD  12  cm  diameter  nozzle  design  and  performance; 

-  Phil  Coleman  on  pre-ionization  issues 

-  Bob  Terry  and  Mike  Frese  on  the  status  of  2-D  RMHD  modeling  in  support  of 
DE  and  DQ  experiments 

The  discussion  was  focused  on  the  following  main  topics,  which  are  listed  below. 

Their  order  corresponds  to  the  community’s  suggestions  for  a  better  understanding  for 
the  progress  of  the  PRS  program  in  the  immediate  future,  as  established  by  consensus 
which  involved  voting  (the  voting  results  are  listed  at  the  end  of  this  report): 

1 .  Determination  of  the  optimum  mass  ratio  of  shells  and  central  jet  (more 
generally,  optimmn  radial/axial  density  distributions)  for  K-shell  yield 
improvement. 

2.  Determination  of  how  much  of  the  K-shell  yield  comes  from  the  atoms 
injected  from  each  plenum  (inner/outer  shell,  central  jet),  the  role  of  inter¬ 
shell  mixing. 

3.  Determination  of  the  fate  and  role  in  the  implosion  of  the  “dark  matter,”  the 
plasma  mass  which  does  not  radiate  in  the  K-shell  spectrum. 

4.  The  use  of  the  lower-Z  gas  in  the  outer  shell  to  improve  yield. 

5.  Status  and  the  need  for  pre-ionization. 

The  main  recommendation  of  the  PRS  Working  Group  to  DTRA  is  to  provide  for 
an  (at  least)  7-week  experimental  campaign  on  Double  Eagle,  with  associated  simulation 
and  analysis,  to  address  the  above  issues  in  the  order  of  their  priority.  There  was  a 
general  agreement  among  the  participants  of  the  Working  Group  that  a  completion  of 
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this  study  with  12-cm  nozzles  is  a  must  before  meaningful  steps  towards  designing  14-16 
cm  diameter  nozzles  for  DQ  and  DE  could  be  taken. 

Some  other  issues  addressed  in  the  discussion: 

-  the  role  of  (enhanced)  resistive  component  in  energy  coupling  to  the  pinch  and 
its  effect  on  the  K-shell  yield; 

-  additional  requirements  to  the  load  design  (including  pre-ionization)  and 
performance  due  to  the  constraints  imposed  by  the  debris  shielding  task; 

-  diagnostics  (fielding  of  new  diagnostics,  restoration  of  some  that  already 
exist). 


H.  Sze  reviewed  the  12  cm  nozzle  design  and  its  characteristics.  Recently 
completed  measurements  of  gas  flow  in  the  12  cm  nozzle,  with  central  jet,  were 
described.  Combining  gas  flow  timing  with  a  central  jet  and  an  appropriate  recess  of  the 
nozzle  with  respect  to  the  cathode  plane  appears  to  be  sufficient  to  control  zippering, 
which  produces  a  tighter  pinch  (and  a  higher  power),  compared  to  other  2-shell 
implosions.  High-quality  Z-pinches  with  short  FWHM  x-ray  pulse  widths  observed  on 
DE  in  recent  experiments  demonstrate  significant  achievement  in  the  technology.  H.  Sze 
indicated  that  although  300+  ns,  12  cm  diameter  are  not  optimal  parameters  for  DE,  this 
facility  is  the  best  the  community  has  for  studies  that  require  good  reproducibility  and 
advanced  diagnostics.  These  results  are  in  a  reasonable  agreement  with  the  TWG  and 
MKQ  scaling  models,  and  translate  into  a  15-20  kJ/cm  prediction  of  the  optimized  K- 
shell  yield  on  DQ,  implying  over  60  kJ  total  yields,  if  the  pinch  lights  up  uniformly  over 
its  length.  (It  was  noted  by  C.  Coverdale  that  the  proportionality  between  the  yield  and 
the  pinch  length  cannot  be  taken  for  granted.  P.  Coleman  reminded  the  group  that 
uniform  K-emission  was  observed  in  some  DE  shots  with  the  1 -2-3-4  nozzle).  Topics 
selected  for  discussion  by  the  team  will  build  on  this  achievement,  to  push  the  radiative 
performance  to  the  limit  and  elucidating  the  relevant  implosion  physics. 

Preionization 


P.  Coleman  started  the  discussion  of  pre-ionization.  Study  of  the  pre-ionization  issue  is 
important  and  it  can  become  increasingly  more  relevant  as  driving  voltage  and  gas  density 
decrease  and  initial  load  radius  and  implosion  time  increase.  It  is  still  not  clear  that  pre¬ 
ionization,  as  presently  conducted,  has  a  consistent  effect  on  implosion  dynamics  and  on 
radiation  yield.  The  effects  may  be  difficult  to  measure  because  of  variations  in  other  parameters 
(e.g.,  current  and  mass  vary  by  at  least  ~5%).  [A  point  was  made,  that  it  has  not  been  shown  that 
the  radiation  fi’om  the  pinch  is  azimuthally  symmetric  -  also  adding  to  measurement 
ambiguities.] 

It  was  agreed  that  the  following  describes  the  state  of  understanding  of  the  pre¬ 
ionization: 

-  DM2  experiments  have  shown  more  consistent  K-shell  output,  when  the  pre¬ 
ionizer  was  used; 

-  Some  good  shots  were  also  obtained  without  the  pre-ionizer; 
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-  Providing  the  azimuthally  uniform  pre-ionization  can  come  into  conflict  with 
the  debris  shield  technical  constraints 

It  was  also  suggested  that  the  prepulses,  associated  with  DE  and  DQ,  may  be  more  than 
sufficient  to  preionize  the  gas:  on  DE,  a  5%  “toe”  current  of 200  kA  flows  through  the  target  ~50 
ns  before  the  main  pulse;  on  Z,  duration  of  the  7%  “toe”  of  1  MA  is  about  the  same.  Even  with  a 
low  level  prepulse  on  DQ,  1  %  of  peak  current  flows  through  the  load  for  ~10  ns.  To  explore 
fiuther  the  prepulse-induced  pre-ionization  may  require  a  bench  test,  which  can  simulate  current¬ 
time  characteristics  of  DQ.  Since  the  main  issue  here  is  the  reproducibility  of  the  positive  effect, 
if  any,  of  the  pre-ionization  on  the  K-shell  yield,  the  results  of  a  considerable  number  of 
otherwise  identical  shots  made  with  and  without  pre-ionization  need  to  be  compared.  This 
translates  into  a  requirement  of  at  least  one  test  week  on  DE,  which  should  be  used  by  the  team 
for  8  and  12  cm  nozzle  shots  with  and  without  pre-ionization,  so  as  to  establish  the  pre-ionizer’s 
effect.  Of  course,  this  study  can  and  should  be  combined  with  the  study  of  other  issues,  like 
robustness  of  the  optimal  mass/density  distribution,  etc.  It  was  agreed  that  shots  should  be 
properly  instrumented  so  that  useful  data  can  be  provided  to  modelers  -  at  this  time  it  is  not  even 
clear  what  pre-ionization  levels  are  being  obtained.  (B.  Weber  has  measured  ~1%  ionization, 
under  specific  conditions.)  R.  Terry  suggested  that  if  the  benefit  of  the  pre-ionizer  is 
inconclusive,  then  “new  technology”  could  be  attempted,  such  as  e-beams  or  lasers. 


Reviewing  the  status  of  the  2-D  modeling,  R.  Terry  highlighted  the  status  of  numerical 
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and  physical  testing  of  Mach2  code,  incorporation  of  the  NRL  radiation  p 

lysics  package  into  it. 

and  compared  the  practical  applicability  of  2-D  codes  Mach2  and  Delta  vs.  the  1-D  RMHD  code 
of  NRL  Code  6720.  The  results  of  his  very  rough  comparison  are  listed  in  the  table. 

M.  Frese  discussed  the  recent  development  of  the  Mach  2  code  and  showed  the  initial 
predictions  for  Quad  12  cm  nozzle  results.  He  qualified  the  results,  stating  that  while 
energy  balance  is  now  satisfied  in  the  code;  other  checks  have  not  been  completed. 

(Present  results  indicate  >  100  kJ  of  K-shell  radiation.)  An  important  result  of  the  Mach 
2  effort  is  ability  to  model  tiie  pinch  phase  of  the  implosion.  A  sample  calculation 
showed  the  behavior  of  several  quantities,  including  the  K-output  rising  up  after  the  main 
pulse.  S.  Chantrenne  reported  a  number  of  simulations  of  the  Ar  shots  on  DE  performed 
by  herself  and  P.  Steen  using  their  version  of  Mach2  code  and  the  simulated  gas  density 
profiles  produced  by  the  Delta  code.  Their  results  for  total  mass  scan  of  the  K-shell  yield 
were  in  good  agreement  with  P.  Coleman’s  experimental  data. 
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Optimal  mass  distribution,  origin  of  K-shell  emission,  mixing,  etc. 


Discussion  of  this  issue  started  with  P.  Coleman’s  remark  that  the  DTRA 
community  is  not  ready  to  design  and  build  a  15  - 16  cm  diameter  nozzle.  The  essential 
physics  needed  to  make  such  a  nozzle  work  as  a  good  radiator  is  not  yet  fully  understood. 
Do  we  need  a  uniform  fill?  A  gentle  gradient  peaked  near  the  axis?  Two  concentric  shells 
imploding  on  a  central  jet?  H.  Sze  noted  that  precisely  for  this  reason  we  need  to  learn  as 
much  as  possible  experimenting  with  the  12-cm  diameter  nozzle  on  a  well-diagnosed 
driver,  DE.  It  is  necessary,  in  particular,  to  observe  the  structme  of  K-  and  L-shell 
radiating  plasma,  to  reliably  measure  these  yields,  as  well  as  the  total  yield  (P.  Coleman). 
Large  radii  and  >2  cm  lengths  are  favored  for  DH.  The  presumed  benefit  of  this 
approach  is  that  the  radiation  output  from  such  a  PRS  could  contain  20  -  40  %  energy  in 
the  K  region  of  the  spectrum.  The  main  issue  is  to  determine  the  optimum  distribution  of 
mass. 

H.  Sze  suggested  that  the  team  proceed  with  optimization  of  mass  distribution  and 
diagnose  the  shell  dynamics  (compression,  mixing),  including  distributions  with  a  central 
jet.  Doping  separate  shells  with  Freon  and,  possibly,  with  H2S  should  answer  question 
with  regard  to  which  mechanisms  optimize  the  yield.  Estimated  number  of  DE  shots  is 
50.  R.  Terry  commented  that  it  is  not  clear  that  multiple  shells  produce  higher  yields, 
compared  to  solid  fills. 

The  mass  optimization  effort,  taken  together  with  measurements  of  the  ratio  of  the 
K-to-total  yield,  will  also  show  how  this  affects  the  DMS  design. 

Non-radiatiag  (“dark”)  mass 

It  is  also  necessary  to  know  the  amoimt  of  initial  mass  that  does  not  participate  in 
K-shell  radiation,  since  it  appears  to  be  a  significant  fraction  of  the  total.  This  issue  is 
related  to  current  sheath  thickness,  that  if  controlled  (and  this  may  be  very  difficult), 
could  produce  better  yields  and/or  better  K/total  ratio.  A  question  remains:  if  100  %  of 
mass  could  be  made  to  radiate  -  how  much  would  be  the  increase  in  the  yield?  A 
possible  channel  for  loss  of  mass  -  through  the  pinch  ends  has  been  partially  tested  (P. 
Coleman),  by  looking  outside  the  anode  with  a  camera  in  the  visible.  The  camera 
showed  no  light  there.  M.  Krishnan  noted  that  the  mass  left  behind  (whose  presence  is 
indicated  by  simulation  and  some  experimental  data,  mostly  from  wire  array  experiments 
on  MAGPIE  and  Z)  should  have  been  seen  on  the  interferograms,  which  has  not  been  the 
case. 

H.  Sze  proposed  2  weeks  of  DE  testing  be  devoted  to  this  problem.  A  key 
diagnostic  would  be  the  L-shell  imaging  and  spectroscopy.  Existing  atomic  data  bases  at 
NRL  would  have  to  be  enhanced  to  properly  interpret  the  data  obtained  in  such  a 
campaign. 


Low-Z  Outer  Shell 


The  reason  for  considering  this  issue  is  that  the  outer,  poorly  radiating  shell,  may 
act  as  an  initial  current  carrier,  which  gets  switched  to  the  Ar  shell,  allowing  it  to 
compress  rapidly.  The  shorter  conduction  in  Ar  may  produce  a  better  pinch,  utilizing  the 
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total  mass.  Previous  Titan  results,  using  low  Z  outer  shell  have,  have  shown,  that  in  200 
ns  implosions,  the  yield  decreases,  as  the  mass  of  the  H2  shell  increases  to  50  %  of  the 
total  mass.  The  same  was  seen  by  Shishlov  in  Tomsk.  He  qualified  his  results  with  a 
possibility  that  it  was  not  clear  what  the  initial  conditions  were;  possibly,  the  fast  flow  put 
hydrogen  gas  in  places  where  it  diverted  the  current  fi'om  the  implosion. 

To  follow  up  this  line  of  thought,  it  was  suggested  that  this  approach  should  be 
tried  with  higher  Z  gasses  (CH,  CO2,  N2, 02) 

Enhanced  energy  coupling. 

Experiments  have  shown  that  the  total  radiation  from  the  pinch,  in  many  cases, 
exceeds  the  JxB  energy  derived  from  the  magnetic  field  surrounding  the  implosion. 
Additional  energy  can  be  derived  fi’om  the  magnetic  energy  stored  in  external 
inductances,  if  the  pinch  has  a  resistive  component.  Such  effects  were  seen  in  Z 
experiments  and  in  Shishlov’s  low  current  (400  kA)  implosions,  where  the  resistive 
component  "was  measured.  It  is  not  clear,  whether  the  increased  coupling  could  lead  to 
increased  emission  of  the  K-shell  lines. 

This  problem  will  need  to  be  examined,  because  the  increase  in  the  total  radiation 
would  add  to  significant  additional  stress  of  the  DMS,  or  alternately,  it  may  increase  the 
K-emission,  producing  higher  fluence  at  the  test  article.  If  the  result  of  enhanced 
dissipation  is  only  UV  emission,  then  we  might  want  to  minimize  it,  in  order  to  mitigate 
the  issue  of  debris  shielding.  A  low-cost  approach  to  addressing  this  problem  (A. 
Velikovich)  would  be  to  improve  the  accuracy  of  voltage  and  inductance  measurements 
on  DE,  in  the  same  way  as  implemented  by  E.  Waisman  on  Z.  The  effective  enhanced 
resistance  should  be  measured  and  found  consistent  with  the  approximate  radial 
trajectory  of  the  implosion,  r(t).  These  measurements  could  be  analyzed  as  piggyback 
experiments  on  any  DE  shots. 


Diagnostics 


The  preferred  diagnostic  for  recording  the  spectra  of  tracer  element  lines  is  the  Johann 
spectrometer  (E.  Yadlowsky).  Detection  of  the  “dark”  or  “lost”  mass  would  require 
interferometry  and  L-shell  spectroscopy.  It  would  be  most  convenient  to  use  radially  resolved 
spectroscopy  with  McDonald’s  spectrographs,  like  the  5-channel  spectrograph  used  at  Sandia. 
Bare  bolometers  need  to  be  fielded  to  measure  the  total  radiation  yield  (P.  Coleman).  Since  the 
voltage  measurements  are  much  less  reliable  on  DE  than  on  Z,  this  need  to  be  improved,  at  least 
by  adding  shielded  B-dot  probes  (M.  Krishnan). 
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Prioritization 


At  the  end  of  the  meeting,  the  above  topics  were  prioritized  by  voting  by  selecting  two  subjects 
most  important  for  the  conummity  to  address.  Each  participant  voted  for  t^  topics  of  greatest 
importance. 


Topic 

Number  of  votes 

Optimizing  mass  ratio 

12 

K-emitting  shell 

9 

Non  radiating  matter 

8 

Low-Z  outer  shell 

2 

Preionization 

1 

Attendees 

NRL  6720:  J.  Apruzese,  R.  Terry,  A.  Velikovich 


TPSD:  H.  Sze,  B.  Failor,  J.  Levine,  S.  Chantrenne,  P.  Steen 
AASC:  M.  Krishnan,  P.  Coleman,  E.  Waisman 


Numerex:  M.  Frese,  S.  Frese 


Sandia:  C.  Coverdale 
Hytech:.  E  Yadlowsky 


Ktech:  D.  Lepell 
NGIT:  I.  Vitkovitsky 
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II.  Zero-  and  One-Dimensional  modeling  of  DQ  Experiments 


A.  Introduction 

In  this  section  the  results  of  our  ID  modeling  of  Decade  Quad  and  Double-Eagle  12  cm 
diameter  nozzle  experiments  are  described.  This  section  begins  with  a  discussion  of  OD 
and  ID  analysis  that  was  provided  to  the  PRS  community  before  the  2003  DQ  shots  were 
performed.  The  OD/ID  analyses  are  then  updated  by  taking  into  account  a  more  realistic 
assessment  of  the  current  that  was  delivered  to  the  load  in  the  DQ  experiments. 

The  one-dimensional  (ID)  models  are  not  capable  of  modeling  all  the  effects  present  in 
multi-dimensional  Z-pinch  plasmas.  However,  by  a  judicious  choice  of  transport 
coefficients  (viscosity,  heat  conductivity,  and  ion-electron  heat  exchange),  a  reasonable 
match  between  calculated  and  measured  K-shell  yields  can  be  made  over  a  broad  range  of 
experimental  conditions  for  a  given  nozzle  configuration.  We  cannot  fully  quantify  the 
extent  that  our  phenomenological  model  captures  the  physics  of  actual  multi-dimensional 
plasma  implosions,  so  to  some  extent  any  insights  gained  by  this  work  are  inherently 
speculative.  However,  for  now  (until  Mach2  is  a  fully  operational  code  containing 
adequate  radiation  physics),  this  is  likely  the  best  modeling  that  can  be  accomplished  that 
includes  the  essential  radiation  physics  of  the  implosion.  For  a  more  complete  description 
of  the  ID  modeling  see  the  2002  Annual  report. 

B.  OD  and  ID  K-yield  assessment  prior  to  2003  DQ  experiments 

The  scaling  of  K-shell  radiation  with  current  and  imploded  mass  can  crudely  be  modeled 
using  OD  snowplow  scaling  models.  These  models  reasonably  model  the  implosion  time, 
peak  current  into  the  load,  and  the  specific  energy  of  the  load.  This  specific  energy  is 
often  represented  by  the  T]*  parameter,  which  is  the  specific  energy  divided  by  the 
minimum  energy  per-unit-mass  needed  to  instantaneously  ionize  the  plasma  into  the  K- 
shell.  However  the  scaling  of  the  K-shell  radiation  with  current,  which  depends  so  much 
on  the  individual  load  and  machine  design,  has  historically  been  very  speculative.  Figure 
1  shows  OD  model  predictions  for  the  2001  argon  experiments  on  Z. 

There  are  two  K-shell  yield  scaling  curves  shown  in  Fig.  1.  Curve  A  is  generated  using 
the  enclosed  “old”  Z  scaling  model  (Fig.  2),  it  is  a  typical  OD  model,  which  has 
historically  been  used  to  predict  K-shell  yield.  This  model  shows  that  a  2  mg/cm  load  is 
capable  of  producing  about  125  kJ/cm  of  K-shell  emission.  Unfortunately,  this  model 
does  not  adequately  represent  the  rapid  fall-off  in  K-shell  yield  with  mass  that  was 
exhibited  in  the  Z  experiments  using  the  4-3-2-1  nozzle.  This  fall  off,  for  which  strong 
evidence  was  also  exhibited  in  recent  DQ  and  long  pulse  Saturn  4-3-2- 1  nozzle 
experiments,  is  likely  to  limit  the  mass  that  can  be  imploded  by  DQ  and  DH  to 
efficiently  produce  K-shell  emission.  The  effect  may  be  even  more  pronounced  for  larger 
radius  loads.  There  is  speculation  that  this  fall  off  is  inherently  an  instability  effect  that 
occurs  in  the  larger  radius  implosions.  This  effect  has  also  been  discussed  in  our  Final 
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reports  to  DTRA  since  1998.  These  reports  attribute  the  fall  off  to  L-shell  radiative 
cooling  rates  being  larger  than  currently  modeled  or  possibly  to  substantial  mixing  which 
enhances  energy  exchange  between  the  hot  core  and  the  rest  of  the  plasma.  Thus,  the 
ionization  levels  needed  to  support  K-shell  emission  would  diminish  faster  with 
increasing  mass  than  presently  modeled.  This  effect  is  crudely  taken  into  account  by 
assuming  that  a  necessary  requirement  for  two  different  loads  to  produce  K-shell 
radiation  at  the  same  efficiency  is  that  they  must  have  approximately  the  same  ionization 
level,  which  is  largely  determined  by  the  ratio  of  plasma  heating  to  cooling.  Taking  the 
heating  rate  at  stagnation  as  being  proportional  to  1/2  mv^/t  and  the  radiative  cooling  rate 
as  proportional  to  massVx,  where  x  is  a  thermalization  time,  this  ratio  can  be  simplified  to 
rj*/mass.  Scaling  curve  B  is  attained  by  reducing  the  K-shell  yield  of  curve  A  by 
(q^/mass)/. 00594,  where  .00594  is  the  ratio  of  q*  (1.7)  to  mass  (286  pg/cm)  that  was 
achieved  in  the  best  4-3-2- 1  nozzle  DQ  argon  experiments  to  date  (shot  #  498  -  28  kJ  K- 
shell  yield),  see  Fig.  3.  Since  this  adjustment  to  the  original  scaling  model  does  a 
reasonable  job  of  replicating  the  observed  4-3-2-1  nozzle  K-shell  yields  from  DE,  DQ, 
and  Z  it  is  a  reasonable  scaling  model  for  interpolating  the  performance  of  4-3-2-1  nozzle 
loads  on  DH.  Likewise  the  model  can  be  used  for  extrapolating  yield  predictions  to 
larger  radius  DQ  and  DH  loads.  However  because  it  is  an  extrapolation  into  an  unknown 
larger-radius  parameter  space,  one  has  to  be  skeptical  about  any  prognostications  derived 
from  this  model.  An  example  of  this  OD  extrapolation  to  DQ  with  the  6-4-2-1  nozzle  is 
shown  if  Fig.  4. 

Preliminary  ID  results  for  DQ  are  also  shown  in  Figure  4.  These  results  are  based  on  a 
consistent  set  of  transport  coefficients  that  not  only  reasonably  replicates  the  existing  4-3- 
2-1  nozzle  results  for  DQ  and  Z  but  also  reasonably  replicates,  or  so  we  thought  at  the 
time  the  calculations  were  made,  the  6-4-2-1  nozzle  DE  experiments.  This  previous 
statement  is  qualified  because  at  the  time  the  calculations  were  made  we  thought  that  the 
inner  to  outer  mass  ratio  was  1:1,  not  1.8:,  as  used  in  the  DE  experiments.  Also,  the  John 
Riordan  version  of  the  open  circuit  voltage  for  long  pulse  Double-Eagle  was  used,  see 
Fig  5.  This  voltage  profile  is  not  accurate  for  the  longer  times  needed  to  model  the  12  cm 
diameter  nozzle  experiments.  To  reasonably  match  this  larger  parameter  space,  the 
enhancement  on  viscosity  had  to  be  increased  to  200.  This  is  a  factor  of  4  larger  than  was 
required  to  match  only  the  4-3-2-1  nozzle  experiments.  This  factor  of  4  produced  little 
effect  on  4-3-2-1  nozzle  calculations  but  it  dramatically  altered  the  12  diameter  results. 
The  mistakes  made  in  modeling  the  DE  experiments  were  discovered  after  the  DQ  12  cm 
diameter  experiments  were  performed.  When  the  more  realistic  300  ns  DE  open  circuit 
voltage  profile,  see  Fig.  5,  and  proper  mass  ratio  are  employed,  it  was  found  that  an 
enhancement  of  250  on  viscosity  (factor  of  5  over  that  used  to  model  just  4-3-2-1) 
experiments  was  required  to  best  replicate  the  12  cm  diameter  DE  experiments  as  well  as 
the  4-3-2- 1  DQ  experiments. 

Comparing  the  ID  to  the  OD  predicted  yields  shown  in  Figure  4  reveals  that  the  ID 
models  predict  substantially  less  K-yield  than  either  our  new  OD  or  old  OD  models  for  the 
larger  mass  12  cm  diameter  loads.  This  is  a  reflection  on  the  fact  that  the  ID  models 
show  that  not  only  are  mass  and  energy  important  in  determining  K-yields  as  do  the  OD 
models,  but  also  the  stagnation  physics;  some  of  which  is  presumably  captured  by  the 
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phenomenological  ID  model.  The  ID  optimal  yield  is  shown  to  occur  at  smaller  mass 
load  than  that  of  the  OD  predictions.  All  three  models  show  similar  yield  behavior  for 
mass  loads  less  than  250  pg/cm.  The  experiments  can  be  assessed  in  terms  of  a  OD  or  ID 
model  perspective  by  comparing  these  results  to  the  experimental  data. 

C.  Experimental  Results 

Table  1  summarizes  the  K-shell  yield  results  from  the  12  cm  diameter  nozzle  DQ  argon 
experiments  as  well  as  a  few  selected  4-3-2-1  DQ  results  and  DE  12  cm  diameter  results. 
Note,  the  efficiency  listed  in  this  table  should  be  interpreted  as  a  relative  measure  of  the 
efficiency  of  converting  the  energy  that  enters  the  diode  into  K-shell  emission.  Since  the 
amount  of  compressional  work  done  on  each  plasma  load  is  not  really  known,  the  energy 
coupled  to  the  load  assuming  a  10:1  compression  is  used  as  a  relative  figure  of  merit  of 
the  electrical  energy  that  enters  the  diode.  In  other  words,  we  really  don’t  know  the 
absolute  efficiency,  but  if  one  load’s  quoted  efficiency  is  10%  and  another  load’s  is  20%, 
then  one  can  say  that  the  20%  conversion  efficiency  load  was  twice  as  efficient  as  the 
10%  load.  Note,  the  efficiency  of  the  12  cm  loads  tended  to  increase  with  the  specific 
energy  as  represented  by  q*  in  Table  1.  The  highest  efficiency  attained  for  the  12  cm 
diameter  nozzle  (13.2  percent)  was  significantly  below  the  best  4-3-2-1  yield  (19.5%).  A 
reasonable  goal  for  the  12  cm  diameter  loads  on  DE  and  DQ  would  be  to  achieve  the 
same  efficiencies  that  were  achieved  using  the  4-3-2-1  nozzle. 

D.  Normalized  OD  and  ID  assessments  of  DQ  12  cm  diameter 
nozzle  experiments 

The  current  that  was  coupled  to  the  load  in  the  DQ  experiments  was  found  to  be 
significantly  less  than  what  was  predicted  by  the  equivalent  circuit  model  for  DQ,  see 
Fig.  5.  In  order  to  analyze  the  results  of  these  experiments  from  a  OD  and  ID  perspective 
it  is  necessary  to  recalculate  the  yields  using  a  circuit  model  that  more  accurately  matches 
the  peak  current  delivered  to  the  load.  Little  sensitivity  was  found  in  terms  of  K-yield  to 
whether  the  6-4-2-1  nozzle  or  the  6-5-3-2  nozzle  were  modeled.  The  results  presented 
here  are  for  the  6-5-3-2  nozzle  with  an  initial  mass  distribution  sirnilar  to  that  measured 
near  the  cathode  plane.  The  ratio  of  inner  to  outer  mass  was  2.2:1  and  the  viscosity  is 
now  enhanced  by  a  factor  of  250.  The  DQ  open  circuit  voltage  shown  in  Fig.  5  was 
multiplied  by  the  factor  0.94  in  order  to  reproduce  the  peak  currents  in  the  OD/ID  models 
that  were  achieved  in  the  experiments,  see  Fig.  6.  These  peak  currents  are  given  as  a 
function  of  load  mass.  The  K-shell  yields  as  calculated  by  the  OD/ID  models  with  peak 
currents  given  by  Fig.  6,  are  compared  with  experiment  in  Fig.  7.  Note,  the  implosion 
times  (not  shown)  were  in  good  agreement  with  the  experiments  but  they  are  sensitive  to 
the  choice  of  initial  mass  distribution.  For  the  purposes  of  these  calculations,  which 
focussed  on  insuring  that  the  same  amount  of  energy  was  available  to  the  theoretical 
loads  as  the  experimental  loads,  no  initial  mass  was  modeled  beyond  6  cm  radius. 

Based  on  a  comparison  of  experimental  and  ID  K-shell  yields,  see  Fig.  7,  it  appears  that 
a  reasonable  job  of  phenomenologically  modeling  some  of  the  large-radius  stagnation 
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physics  with  our  ID  model  was  accomplished.  While  the  ID  results  are  ostensibly  only 
for  the  6-5-3-2  nozzle  with  a  2.2  :  1  inner-to-outer  mass  ratio  .there  is  not  much  model 
sensitivity  between  the  6-5-3-2  and  6-4-2-1  nozzles.  Actually  the  ID  model  predicts 
slightly  more  yield  for  shots  547  and  547,  which  employed  the  6-4-2- 1  nozzle,  than 
shown  in  Fig.  7.  This  is  primarily  because  they  were  shorter  length  pinches,  3.7  cm  as 
opposed  to  3.9  cm  for  the  Titan  experiments.  The  ID  and  experimental  yields  and  trends 
(increasing  yield  with  decreasing  mass)  are  in  excellent  agreement  for  mass  loads  greater 
than  300  pg/cm.  The  only  significant  deviation  between  the  model  and  the  experiment 
was  for  shot  #548,  which  was  the  lowest  mass  load  (273  pg/cm).  Given  that  the  ID 
model  shows  optimal  yields  for  load  mass  less  than  350  pg/cm  and  that  there  is  a  paucity 
of  experimental  data  in  this  mass  region  (only  one  shot)  -  we  should  be  encouraged  to 
thoroughly  explore  this  region  experimentally. 

E.  Tradeoff  between  stability  and  energy  coupling 

The  2004  DE  and  DQ  large  radius  experiments  are  going  to  systematically  explore  the 
sensitivity  of  K-shell  yield  to  the  initial  mass  distribution.  Here  we  examine  this  issue 
with  the  ID  model  that  has  been  benchmarked  to  agree  with  the  2.2:  1  mass  ratio  Double- 
Eagle  experiments.  In  terms  of  a  ID  model  perspective,  it  is  generally  true  that  the  more 
energy  coupled  to  the  load  -  the  more  K-shell  emission  will  be  calculated.  This  is 
because  the  ID  models  cannot  account  for  effects  due  to  instabilities.  For  this  reason  the 
ID  models  will  predict  that  the  distribution  with  the  smallest  mass  ratio  will  produce  the 
most  K-shell  emission.  In  this  section  we  display  the  calculated  yields  for  Double  Eagle 
and  Decade  Quad  for  the  6-5-3-2  nozzle  in  which  the  mass  distribution  is  varied  from 
2.2:  1  inner-to-outer  ratio  to  a  ratio  of  1  :  2.  Presumably  the  2.2  :  1  distribution  is  more 
stable  than  the  1  :  2  distribution  because  more  mass  is  located  close  to  the  axis  of  the 
pinch.  These  calculated  yields  as  a  function  of  load  mass  and  distribution  are  shown  in 
Figs.  8  and  9  for  DQ  and  DE,  respectively.  Insights  into  this  tradeoff  between  stability 
and  energy  coupling  will  be  obtained  by  comparing  these  results  with  the  upcoming 
experimental  results. 


Summary 

Once  benchmarked  to  the  4-3-2- 1  nozzle  experiments  performed  on  DE,  DQ  and  Z  the 
OD  model  bases  all  of  its  K-shell  yield  predictions  on  the  load  mass  and  energy  coupled 
to  the  load  as  calculated  by  a  10:1  compression  ratio.  This  model  can  be  used  to 
reasonably  interpolate  yield  behavior  for  the  4-3-2-1  nozzle  on  Z,  DE,  DQ  and  DH. 
When  the  OD  model  is  extrapolated  to  predict  the  performance  of  the  6-5-3-2  nozzle  shot 
#5161  experiment  on  DE,  see  Fig.  9,  it  also  does  a  reasonable  job.  However,  when  the  OD 
model  is  extrapolated  to  predict  the  performance  of  DQ  12  cm  diameter  loads.  It 
substantially  over  predicts  the  yields  for  every  experimental  mass  load  tested,  see  Fig.  7. 
From  a  OD  perspective  one  might  be  tempted  to  speculate  that  because  the  OD  model 
systematically  over  predicts  the  DQ  yields,  but  not  DE’s,  that  the  DQ  load  performance 
was  somehow  degraded  below  that  of  Double-Eagle. 
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The  ID  model  uses  a  different  phenomenological  approach  than  the  OD  model  in  order  to 
match  the  experimental  data  base.  Rather  than  fixing  a  compression  ratio  or  ratio  of 
thermalization  to  radiation  rates  as  is  done  in  the  OD  modeling,  the  ID  model  uses  a 
judicious  choice  of  transport  coefficients  to  benchmark  to  previous  experiments.  Thus 
the  ID  model  bases  its  results  on  a  self-consistent  coupling  of  load  mass,  energy 
coupling,  and  implosion/radiation  dynamics.  Here  the  ID  model  was  benchmarked  to  the 
data  base  of  4-3-2- 1  experiments  as  well  as  the  shot  5161  DE  6-5-3-2  experiment. 

Similar  to  the  OD  model,  the  ID  model  can  be  used  to  reasonably  interpolate  yield 
behavior  for  the  4-3-2-1  on  Z,  DE,  DQ,  and  DH.  By  construction  the  ID  model  does  a 
reasonable  job  of  representing  DE  shot  Slbl’s  yield  behavior.  However,  unlike  the  OD 
model,  the  ID  model  also  reasonably  predicts,  with  one  exception,  the  yield  behavior  of 
the  DQ  experiments  (Fig.  7).  The  one  exception  was  DQ  shot  #548,  which  was  the  lowest 
mass  load  (273  p,g/cm)  experiment.  Thus,  from  a  ID  perspective  one  can  say,  in  direct 
contrast  to  the  OD  perspective,  that  the  12  cm  diameter  loads  performed  nearly  the  same 
on  DQ  as  DE. 

The  ID  model  predicts  that  the  highest  DQ  K-yields,  based  on  past  drive  currents,  will 
occur  for  load  masses  less  than  350  |ig/cm  and  that  yields  in  excess  of  9  kJ/cm  are 
attainable.  Unfortunately,  the  yield  prediction  for  shot  #548  was  10  kJ/cm  and  only  7.2 
kJ/cm  was  achieved.  Given  the  power  flow  problems,  and  that  this  was  the  only  shot 
taken  in  the  range  predicted  to  give  the  best  yields,  we  should  not  be  overly  discouraged 
by  this  result.  On  the  other  hand  it  is  discouraging  that  it  required  5.5  MA  to  achieve  28 
kJ  yield  with  the  12  cm  diameter  nozzle  while  only  4.6  MA  was  needed  to  achieve  the 
same  yield  with  the  8  cm  diameter  nozzle  (approximately  same  mass).  Until  it  is 
experimentally  demonstrated  that  the  deleterious  effects  on  K-yield  due  to  increased 
radius,  such  as  instabilities,  can  be  overcome  by  the  larger  drive  currents  attainable  at 
large  radius,  there  is  no  compelling  reason  to  pursue  >  12  cm  diameter  loads. 

The  2004  DQ  experiments  will  address  this  load  design  issue  by  mapping  out  the  tradeoff 
between  stability  (higher  mass  ratio)  and  more  current/energy  coupling  (lower  mass  ratio) 
to  search  for  a  parameter  regime  in  which  substantially  more  K-shell  yield  is  produced  at 
12  cm  diameter  than  8  cm.  The  ID  model  is  little  help  in  this  regard  in  that  using  the  ID 
model  to  make  yield  predictions  for  mass  distributions  other  than  the  benchmarked  2.2:1 
distribution  is  suspect  (Fig.  8).  For  this  task  our  2D  model,  which  is  currently  being 
developed,  is  better  suited. 
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Fig.  1  Calculated  and  exerimenfal  k-shell  yields  as  a  function  of 
mass  for  argon  loads  on  the  Z  machine.  Experimental  yields  are 
from  2001  experiments  with  the  4-3-2- 1  nozzle. 


Z  Scaling  of  K-shell  Emission 

(old  scaling  -  taken  from  Dense  Z-Pinches  4th  Intemational  Conf.  pp.  197) 
The  origin^  published  scaling  was  slightly  altered  to  do  a  better  job 
of  scaling  in  the  low  t)*  reginie.  The  fonn  currently  used  is  given  below: 

Energy  (j/cm)  =  (kinetic  energy  +  internal  energy)  (j/cm) 

n*  =  Energy  •  I.l  -IC^ 

Mass  (g/ctn)  • 

Z  is  the  atomic  nmnber 

o  =  2J6  •  10~^^  »  maximum  (1.0,  t)*^/(n*+12)) 
exp(-20.6/Z-9) 

P  =  minimum  of  0.25  •  Energy 
and 

0.3  -Energy*  Mass(g/cm) 

— 

K-Yield  (j/cm)  =  P  forn*>  1.75 

K-Yield  (j/cm)  =  p  •  (t)*  -  0.75)  for  1.0  <  t)*  <  1.75 

K-Yield  (j/cm)  =  0.0  f6ri)*<  1.0 

Nme,  a  leaso^le  approximation  for  the  energy  in  (J/cm)  is  given  by 
Ip2»1800.  which  is  close  to  the  coupled  energy  (10:1)  compression 
for  most  machines.  The  peak  current  Jp  ,  is  in  units  of  MAs. 

Fig.  2 
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New  Z  -scaling  -  phenomenologically 
matched  to  4-3-2- 1  nozzle  results 

Replace  old  p  by  p  x  r|* 

Mass  (pg/cm)  x  .00594 

and  recalculate  K-yield 


The  maximum  efficiency  is  limited  to  18  percent  at  low 
and  it  is  allowed  to  rise  to  27  percent  at  high 


K-yield  =  min(K-yield,  0.18  x  energy  x  xcon) 


xcon  =  min(  1.5,  (1.+  .6  x  r|*/5.0)  ) 


Fig.  4  OD  and  ID  calculated  K-shell  yields  for  DQ  with 
a  12  cm  diameter  nozzle  and  assumed  inner  to  outer 
mass  ratio  of  1 : 1  .  These  calculations  were  performed 
prior  to  the  2003  DQ  experiments. 


14 


Open  Circuit  Voltage  (MV) 


o 

D 


2 


1 


0  100  200  300  400 

Time  (ns) 


Fig.  5  Double-Eagle  and  Decade  Quad  open  circuit  voltage 
profiles  and  circuit  parameters. 
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Current  (MA) 


Fig.  6  Calculated  and  experimental  peak  currents  as  a  function 
of  load  mass  for  6-5-3-2  nozzle  on  DQ.  The  modeled  and 
experimental  inner-to-outer  mass  ratio  is  2.2  :  1. 


16 


Argon  K-Yield  (kJ/cm) 
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Fig.  7  OD  and  ID  calculated  K-yields  as  a  function  of  mass 
for  6-5-3-2  nozzle  with  2.2  :  1  inner-to-outer  mass  ratio.  The 
DQ  12  cm  nozzle  results  are  also  shown  for  comparison.  Only 
shots  548,  549  and  544  had  a  2.2  :  1  mass  ratio. 
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Argon  K- Yield  (kJ/cm) 


Fig.  8  ID  calculated  K-yields  for  6-5-3-2  nozzle  as  a  function 
of  mass  and  inner-to-outer  nozzle  mass  ratio  on  Decade  Quad. 
Experimental  results  are  also  shown.  Shots  548,  549  and  544 
were  6-5-3-2  experiments  with  2.2:1  mass  ratio.  Shots  546  and 
547  were  6-4-2-1  experiments  with  1.4:1  and  1.8:1  mass  ratios, 
respectively. 


100  150  200 


Mass  (pig/cm) 


Fig.  9  ID  and  OD  calculated  K-shell  yields  as  a  function  of  mass 
load  and  mass  distribution.  The  6-5-3-2  nozzle  is  modeled  in  these 
calculations. 
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Table  1 


Shot 

Mass 

(^ig/cm) 

length 

(cm) 

1  peak 

(MA) 

K-yield 

(kj) 

T]* 

8 

(%) 

498  (tpsd-8) 

286 

3.8 

4.64 

28.7 

1.45 

19.5 

543  (tpsd-8) 

443 

4.0 

5.06 

17.4 

1.11 

9.4 

544  (tpsd-12) 

348 

3.9 

5.65 

26.1 

1.76 

11.6 

545  {tpsd-12) 

348 

3.3 

5.78 

26.2 

1.84 

13.2 

546  (aasc-12) 

425 

3.7 

6.14 

18.9 

1.70 

7.5 

547  (aasc-12) 

450 

3.7 

6.15 

21.5 

1.61 

8.5 

548  (tpsd-12) 

273 

3.9 

5.53 

28.1 

2.15 

13.1 

549  (tpsd-12) 

348 

3.9 

5.83 

28.1 

1.88 

11.8 

5161(DE-12) 

185 

3.8 

3.47 

9.0 

1.25 

10.9 

4428(DE-4321) 

113 

4.0 

3.7 

12.0 

2.33 

12.5 

5134(DE-4321) 

3.8 

3.8 

18.0 

18.2 
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III.  Enhanced  energy  coupling  and  x-ray  emission  in  Z-pinch 

plasma  implosions 

I.  INTRODUCTION 

The  short-circuit  current  risetime  of  the  Saturn  generator  at  Sandia  National  Labora¬ 
tories  was  recently  increased  from  ~  60  to  ~  180  ns  in  order  to  investigate  the  effect  of 
current  risetime  on  the  implosion  dynamics  of  Z-pinches.  The  effect  was  large.  In  both 
aluminum  and  tungsten  wire  experiments,  significant  increases  in  total  x-ray  output  were 
observed  [1-3]  compared  to  the  earlier  outputs  that  were  seen  in  short  current-risetime  Sat¬ 
urn  experiments.  Present  magnetohydrodynamic  (MHD)  calculations  are  unable  to  generate 
such  yields  because  the  amount  of  x-ray  energy  radiated  in  the  experiments  is  significantly 
larger  than  the  sum  of  the  j  x  B  work  and  Ohmic  heating  calculated  using  the  Saturn  circuit 
model.  Such  energy  discrepancies  between  Z-pinch  calculations  and  experiment  have  existed 
for  some  time  [4-7],  but  not  to  the  degree  seen  in  these  experiments.  The  energy  inputs 
calculated  for  these  experiments  in  a  one-dimensional  (1-D)  MHD  are  in  deficit  by  factors 
of  2  to  4.  Such  energy  discrepancies  pose  serious  theoretical  challenges  even  as  they  offer 
new  possibilities  for  increasing  x-ray  extraction  efficiencies. 

In  early  work  on  imploding  Z  pinches,  it  was  observed  that  the  total  radiation  yield 
could  exceed  the  0-D  or  1-D  estimates  of  radial  kinetic  energy  by  roughly  50%.  Giuliani, 
et  al.  [8]  suggested  that  enhanced  or  anomalous  resistivity  could  explain  this  discrepancy. 
Peterson,  et  al.  [9]  subsequently  showed  that  2-D  MHD  calculations  could  increase  the 
coupled  energy  by  more  than  a  factor  of  two  over  the  0-D  and  1-D  kinetic  energy  couplings 
by  accounting  for  mass  distributions  in  the  radial  and  axial  directions  that  lead  to  reduced 
peak  radial  kinetic  energy  but  significantly  more  pdV  work.  Initial  2-D  calculations  of  the 
long  current-pulse  Saturn  experiments  [10],  however,  indicate  that  they  cannot  reproduce 
the  high  energy  coupling  and  hence  the  x-ray  output  seen  in  these  experiments.  Further 
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evidence  that  significant  energy  coupling  occurs  through  other  than  radial  kinetic  energy 
generation  and  pdV  work  comes  from  a  single  wire  experiment  described  in  Ref.  ( [11]).  The 
total  energy  radiated  in  this  experiment  was  found  to  be  comparable  to  that  measured  in  a 
set  of  imploding  array  experiments,  even  though  the  recorded  pinch  sizes  were  comparable 
in  all  cases. 

A  theoretical  mechanism  was  recently  proposed  in  which  energy  coupling  to  a  Z-pinch  is 
increased  through  the  generation  and  storage  of  additional  amounts  of  magnetic  field  energy 
in  the  Z-pinch  plasma.  The  dissipation  of  this  energy  then  increases  the  energy  coupling  to 
the  plasma  [6,7].  It  was  hypothesized  that  tubes  of  magnetic  flux  build-up  within  a  pinch 
as  an  additional  effect  of  Rayleigh-Taylor  or  m=0  instability  growth  at  the  surface  of  the 
plasma.  If  tubes  of  magnetic  flux  were  to  form  and  be  trapped  within  the  plasma,  the 
currents  that  produce  them  would  decouple  from  the  circuit  current  and  effectively  increase 
the  energy  coupling  between  the  Z-pinch  and  the  generator  transmission  line.  In  this  paper, 
however,  a  different  approach  to  enhanced  energy  coupling  is  investigated.  First,  the  problem 
of  coupling  transmission  lines  to  Z-pinch  plasmas  is  analyzed.  In  a  1-D  hydrodynamics 
calculation,  this  analysis  relates  the  strength  of  the  energy  coupling  to  the  plasma’s  surface 
resistivity.  Hence,  enhanced  coupling  occurs  because  of  increased  surface  resistivity. 

A  mechanism  for  increasing  the  electron  collision  frequency  at  the  surface  of  a  pinch  was 
proposed  in  Ref.  [12];  namely,  that  currents  flowing  in  the  low  density  plasma  at  the  pinch 
surface  acquire  drift  velocities  that  can  greatly  exceed  the  sound  speed.  This  event  leads  to 
the  build-up  of  micro-instabilities  that  increase  the  plasma’s  resistance  to  the  current  flow 
and  that  inhibit  further  growth  in  the  drift  velocity.  How  this  mechanism  operates  in  detail 
has  yet  to  be  worked  out.  Consequently,  the  approach  adopted  in  this  paper  was  to  increase 
energy  coupling  by  phenomenologically  increasing  the  strength  of  the  surface  resistivity.  The 
kind  of  increases  that  could  be  applied  were  guided  entirely  by  the  x-ray  data  obtained  in 
the  Saturn  experiments. 

This  data  was  obtained  from  three  identically  conflgured  aluminum  wire  experiments 
and  is  presented  in  section  II.  While  the  three  aluminum  shots  were  ostensibly  identical,  the 
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x-ray  data  they  produced  was  not.  Large  shot-to-shot  variations  were  seen  indicating  (1) 
possible  small  differences  in  the  wire  loads  or  in  the  Saturn  machine  operation,  or  (2)  a  high 
level  of  randomness  or  sensitivity  to  the  wire  breakdown  conditions  in  the  experiments,  or 
(3)  a  high  level  of  sensitivity  to  the  onset  of  the  resistivity  changes.  In  one  shot,  the  total 
radiated  yield  was  a  factor  of  two  larger  than  in  the  other  two  shots,  and,  in  all  three  shots, 
significantly  more  energy  was  radiated  than  could  be  accounted  for  by  conventional  MHD 
resistivity  theory. 

In  section  III,  the  problem  of  coupling  circuit  equations  to  fluid  equations  is  investigated. 
Under  the  assumption  that  the  fluid  dynamics  is  one-dimensional,  we  show  that  an  analytic 
expression  for  the  circuit  resistance  and  inductance  is  derivable  in  terms  of  surface  fluid 
quantities.  The  resistance  is  found  to  depend  on  both  the  surface  resistivity  and  the  surface 
current  density  of  the  pinch.  Because  increases  in  resistivity  increase  the  rate  of  current 
diffusion,  the  current  density  at  the  surface,  on  which  the  energy  coupling  also  depends,  is 
lowered.  Thus,  the  energy  coupling  problem  is  nonlinear,  and  MHD  calculations  are  needed 
to  quantitatively  assess  how  increases  in  circuit  resistance  actually  affect  the  fluid  dynamics. 
These  calculations  are  described  in  section  IV.  Two  of  the  problems  investigated  in  sections 
IV  concern  (1)  the  time  history  of  circuit  resistance  needed  to  approximately  model  the 
Saturn  experiments  and  (2)  the  sensitivity  of  the  calculations  to  variations  of  this  time 
history. 

The  Saturn  experiments  manifest  enhanced  energy  coupling  because  they  produced  en¬ 
hanced  x-ray  outputs.  Thus,  fluid  calculations  must  also  demonstrate  this  capability.  Some 
of  the  modeling  issues  that  this  energy  channeling  raises  are  investigated  in  section  IV.  It  is 
found  that,  under  certain  dynamical  assumptions,  the  calculations  can  convert  substantially 
all  of  the  input  energy  into  x-rays,  but  not  otherwise.  Moreover,  the  x-ray  pulses  observed 
in  the  Saturn  experiments,  both  as  to  their  energy  content  and  the  presence  or  absence  of 
sizable  amounts  of  emission  following  the  main  x-ray  pulse,  are  found  to  correlate  well  with 
the  assumed  behavior  of  the  circuit  resistance.  Thus,  the  calculations  serve  to  identify  a 
variety  of  important  issues  related  to  the  modeling  of  Z-pinch  dynamics. 
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Finally,  a  formula  for  enhanced  resistivity  produced  by  the  growth  and  saturation  of 
a  hybrid  drift-wave  instability  [12]  is  tested  in  section  V  by  post-processing  the  principal 
calculation  described  in  section  IV.  In  this  one  case,  the  resistivity  increases  predicted  by 
the  formula  do  not  support  the  resistivity  assumptions  that  were  made  in  the  calculations. 
By  a  significant  amount,  they  are  neither  large  enough  nor  do  they  occur  early  enough  to 
offer  support  to  the  resistivity  assumptions  that  were  built  into  the  calculations.  This  work 
is  then  summarized  in  section  VI. 

II.  EXPERIMENTAL  DATA 

Saturn  generates  sufficient  current  to  drive  large  wire-number  aluminum  loads  with  small 
interwire  gap  spacings  from  radii  of  two  or  more  centimeters  when  it  is  operated  in  a  long 
current-risetime  mode.  Small  gap  spacings,  in  principle,  allow  plasma  shells  to  form  prior  to 
wire  implosion,  and  these  shells  implode  with  greater  azimuthal  symmetry  than  when  gaps 
remain  present  during  the  run-in.  For  the  aluminum  wire  shots  performed  on  Saturn  in  this 
long  current-risetime  mode  of  operation,  a  gap  spacing  that  produced  the  most  symmetric 
shell,  and  thus  the  most  x-ray  output,  was  first  empirically  determined.  This  was  done 
by  selecting  a  wire  number  that  maximized  the  total  x-ray  power  output  from  a  series  of 
aluminum  wire  experiments  in  which  only  the  wire  number  and  wire  diameter  of  the  arrays 
were  varied  while  the  array  mass  was  kept  fixed  [2,3].  In  these  shots,  the  arrays  were  2  cm 
long,  4  cm  in  diameter,  each  contained  approximately  616  jug/cm,  and  the  wire  numbers 
were  32,  56,  70,  126,  180,  and  282.  The  highest  x-ray  power  was  obtained  from  the  180  wire 
load.  This  load  as  chosen  for  our  1-D  analysis  since  (1)  the  pinch  had  the  largest  cmpression 
and  hence,  it  was  the  closest  to  1-D,  i.e.,  it  had  the  minimum  of  ire  array  dynamics  and  (2) 
the  short  x-ray  pulse  minimized  errors  in  inferring  an  implosion  time.  However,  this  power 
was  highly  variable.  When  the  180  wire  shot  was  repeated  2  more  times,  large  variations  in 
the  total  x-ray  power  were  observed  (Fig.  1).  The  shots  had  nearly  identical  implosion  times 
(±2  ns),  but  they  x-ray  pulses  are  shifted  in  Fig.  1  for  clarity.  However,  in  shot  2693,  the 
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peak  x-ray  intensity  was  twice  that  measured  in  the  other  two  shots,  and  the  power  pulse 
had  a  substantial  late-time  tail.  Shot  2706  also  had  a  late-time  tail;  however,  shot  2636  did 
not.  Similar,  but  smaller,  variations  in  time  were  evident  in  the  recorded  aluminum  K-shell 
emissions,  which  are  shown  in  Fig.  2. 

Running  time  integrations  of  the  Fig.  1  powers,  which  are  shown  in  Fig.  3,  produced 
equally  large  variations  in  the  total  energy  radiated,  both  during  and  after  the  main  pulse. 
Twice  the  energy  was  radiated  in  shot  2693  than  in  the  other  two  shots,  but  half  of  this 
energy  was  late-time.  Shot  2706,  on  the  other  hand,  had  half  the  energy  in  its  main  pulse  as 
shots  2693  and  2636  did,  although  eventually  the  same  total  amount  of  energy  was  radiated 
as  in  shot  2636.  Approximately  the  same  amount  (~  500  kJ)  of  energy  was  radiated  in  the 
main  x-ray  pulse  of  shot  2636  as  was  seen  in  the  total  emissions  from  shot  2693,  but  shot 
2636  produced  essentially  no  late-time  emission. 

Initially,  when  an  analysis  of  Saturn  shots  2636,  2693,  and  2704  was  carried  out  using  a 
conventional  1-D  magnetohydrodynamics  (MHD)  calculation  based  on  classically  calculated 
transport  coefficients  [13]  and  the  Saturn  circuit  model,  two  to  four  times  more  energy  was 
found  to  be  radiated  in  the  experiments  than  had  been  supplied  in  the  calculation.  The  large 
energies  that  were  radiated  in  the  experiments  imply  that  even  larger  amounts  of  energy  had 
coupled  to  the  pinch.  As  our  analysis  will  show,  the  pinch  had  to  have  a  higher  impedance 
than  the  conventional  1-D  calculation  had  predicted. 

A  direct  measurement  of  impedance  in  experiments  conducted  on  the  Tomsk  IMRJ- 
IV  generator  [14]  had  confirmed  that  large  increases  in  resistivity,  and  therefore  in  energy 
coupling,  were  taking  place  late  in  the  pinch  implosion  as  it  began  to  assemble  on  axis. 
The  calculations  described  in  section  IV  of  this  paper  are  guided  by  this  observation  as 
well  as  by  the  x-ray  data  in  Figs.  1-3  from  shots  2636,  2693,  and  2704.  Large  increases  in 
the  resistance  of  the  plasma  as  seen  by  the  circuit  were  introduced  into  these  calculations 
slightly  in  advance  of  pinch  assembly,  and  the  turn-on  time  was  taken  to  be  relatively 
short  in  accord  with  the  Tomsk  experiments  and  with  the  Saturn  x-ray  data.  How  these 
resistance  increases  were  made  in  the  calculation  is  explained  in  the  following  section,  where 
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the  problem  of  coupling  a  set  of  1-D  fluid  equations  to  a  transmission  line  is  examined  and 
a  formula  relating  the  circuit  resistance  to  the  surface  values  of  the  plasma  resistivity  and 
current  density  is  derived. 

III.  ENERGY  COUPLING  THEORY 

The  coupling  of  pulsed-power  generators  to  Z  pinches  is  usefully  addressed  by  separately 
considering  the  dynamics  of  the  generator,  the  diode  region,  and  the  Z  pinch.  Maxwell’s 
equations  imply  that  transmission  line  energy  couples  to  a  pinch  at  its  outer  boundary.  The 
location,  b{z,t),  of  this  boundary  as  a  function  of  time  is,  therefore,  an  important  coupling 
parameter  even  though  it  may  not  be  an  easily  definable,  physical  variable.  Most  theoretical 
analyses,  however,  assume  it  to  be  well  defined.  In  a  non-fluid,  slug  model  description  of 
Z-pinch  implosions  [15],  for  example,  the  pinch  boundary  is  the  only  variable  employed.  It 
is  also  an  important  variable  in  the  definition  of  self-similar  solutions,  which  are  used  in  a 
variety  of  approximate  hydrodynamics  descriptions  of  Z-pinch  dynamics  [16].  It  was  also  the 
important  variable  in  a  recent  analysis  of  the  free  surface  dynamics  of  an  ideal  fluid  consisting 
of  imploding  wire-like  plasmas  [17].  It  was  used  to  describe  the  radiative  collapse  of  hydrogen 
and  helium  Z  pinches  [18],  and  it  is  the  key  variable  in  the  original  work  of  Bennett  [19],  Pease 
[20],  and  Braginskii  [21]  in  their  early  treatment  of  Z-pinch  equilibria.  It  is  less  well  defined 
in  more  recent  treatments  of  this  problem,  such  as  the  work  on  gas-imbedded  Z  pinches  [22] 
or  on  Bennett  equilibria  as  it  is  modified  once  microturbulence  limits  are  set  on  the  current 
flow  [12].  Most  importantly,  the  outer  pinch  boundary  is  sharply  defined  (minus  the  presence 
of  coronal  plasmas)  in  all  MHD  computer  calculations  of  Z-pinch  plasma  implosions,  and  it 
Tvill  be  assumed  to  be  equally  well-defined  in  the  analysis  that  follows. 

A.  Generator 

The  Saturn  generator  is  a  voltage-source  driven  transmission  line  in  which  electromag¬ 
netic  energy  flows  to,  couples  to,  and  reflects  from  the  Z-pinch  load.  Generally,  however,  the 
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generator  is  modeled  as  a  lumped  circuit  in  which  a  voltage,  V{t),  with  a  given  (measured) 
time  history  drives  a  time  varying  inductance,  L{t),  and  a  time  varying  resistance,  R{t): 

=  V(t)  (1) 

The  inductance  and  the  resistance  are  sums  of  time-independent  line  impedances  and  time- 
dependent  load  impedances:  L{t)  =  Lune  +  Lioadit)  and  R{t)  =  Rune  +  Rioadit)- 

The  energy  equation  for  this  simple  circuit  equation  is  found  by  multiplying  Eq.  (1) 
through  by  I: 

VI  =  +  t/f”")  +  {ikine  +  -Rtod  + 

m  J  (£/J“  +  £/|“)  +  Ru„P  +  (2) 

Here,  f/g”®  =  {l/2)LiineP  and  =  {\/2)LioadI^  are  the  magnetic  field  energies  stored  in 
the  transmission  line  and  in  the  diode  region  between  the  pinch  and  return  current  circuitry 
respectively.  This  equation  describes  how  once-stored  generator  energy  is  converted  into  the 
growth  of  line  and  diode  magnetic  field  energy.  Ohmic  heating  losses  in  the  line,  and  a  power 
flow  into  the  Z-pinch  load,  as  described  by 

P""*  =  +  R<oJ)l^-  (3) 


B.  Diode 

The  transmission  line  terminates  at  the  Z-pinch  diode.  In  principle,  therefore,  the  diode’s 
impedance  is  difiicult  to  compute  accurately,  and  it  must  be  determined  empirically.  How¬ 
ever,  if  it  is  assumed  that  the  diode  has  perfect  cylindrical  symmetry  and  perfectly  conduct¬ 
ing  walls,  that  the  region  between  the  pinch  and  the  return  current  is  a  vacuum,  and  that 
displacement  currents  within  the  diode  can  be  neglected,  then  an  approximate  evaluation 
of  the  electromagnetic  fields  in  the  diode  can  be  made  [23].  These  fields  are  needed  to 
determine  the  flow  of  electromagnetic  energy  into  the  diode  and  into  the  pinch,  and  thus. 
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they  relate  to  the  load  impedance.  Under  the  above  assumptions,  Maxwell’s  equations  can 
be  solved  subject  to  the  boundary  conditions  that  the  current  in  the  Z  pinch,  /(t),  and  the 
voltage  drop  across  the  entry  surface  to  the  diode,  Vp{t)  =  —  J^drEr{r,z  =  £),  are  known 
functions  of  the  time.  The  B-field,  in  this  case,  has  only  an  azimuthal  component,  Bg,  which 
is  determined  by  applying  Stokes’  theorem: 


B  =  Bgug  =  ^^^ug,  for  b  <r  <  R  and  0  <  z  <£,  (4) 

cr 

where  the  vacuum  diode  geometry  is  defined  by:  b,  the  outer  boundary  of  the  pinch,  R, 
the  radius  of  the  diode  region  and  of  the  return  current  path,  and  i,  both  the  length  of 
the  plasma  and  the  diode  region.  Maxwell’s  equations  then  require  that  the  E-field  has  a 
component  along  the  z-direction  of  the  pinch,  and  one  along  the  radial  direction,  Ur'. 
E  =  Erilr  +  E^u^,  where  (from  Ref.  [23]) 


E,  =  - 


<0, 


(5) 


and  /  =  dtl.  Since  the  walls  of  the  return  current  path  have  infinite  conductivity,  the  voltage 
drop  across  the  entrance  to  the  diode  region  also  appears  as  a  voltage  drop  across  the  load. 
The  full  voltage  drop,  j dzE^{r  =  b,z)  across  the  load,  however,  includes  a  back  EMF  and 
is  given  by 


2/ 

-^\n(Rlb)i  =  Vp-LioaJ, 


(6) 


where 


^ioad  (7) 

The  neglect  of  the  displacement  current  is  based  on  the  assumption  that  the  electric  field 
energy  density,  =  E^/{8Tr),  is  negligible  relative  to  u^°'^  =  Bg/{8Tr)  in  the  diode.  In 
this  case,  energy  conservation  is  expressed  locally  by  dtUg°'^  4-  V  •  S  =  0,  where  S  is  the 
Poynting  vector,  S  =  (c/(47r))E  x  B,  and  an  integration  of  this  equation  over  the  volume  of 
the  diode  region  then  yields  a  global  equation  for  energy  conservation  of  the  form. 
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pdiode  _  Q^jjdiode  _j_  pload  pwall 


where  is  the  energy  flux  into  the  diode  across  its  entry  surface  si  z  =  t 


/  surface 


dAS  -Uz 


=  -  2wrdr-^(Er{r,z  =  £)Be{r)^  =  Vpl, 


=  271  fdz  drrut^'^^  =  (10) 

JQ  Jb(t)  I 

Power  losses  into  the  return  current  walls,  P*""”,  are  found  by  integrating  the  Poynting 
flux  over  the  walls  of  the  return  current  can.  They  are  zero  in  this  analysis  (because  perfectly 
conducting  walls  were  assumed  in  it)  but  not  necessarily  in  experiments.  Power  flow  into  the 
Z  pinch,  has  two  parts,  one  coming  from  an  integration  of  the  Poynting  flux  across  the 
surface  of  the  pinch,  P^,  and  the  other,  from  the  magnetic  field  energy  density,  or  pressure, 
acting  on  the  moving  boundary  of  the  pinch: 


=  2-Kbk^^~-  -  -  +  P^  =  pL 

Stt  2  dt 


where  b  =  dtb.  cannot  be  determined  from  this  diode  analysis  alone.  The  fields  of 

Eqs.  (4)  and  (5)  satisfy  Eq.  (8)  tautologically.  To  determine  one  must  turn  to  an 

analysis  of  the  plasma  fluid  equations,  which  govern  the  main  characteristics  of  an  idealized 
Z-pinch  dynamics,  i.e.,  of  the  dynamics  that  begin  after  the  wires  have  broken  down  and 
plasma  has  fully  formed. 


C.  Z-pinch  Fluid  Theory 


Z  pinches  have  been  described  by  one-,  two-,  or  three-dimensional  fluid  dynamic  mod¬ 
els,  which  must,  in  each  case,  be  coupled  to  a  generator  circuit  equation  self-consistently. 
Plasmas  contain  four  kinds  of  fluid  energy:  kinetic,  uke  =  ion  thermal,  ulf^  — 
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{Z/2)nikBTi,  electron  thermal,  =  {'i/2)nekBTe,  and  an  internal  energy  of  ionization  and 
excitation,  «,on-  The  full  set  of  fluid  equations  describe  the  various  ways  that  these  energies 
are  either  partitioned  within  the  plasma  or  radiated  away.  Their  coupling  to  a  circuit  equa¬ 
tion  is  carried  out  so  as  to  maintain  power  flow  continuity  and  overall  energy  conservation. 
However,  while  fluid  dynamic  models  are  much  more  detailed  than  a  0-D  model;  they  are 
nevertheless  idealized  and  approximate  as  well,  since,  as  noted,  they  do  not  treat  the  early- 
time  dynamics  by  which  solid  wires  explode,  break  down,  and  become  plasmas  or  by  which 
gases  break  down.  In  addition,  computer  solutions  to  fluid  equations  are  resolution  limited. 

Equations  that  describe  the  build-up  and  decay  of  electron  and  ion  energy  densities  are 
derivable  directly  by  taking  energy  moments  of  the  electron  and  ion  kinetic  equations  [13]: 


dtUTot  +  V  •  (wfojVe  -f  Pe  •  Ve  -b  Qe) 


=  -ngCE  •  Ve  +  Qe  +  Rei  •  Ve  -b  Q*"®'  -  Pbrems,  (12) 

+  V  •  (wkfVt  +  Pi  •  Vi  -b  qi) 

=  UiZeE  •  Vi  -b  Qi  +  Rie  •  Vi,  (13) 

where  =  ^PeV^  -b  is  the  sum  of  electron  kinetic  and  thermal  energies,  Vg  is  the 
electron  fluid  velocity,  Pg  =  Peh  +  Tig  is  the  electron  pressure  tensor,  qg  is  the  electron  heat 
flux  vector,  E  is  the  electric  field,  Qg  is  the  elastic-collision  energy  transfer  rate  between 
electrons  and  ions,  Rgi  is  the  elastic-collision  momentum  transfer  rate  between  electrons 
and  ions,  Q*”®'  is  the  inelastic-collision  energy  transfer  rate  between  electrons  and  ions,  and 
Pbrems  is  the  bremsstrahlung  energy  loss  rate.  The  quantities,  v,-,  Pj  =  pjl2  +  H,-,  q,-, 
Qi,  and  R,e  are  similarly  defined  for  the  ions.  The  term,  Qy'®^  describes  the  rate  of  increase 
in  the  ion  excitation  and  ionization  energy  density,  Ujon,  which  is  given  by 


diUion  +  V  •  (UionVj)  =  Qf  -  Pi  -  Prr- 


Here,  pi  and  Prr  are  line  and  recombination  radiation  energy  density  loss  rates,  respectively. 
The  fluid  equations  of  a  plasma  are  defined  in  terms  of  locally  defined  center  of  mass 
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and  relative  (or  drift)  velocities,  v/  and  respectively,  in  terms  of  which  the  electron  and 
ion  fluid  velocities  are  given  by 

Ve  =  V/  +  =  V/  +  Vd,  Vi  =  V/  -  ^Vd  =  Vf.  (15) 

where  p  =  pe  pi  =  pi.  An  equation  describing  total  energy  conservation  can  then  be 
derived  in  terms  of  Vf  and  Vd  by  adding  Eqs.  (12)  and  (13)  together  and  using  the  condition, 
Qe  +  Rei  •  Vg  =  —Qi  —  Rie  *  Vj,  for  energy  conservation  in  elastic  collisions: 

dtUTot  +  V  •  (uTot^Tf  +  P  •  V/  +  q  +  U^ot'^d  +  Pe  '  Vd) 

=  j  •  E  -  Prad,  (16) 

where,  urot  =  't^Tot  +  “rot  +  “ioni  9  =  9e  +  ^i,  j  =  —rieevd  is  the  electric  current  density,  and 
Prad  =  Pi  +  Ptt  +  Pbrems  Is  the  total  radiation  loss  rate. 

The  plasma  couples  to  the  electromagnetic  field  by  way  of  the  j-E  term  in  Eq.  (16),  which 
also  appears  in  the  energy  equation  derived  from  Maxwell’s  equations.  Since  displacement 
currents  are  small  in  a  plasma,  the  left  hand  side  of  this  equation  has  the  same  form  as  the 
energy  equation  in  the  diode: 

5tug  +  V.SP'  =  -j.E,  (17) 

where  Ug  =  B^/ (Stt)  is  the  magnetic  field  energy  density  stored  within  the  plasma.  Finally, 
the  addition  of  Eq.  (16)  to  Eq.  (17)  produces  an  equation  expressing  local  total  energy 
conservation: 

dt  (uTot  +  ufj  +  V  •  (uTot^f  +  P  •  V/  +  q  +  Mf^Vd  +  Pg  •  Vd  +  S^')  =  -Prad-  (18) 

The  coupling  of  the  plasma  fluid  to  the  generator  is  now  accomplished  by  integrating  Eq.  (18) 
over  the  volume  of  the  plasma.  The  resulting  equation  connects  the  build-up  of  energy  within 
the  plasma  to  the  flow  of  energy  into  the  plasma  from  the  diode, 

P^=J^{U?,t+Ut)+Pr^,  (19) 
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and  thus  to  the  flow  of  energy  into  the  diode  from  the  generator  when  combined  with  Eq.  (8), 
where 


ri  fb{t) 

U^ot  =  27:/  dz  I  rdruTot  (20) 

Jo  Jo 

and  Ug  and  Prad  are  similarly  defined  integrations  over  Ug  and  Prod  respectively. 

When  the  surface  of  the  integration  volume  is  taken  to  lie  just  inside  of  the  plasma 
boundary,  the  fluid  variables  are  nonzero  and,  in  particular,  the  electric  field  satisfies  an 
Ohm’s  law,  which,  in  a  cylindrically  symmetric  plasma,  is  given  by 

E.  = ^  .  (21) 

C  rigC  I  Tlee  )  r=b-e 

Then,  because 

Z  ^  •  (“TotV/)  }  =  (22) 

and  because  the  current  flows  parallel  to  the  surface  of  the  integration  volume,  whose  normal 
is  in  the  radial  direction,  it  follows  from  Eqs.  (18)  and  (19)  that 

ffiU?o.  +  US)+Pr.^=  (23) 

Taken  together,  therefore,  Eqs.  (3),  (21),  and  (23)  imply  that 

P"^  =  (j^  +  =  (24) 

-  lim  ( 2jr66£^  +  2nb£(qr  +  (P  -  v.).  -  -^E^Bg)  )  I 

«-»o+ 1  St:  '■  47:  'J  lr=6-£ 

Eq.  (25)  can  be  used  to  determine  Rioad  and  Liaad  provided  that  the  plasma  boundary  is 
located  where  pressure  forces  make  a  negligible  contribution  to  the  equation.  On  physical 
grounds,  gr  is  strictly  zero  at  the  plasma  outer  surface;  however,  one  cannot,  in  principle,  set 
the  surface  pressure  to  zero  unless  one  sets  rie  =  0  and  nj  =  0.  However,  the  current  density 
would  then  go  to  zero  unless  one  allows  Vd  at  the  surface  to  become  infinitely  large.  In 


Ref.  [12],  it  is  this  behavior  that  argues  for  the  growth  of  micro-instabilities  and  anomalous 
resistivity  near  to  and  at  the  plasma  surface.  This  behavior  also  illustrates  one  of  the 
ambiguities  associated  with  producing  a  physically  meaningful  definition  of  an  outer  Z-pinch 
boundary.  However,  on  setting  (P  •  v/)r  =  0,  one  finds  first  that 


1  dLioad  r2  I 

2  dt 


Rioadl^  =  -  lim  { 27rbbip-  -  ^E,Be)  ] 

«-+0+l.  StT  47r  r=b-e 


and  then,  on  substituting  Eqs.  (4)  and  (21)  into  this  equation,  one  infers  that 

^Bload  ^b 

dt  c^b  dt 


Rloadl^  —  "^T^bi 


{  —  —i 

1 2'Kb  Ttge  V 


a±—^  +  ^AdriksTg) 
n^e 


Eq.  (26)  for  Lioad  is  consistent  with  Eq.  (7),  while  Eq.  (27)  can  be  used  to  express  Rioad 
in  a  more  physically  meaningful  form  through  the  definition,  first  of  all,  of  two  areas,  Aj 
and  Ath,  and  a  thermal  velocity,  Vth'- 

J  j 

^  jz{r  =  by  ^  jth{r=^by  ^  (28) 

and  then,  of  a  thermal  current  density, 

jth  =  rieevth  =  (29) 

dimensionless  transport  coefficients,  dj.  and  by  Ref.  [13] 

mrie ..  ^  A  nJzVTp 

= - ax,  Pa  =  n^pA,  «x  = - ^rgKj.,  (30) 

Tg  771 

and  a  dimensionless  temperature  gradient, 

^tr  —  rp  d/Re  •  (31) 

J-e  r=b 


One  then  finds  that 
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where  the  resistivity,  rjo,  is  defined  by 

r}o  =  m/{nee^Te).  (33) 

The  above  expression  for  Rioad  has  the  added  benefit  that  it  guarantees  the  continuity 
of  across  the  pinch  boundary.  Thus,  one  finds,  on  combining  Eqs.  (4),  (5),  (7),  and  (21), 
that 

E,{r  =  b  +  €)  =  ]{V,-Lioadi) 


- 

c\cb 


)  Uet 


a±—  +M 

L  flgC 


=  E^ir  =  6  -  c).  (34) 


r=:b—( 


Using  Eq.  (26),  this  equation  can  be  rewritten  as 

v;  =  +  —Ui—  -  MiksT.) 

dt^  >  n.e\  n.e 


r=b 


However,  by  definition  of  Vp,  the  circuit  equation  [Eq.  (1)]  can  also  be  rewritten  as 

dl 


which  identifies  Vp  as 


V{t)=^Lii„e-^+ReineI  +  Vp, 


Vp  —  {Etoadl^  "h  RtoadE 


(35) 


(36) 


(37) 


The  comparison  of  Eq.  (35)  with  Eq.  (37)  provides  an  expression  for  Rtoad  that  is  identical 
to  Eq.  (33). 

Finally,  a  third  expression  for  Rtoad  can  be  derived,  which  provides  an  important  check  on 
the  accuracy  of  Z-pinch  fluid  dynamics  calculations  when  it  is  compared  to  the  calculation 
of  Rtoad  from  Eq.  (33).  To  obtain  this  expression,  one  integrates  Eq.  (17)  over  the  volume 
of  the  Z-pinch  plasma  to  find  [using  Eq.  (11)]  that 

f6(t) 


dU^ 


-  =  -2ne  /  rdrj.E,. 

Jo 


dt  JO 

One  can  then  substitute  Eq.  (3)  into  Eq.  (36)  and  solve  for  Rtoad’ 

dLtoad 


2  dt 


(38) 


(39) 


In  this  equation,  the  volume  integral  over  j^E^  =  j  •  E  includes  both  vy  •  j  x  B  work  and 
Ohmic  heating  by  virtue  of  Ohm’s  law  [Eq.  (21)]. 
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IV.  MHD  MODELING 


The  number  of  approximations  built  into  the  above  analysis  indicates  the  phenomenolog¬ 
ical  nature  of  Z-pinch  modeling  (even  when  computers  are  used  to  solve  sets  of  nonlinearly 
coupled,  energy  conserving  circuit  and  fluid  equations).  It  is  not  entirely  surprising,  there¬ 
fore,  when  such  calculations  contain  insufficient  resistive  coupling  to  reproduce  the  energy 
coupling  inferable  from  the  Saturn  experiments.  Based  on  the  above  analysis,  this  cou¬ 
pling  deficiency  can  be  overcome  in  1-D  modeling  through  the  introduction  of  an  additional 
(anomalous)  surface  resistivity; 


\r=:b 


^  ^XVo  .  “t"  VAnomalous  , ' 
r=o  r=o 


(40) 


The  magnitude  and  temporal  characteristics  of  rjAnomaious  within  the  plasma,  however,  are 
not  determined  by  coupling  considerations;  they  must  be  inferred  indirectly,  if  possible, 
from  the  influence  that  ^Anomalous  has  on  the  stagnation  dynamics  of  a  Z  pinch  and  from  the 
ability  of  this  dynamics  to  replicate  the  x-ray  output  trends  seen  in  experimental  data.  In 
this  section,  some  results  are  presented  from  a  phenomenological  study  of  the  Saturn  data 
in  Figs.  1-3  within  the  context  of  1-D  MHD  calculations. 

This  study  addressed  such  issues  as:  how  much  energy  actually  couples  to  a  pinch  when 
additions  to  the  circuit  resistance  are  made,  how  large  an  increase  in  circuit  resistance  is 
needed  in  order  to  double  the  amount  of  coupled  energy,  when  and  how  rapidly  can  this 
resistance  be  applied  during  the  course  of  the  implosion,  how  sensitive  is  the  calculation  to 
variations  in  the  enhanced  resistivity,  and,  finally,  what  rate  of  energy  transport  is  required 
within  a  plasma  for  enhanced  energy  inputs  to  show  up  as  enhanced  x-ray  outputs?  These 
questions  were  investigated  within  the  context  of  a  1-D  radiative  hydrodynamics  calculation 
by  (1)  specifying  a  functional  form  for  the  anomalous  circuit  resistance,  (2)  varying  the 
parameters  of  the  functional  form,  and  (3)  comparing  the  calculated  x-ray  emission  to  the 
experimental  x-ray  pulses  shown  in  Fig.  1. 

One  learns  by  this  procedure  that  the  circuit  resistance  cannot  be  sizably  increased 
until  the  pinch  nears  assembly  since  otherwise  observable  amounts  of  soft  x  rays  are  emitted 
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earlier  than  when  they  are  seen  experimentally.  This  finding  is  consistent  with  the  impedance 
measurement  reported  in  Ref.  [14].  One  also  finds  limits  on  the  rate  at  which  energy  can 
be  coupled  effectively  to  the  plasma  because  of  the  large  impact  that  this  rate  has  on  the 
hydrodynamics.  These  limits  are  set  by  both  the  rate  at  which  energy  is  transported  into 
the  plasma  and  on  the  rate  at  which  it  is  radiated  from  the  plasma.  Only  when  these  two 
rates  are  sufficiently  large  does  the  calculated  x-ray  output  move  in  the  direction  of  the  data 
shown  in  Figs.  1-3. 

The  calculations  described  in  this  section  were  computer  intensive;  the  aluminum  model 
contained  136  energy  levels,  22  K-shell  states,  77  L-shell  states,  125  K-shell  lines,  and  202 
L-shell  lines,  and  the  calculations  employed  a  tabular  lookup,  escape  probability  method  to 
carry  out  both  the  kilovolt  and  subkilovolt  photon  transport  [24].  Nevertheless,  77  states 
do  not  provide  an  accurate  description  of  the  multiplet  structure  of  the  L-shell  ionization 
stages;  thus,  L-shell  emissions  are  not  as  accurately  transported  within  the  plasma  as  they 
should  be.  The  transport  of  soft  x  rays  is  more  diffusive  than  that  of  kilovolt  x  rays,  and  thus 
escape  probability  methods  are  more  appropriate  for  the  kilovolt  emissions  than  subkilovolt. 
For  this  reason,  we  introduced  two  parameters  into  our  calculations  in  order  to  be  able  to 
phenomenologically  increase  and  vary  the  rates  of  soft  x-ray  energy  transport  and  partially 
compensate  for  the  missing  atomic  structure.  A  multiplier,  o;,  on  the  electron  heat  flux 
was  introduced  as  a  way  of  approximately  simulating  the  influence  of  soft  x-ray  diffusive 
transport  on  the  pinch  dynamics,  and  a  multiplier,  P,  was  included  in  the  calculation  of 
L-shell  emissivities  to  compensate  for  missing  distributions  of  L-shell  oscillator  strength. 

In  all,  four  modifications  to  the  MHD  were  made  in  this  investigation  of  enhanced  energy 
coupling.  The  load  impedance,  the  rate  of  diffusive  energy  transport,  and  the  rate  of  soft 
x-ray  emission  were  all  increased  as  described  above.  In  addition,  the  value  of  the  anomalous 
resistivity  was  specified  throughout  the  plasma.  For  lack  of  specific  theoretical  guidance, 
we  assumed  a  fixed,  gradual,  linear  reduction  by  roughly  two  orders  of  magnitude  from  its 
surface  value  in  the  value  of  rjAnomaious  throughout  the  plasma  interior.  By  this  assump¬ 
tion,  resistivity  levels  remained  anomalously  high  throughout  the  plasma,  which  produced 
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increases  in  the  rate  of  magnetic  field  and  current  diffusion  and,  therefore,  in  the  amount 
of  Ohmic  heating  interior  to  the  plasma.  The  impact  that  this  resistivity  assumption  had 
on  x-ray  conversion  efficiencies  was,  therefore,  positive,  since  it  aided  the  energy  transport 
into  the  plasma,  i.e.,  steepening  of  the  resistivity  gradient  at  the  surface  of  the  pinch  serves 
to  reduce  the  degree  of  current  diffusion  and  Ohmic  heating  within  the  plasma.  Energy 
must  then  be  transported  by  secondary  means  into  the  plasma  interior  (by  pressure  forces, 
electron  heat  conduction,  or  soft  x-ray  transport)  in  order  to  be  radiated  efficiently.  The 
impact  that  varying  the  resistivity  profile  has  on  the  calculations  was  not  addressed  in  this 
work.  The  resistivity  was  reduced  from  zone  to  zone  as  described  above,  and  this  contouring 
was  then  kept  the  same  in  all  of  our  calculations. 

By  applying  the  above  phenomenology,  we  were  able  to  replicate  several  important  fea¬ 
tures  of  the  Figs.  1-3  Saturn  data.  For  example,  when  the  resistance  time  history  shown  in 
Fig.  4  was  employed  in  a  100  zone  1-D  calculation,  it  produced  the  total  x-ray  output  pulse 
shown  in  Fig.  5.  This  pulse  approximates  both  the  early-  and  late-time  features  of  the  x-ray 
pulses  shown  in  Fig.  1.  The  heat  flux  multiplier,  a,  had  a  value  of  600  and  a  value  of 
3  in  this  calculation.  The  sensitivity  of  the  calculation  to  and  the  validity  of  the  a  value 
used  are  discussed  later  in  this  section.  The  resistance  in  Fig.  4  is  composed  of  two  parts. 
One  (labeled  Braginskii  in  the  figure)  was  calculated  from  the  surface  value  of  the  classical 
resistivity,  a±r)o,  using  Eq.  (32),  while  the  other  was  a  prescribed  input  to  the  calculation. 
The  prescribed  resistance  was  taken  to  rise  exponentially  with  a  risetime  of  5  ns,  to  peak  at 
172  ns  at  a  value  of  0.658  and  to  fall  exponentially  with  a  falltime  of  30  ns.  The  surface 
value  of  TjAnomaious  that  was  used  in  the  fluid  equations  was  calculated  from  the  prescribed 
resistance  using  Eqs.  (32)  and  (40). 

Because  the  circuit  resistance  had  an  extended  decay,  a  substantial  amount  of  energy 
coupled  to  the  plasma  following  the  main  x-ray  pulse,  and  this  late-time  coupling  produced 
the  foot  (or  tail)  to  the  main  pulse  seen  in  Fig.  5.  The  emission  spike  at  the  end  of  the  tail 
was  produced  by  a  secondary,  late-time  implosion,  an  event  often  seen  in  1-D  calculations.  In 
this  instance,  the  presence  of  this  spike  correlates  with  the  more  subdued  late-time  spiking 
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seen  in  the  x-ray  pulses  of  the  (3-D)  shots  2693  and  2706  (and  noted  in  Fig.  1).  In  this 
100  zone,  fiducial  calculation,  671  kJ  coupled  to  the  pinch  in  245  ns;  369  kJ  of  this  coupled 
energy  was  emitted  in  the  main  x-ray  pulse,  and  292  kJ,  in  the  tail.  Overall,  therefore, 
661  kJ  of  the  671  kJ  input  was  radiated.  It  was  divided  into  143  kJ  of  kilovolt  and  518 
kJ  of  subkilovolt  emissions.  Thus,  in  the  245  ns  time  interval  over  which  the  calculation 
extended,  substantially  all  of  the  energy  that  coupled  to  the  pinch  was  radiated,  and  it 
had  to  be  in  order  to  be  comparable  to  the  total  x-ray  output  measured  in  shots  2636  and 
2706  (Fig.  3).  For  the  calculation  to  produce  this  result,  however,  a  nearly  three  order  of 
magnitude  increase  in  the  circuit  resistance  was  needed,  in  agreement  with  the  resistance 
increases  predicted  by  the  0-D  calculations  described  in  Ref.  [25].  Finally,  we  note  that, 
while  a  comparable  amount  of  energy  coupled  to  the  pinch  in  this  calculation  as  in  the 
experiments,  there  was  much  less  K-shell  emission  (by  a  factor  ~  2.5)  in  the  experiments 
than  in  the  calculation.  Moreover,  the  calculated  K-shell  pulse  had  less  late-time  structure 
than  was  seen  experimentally  in  shots  2693  and  2706  (Fig.  2). 

Fig.  4  contains  two  plots  of  Rtoadit).  The  one  labeled  surface  calculation  was  obtained 
from  Eqs.  (32)  and  (40)  with  r^AnomcUous  derived  from  the  prescribed  part  of  Rioad-  The  dotted 
curve,  labeled  volume  calculation,  was  calculated  from  Eq.  (39).  When  these  two  calculations 
are  in  close  agreement  at  a  given  time,  they  indicate  that  an  accurate  spatial  integration 
of  the  magnetic  field  diffusion  equation  has  been  carried  out  at  that  time.  This  accuracy 
depends  on  the  zoning  in  the  calculation,  and  deviations  between  the  curves  temporarily 
occurred  in  the  calculation  late-time  because  irregular  spacings  between  zones  had  been 
generated  by  the  Lagrangian  fluid  dynamics.  Thus,  the  quality  of  the  agreement  between 
the  surface  and  volume  calculated  values  of  Rtoad  (which  was  generally  excellent  for  all  of 
the  calculations  described  in  this  section)  can  be  degraded  when  the  number  of  zones  used 
in  the  calculation  is  too  small  or  when  their  spacing  becomes  too  irregular. 

The  above,  fiducial,  calculation  was  constructed  so  that  the  MHD  dynamics  would  couple, 
transport,  and  radiate  energy  at  comparable  rates,  and  retard  expansion  of  the  plasma 
following  stagnation.  Under  these  conditions,  the  plasma  continues  to  radiate  while  the 
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coupling  remains  strong,  producing  a  tail  on  the  x-ray  pulse  as  seen  in  Fig.  5.  On  the 
other  hand,  when  the  enhanced  resistance  is  taken  to  fall  as  fast  as  it  rises  (by  replacing  the 
resistance  of  Fig.  4  by  a  comparably  fast  rising  and  falling  Gaussian  pulse,  for  example),  only 
a  single  x-ray  pulse  is  generated,  similar  to  the  one  observed  in  shot  2636  (Fig.  1).  X-ray 
emission  cuts-off  once  the  energy  coupling  enhancements  are  terminated.  In  this  case,  the 
inductive  coupling  between  pinch  and  diode  is  reversed  as  the  pinch  maintains  an  outward 
expansion. 

The  x-ray  emission  generated  in  this  fiducial  calculation  depended  critically  on  four  pa¬ 
rameters:  the  magnitude  of  the  heat  flux  multiplier,  the  magnitude  of  the  L-shell  emissivity 
multiplier,  the  maximum  value  of  Rioad,  and  the  time  relative  to  plasma  assembly  at  which 
Rioad  achieves  its  maximum  value.  Consequently,  we  determined  the  sensitivity  of  the  cal¬ 
culation  to  each  of  these  parameters  by  varying  each  parameter  in  turn  about  the  value  it 
had  in  the  fiducial  calculation.  The  effects  of  these  single  parameter  variations  on  energy 
coupling  and  x-ray  output  are  presented  in  Figs.  6-15. 

In  order  to  determine  the  dependence  of  x-ray  production  on  the  strength  of  diffusive  en¬ 
ergy  transport,  a  physically  based  range  over  which  a  can  be  varied  needed  to  be  determined. 
This  was  done  by  comparing  the  heat  flux,  Qrad  —  —i^R^rTe,  due  to  blackbody  radiation  to 
the  electron  heat  flux,  =  —KgdrTe,  where  /Cg  =  k^[nJiyTeTe) / rrie  and  kr  =  \QaTl/{ZpKjt. 
Here,  kg  is  a  dimensionless  transport  coefficient  [13],  p  is  the  plasma  density,  a  is  the  Stefan- 
Boltzmann  constant,  and  Kr  is  the  Rosseland  mean  opacity  per  unit  mass  density.  The 
Rosseland  mean  for  aluminum  was  recently  calculated  [26].  It  has  values  ranging  between 
10^  and  10'*  cm^/g  at  the  plasma  densities  and  temperatures  that  were  present  in  our  calcula¬ 
tions  prior  to  Z-pinch  assembly.  When  employed  in  our  calculations,  these  temperature  and 
density  values  predicted  blackbody  heat  fluxes  that  were  more  than  four  orders  of  magnitude 
larger  than  the  electron  heat  flux.  Consequently,  we  varied  a  up  to  values  of  2000. 

In  all  the  calculations  in  which  a  was  varied,  the  energy  that  coupled  to  the  pinch  in  200 
ns  was  only  slightly  affected,  ranging  narrowly  between  500  and  550  kJ.  The  ability  of  the 
plasma  to  radiate  this  energy,  however,  depended  sensitively  on  its  ability  to  conduct  the 
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coupled  energy  into  its  interior.  The  three  energy  curves  shown  in  Figs.  6  and  7  demonstrate 
this  point.  The  curve  labled  180  ns  in  Fig.  6  shows  the  total  x-ray  energy  radiated  in  the 
main  pulse,  which  terminated  at  180  ns  in  these  calculations.  The  other  curves  in  Figs.  6  and 
7  show  how  the  total  and  K-shell  x-ray  outputs  varied,  respectively,  over  a  200  ns  period. 
At  first,  the  total  emission  rises  as  more  energy  is  conducted  into  the  plasma,  but  then  it 
falls  because  sufficient  energy  has  been  conducted  inwards  to  ionize  the  plasma  into  He-like 
ground  states.  At  this  point,  L-shell  emissions  turn  down  while  K-shell  emissions  continue 
to  rise.  Once  these  K-shell  emissions  slow  their  rise,  saturate,  and  then  fall  for  a  >  600, 
L-shell  emissions  can  return  to  their  peak  values. 

The  influence  of  the  L-shell  emission  rate  on  x-ray  outputs  is  shown  in  Figs.  8  and  9. 
The  energies  plotted  are  emissions  that  were  calculated  over  the  first  200  ns  of  the  pinch 
dynamics.  Up  to  a  level  of  2,  both  the  K-  and  L-shell  emissions  increase  with  jS  because  an 
increased  cooling  rate  moderates  the  fluid  dynamics  during  plasma  assembly  and  inhibites 
the  subsequent  regeneration  of  kinetic  energy.  However,  as  ^  becomes  larger  than  2,  K-shell 
emissions  decrease  rapidly  because  a  growing  amount  of  L-shell  cooled  plasma  is  unable  to 
ionize  into  the  K-shell.  Throughout  these  variations  in  L-shell  cooling  rate,  energy  coupling 
to  the  plasma  remained  relatively  unchanged. 

Both  the  shape  and  the  magnitude  of  the  measured  x-ray  pulses  place  severe  restrictions 
on  the  time  history  of  the  load  resistance  that  is  needed  to  replicate  these  pulses.  When  too 
much  heating  occurs  too  early  in  the  implosion,  the  plasma  emits  x-rays  prematurely  and 
produces  a  pulse  with  an  early-time  foot,  which  was  not  seen  in  the  experiments.  Moreover, 
too  much  early-time  heating  causes  the  plasma  to  overheat  during  run-in,  reducing  its  com¬ 
pressibility  during  assembly  and,  in  turn,  its  x-ray  conversion  efficiency.  Thus,  if  large  and 
rapid  increases  in  pinch  resistance  occurred  in  the  experiments,  these  increases  must  have 
taken  place  near  to  the  time  of  pinch  assembly. 

Given  these  considerations,  two  sets  of  calculations  were  carried  out  to  determine  the 
sensitivity  of  energy  coupling  and  the  accompaning  x-ray  production  to  the  size  and  timing 
of  circuit  resistance  increases.  The  results  of  these  calculations  are  shown  in  Figs.  10-14.  In 
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these  calculations,  the  shape  of  the  resistance  pulse  was  kept  the  same  as  shown  in  Fig.  4, 
i.e.,  the  rise  and  fall  times  were  kept  unchanged,  while  variations  were  made  to  its  peak 
value  and  to  the  time  at  which  it  reached  this  value,  while  a  and  ^  were  held  at  values  of 
600  and  3  respectively. 

In  one  set  of  calculations,  the  maximum  value  of  the  load  resistance,  was  varied 

while  the  rise  and  fall  times  and  the  time  to  peak  resistance  remained  unchanged.  The  energy 
coupling  and  radiative  outputs  calculated  during  the  first  200  ns  of  the  pinch  dynamics  are 
shown  in  Figs.  10-12.  The  total  energy  coupled  to  the  plasma,  which  is  shown  in  Fig.  10, 
has  been  broken  into  an  inductive  and  a  resistive  part.  The  inductive  part  is  defined  by 

(41) 

and  the  resistive  part  by 

B,„(t)  =  f‘  dt' (42) 
Jo 

Initially  both  of  these  couplings,  and  thus  the  total  coupling,  increase  with  However, 

once  R^ad  becomes  larger  than  0.3fl,  the  inductive  coupling  falls  at  roughly  at  the  same 
rate  as  the  resistive  coupling  rises,  causing  the  total  coupled  energy  to  remain  relatively  flat. 
Beyond  a  peak  of  0.60,  the  rise  and  fall  rates  of  Eres  and  Ein^  respectively,  increase,  and  this 
change  correlates  with  the  decline  in  calculated  x-ray  emissions,  both  total  and  K-shell  that 
are  shown  in  Figs.  11  and  12.  This  behavior  indicates  that,  when  the  amount  of  resistive 
heating  is  too  large,  the  plasma  cannot  radiate  coupled  energy  fast  enough  to  prevent  its 
conversion  into  kinetic  energy. 

Finally,  the  sensitivity  of  coupled  and  radiated  energy  to  the  time  at  which  the  load 
resistance  reaches  its  maximum  value  is  illustrated  in  Figs.  13  and  14.  In  these  calculations, 
R^d  was  held  at  a  value  of  0.6530,  and  the  time  to  peak  resistance,  to  +  At,  was  varied 
axound  its  Fig.  4  value  in  which  to  =  172ns.  The  total  and  resistive  energies  that  were 
coupled  in  these  calculations  over  200  ns  are  plotted  in  Fig.  13  as  a  function  of  At.  Fig.  14 
contains  similar  plots  of  the  calculated  total  and  K-shell  radiated  energies.  As  noted  above. 
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these  figures  demonstrate  that  too  much  heating  too  soon  in  the  implosion  reduces  the 
plasma  compressibility  and  the  radiative  cooling  at  assembly.  Extra  energy  inputs  convert 
instead  into  kinetic  energy.  Hence,  the  radiative  output  rises  rapidly  as  >  0,  and  K-shell 
emission  is  maximized  at  At  =  0.  When  At —  2  ns,  the  circuit  resistance  peaks  at  the  time 
of  peak  emission,  causing  the  K-shell  emission  to  fall  to  half  of  its  maximum  value.  We  note 
finally  that  the  extreme  sensitivity,  which  is  seen  in  all  of  the  above  calculations,  to  small 
changes  in  the  size,  shape,  and  timing  of  the  enhanced  circuit  resistance  is  consistent  with 
the  large  shot-to-shot  variations  seen  experimentally  (Fig.  1)  in  the  x-ray  data. 


V.  PLASMA  COLLISIONALITY 

In  Ref.  [12],  a  resistivity  formula  based  on  the  growth  and  saturation  of  lower-hybrid 
drift  wave  micro-instabilities  at  the  surface  of  a  pinch  was  proposed  and  its  effects  on  the 
equilibrium  states  and  the  radiative  collapse  of  Z-pinches  was  investigated.  This  formula 
was  written  in  the  form  of  a  classical  resistivity  with  an  enhanced  collision  frequency,  u*, 
produced  by  an  electron  drift-wave  interaction: 

where  Uce  and  Ud  are  the  electron  and  ion  plasma  frequencies  respectively,  =  je/inee) 


is  the  electron  drift  velocity,  and  is  the  average  ion  thermal  speed.  This 

resistivity  is  temperature  insensitive,  and  it  increases  as  the  surface  density  decreases  in 
contrast  to  Spitzer  resistivity,  which  is  essentially  density  independent  and  which  decreases 
with  temperature  increases.  As  electromagnetic  energy  flows  into  and  couples  to  a  pinch,  its 
surface  temperature  rises  and  its  surface  density  falls.  As  these  changes  take  place,  Spitzer 
resistivity  decreases,  while  drift-wave  resistivity  increases. 

Because  a  drift-wave  resistivity  has  properties  in  common  with  the  surface  resistivity 
employed  in  the  above  calculations,  it  was  of  interest  to  compare  the  two  to  determine 
how  closely  they  would  match  up  in,  for  example,  the  above,  fiducial  calculation.  This 
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comparison  is  made  in  Fig.  15,  which  shows  the  time  evolution  of  four  resistivities.  The 
resistivity  at  the  plasma  outer  boundary,  labeled  r]r=b  in  Fig.  15,  was  calculated  from  Reoad- 
r]r=b  =  Rtoadl/i^jzif  =  ^))-  The  curve  labeled  7)r=ro  shows  the  resistivity  at  the  nominal 
inner  boundary  of  the  shell-like  imploding  plasma.  It  is  shown  to  illustrate  the  resistivity 
contouring  that  was  employed  in  the  calculation.  The  classical  resistivity,  labeled  r^Bragi  was 
calculated  at  the  plasma  outer  boundary  using  the  formula,  qBrag  =  lr=6  [13].  Finally, 
the  resistivity,  labeled  t]*,  was  calculated  at  the  outer  boundary  of  the  plasma  using  the 
drift- wave  formula  (Eq.  (43)).  As  expected,  it  rose  more  than  two  orders  of  magnitude; 
however,  its  rise  occurred  late  in  time,  well  after  the  rise  in  r)r=b  that  was  needed  to  produce 
the  x-ray  output  shown  in  Fig.  5. 


VI.  SUMMARY 

Because  both  Maxwell’s  equations  and  the  fluid  equations  describe  local  flows  of  energy, 
Z-pinch  to  pulsed-power  generator  coupling  is  determined  by  the  continuity  of  the  transmis¬ 
sion  line  energy  flow  across  the  outer  boundary  of  the  pinch.  Increases  in  Z-pinch  resistance, 
therefore,  are  related  to  increases  in  the  plasma’s  surface  resistivity  (and  to  its  surface  area). 
Resistivity  increases  can  occur  when  surface  conditions  are  set  up  during  a  pulsed-power  dis¬ 
charge  that  stimulate  the  growth  of  plasma  micro-instabilities  (which  increase  the  plasma’s 
surface  collisionality)  as  well  as  the  growth  of  Rayleigh-Taylor  or  other  implosion  instabilities 
(which  are  calculated  to  increase  a  pinch’s  surface  area).  In  this  paper,  the  impact  on  the 
fluid  dynamics  of  augmenting  the  surface  collisionality  was  investigated  in  order  to  increase 
the  energy  coupling  and  to  evaluate  how  well  such  augmentations  are  able  to  replicate  the 
x-ray  behavior  seen  in  a  set  of  identical  Saturn  experiments.  A  phenomenological  approach 
to  modeling  this  coupling  was  adopted  for  lack  of  a  first-principles  theory  that  accurately 
describes  how  micro-instability  growth,  for  example,  affects  a  Z-pinch’s  surface  resistivity. 

Four  things  were  learned  from  this  1-D  analysis  of  the  Saturn  experiments.  One,  three 
orders  of  magnitude  increases  in  resistance  were  needed  to  produce  comparable  amounts  of 
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energy  coupling  as  were  seen  experimentally.  Two,  the  rate  at  which  energy  is  transported 
within  the  plasma  and  the  rate  at  which  it  is  radiated  had  to  be  comparable  to  the  energy 
coupling  rate  for  the  added  energy  inputs  to  end  up  as  added  x-ray  outputs.  Three,  for  the 
calculations  to  approximately  replicate  the  time  history  of  a  given  measured  x-ray  pulse,  only 
a  highly  restricted  range  of  resistance  enhancements  could  be  employed  in  these  calculations. 
In  particular,  the  presence  or  absence  of  a  late-time  tail  depended  on  the  decay  rate  of  the 
enhanced  resistance.  And  four,  the  extreme  sensitivity  of  the  calculations  to  variations  in 
(1)  the  rates  of  energy  transport  and  loss,  (2)  the  peak  value  of  the  enhanced  resistance, 
and  (3)  the  onset  time  of  the  rapid  resistance  enhancements  were  consistent  with  the  large 
shot-to-shot  variations  seen  in  the  Saturn  experiments. 

Our  calculations  showed  that  significant  changes  in  the  fluid  dynamics  occur  when  the 
resistive  coupling  to  a  Z-pinch  becomes  competitive  with  the  inductive  coupling.  As  the  re¬ 
sistive  coupling  is  increased,  Ohmic  heating  in  the  outer  regions  of  the  plasma  is  increased. 
If  this  heat  energy  is  not  transported  into  the  inner  regions  of  the  plasma  fast  enough, 
the  energized  surface  simply  blows  away  in  a  1-D  calculation.  The  analysis  of  the  Saturn 
aluminum-wire  shots  presented  in  this  paper  has  suggested  that  L-shell  emissions  have  a 
significant  influence  on  energy  transport  within  the  plasma.  It  is  already  known  experi¬ 
mentally  that  they  have  increasing  effect  on  K-shell  emissions  as  the  atomic  number  of  the 
Z-pinch  is  increased  [31].  However,  as  a  caveat,  kinetic  energy  generation  is  diflferent  in 
one  dimensional  fluid  dynamics  than  in  two-  or  three-dimensions.  Recent  work  on  the  early 
time  dynamics  of  wire  arrays  has  shown  that  plasma  does  not  blow  away  when  there  are 
gaps  between  the  wires.  It  flows  into  the  gaps  instead.  If  similar  two  dimensional  behavior 
persists  late-time,  it  may  play  a  role  in  the  energy  coupling  problem.  In  this  case,  a  gap 
dependent  coupling  should  be  determinable  experimentally. 

Other  issues  need  to  be  addressed  before  a  more  definitive  understanding  of  the  Saturn 
experiments  can  be  worked  out.  More  needs  to  be  learned  about  the  strength  of  surface 
micro-instabilities  and  about  their  influence  on  curent  flow.  If  enhanced  energy  coupling 
is  related  to  the  multi-dimensional  character  of  Z-pinch  implosions,  it  should  be  detectable 
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in  both  wire  array  and  gas  puff  experiments.  On  the  other  hand,  if  it  is  more  prominent 
in  wire  array  than  in  gas  puff  experiments,  it  may  be  related  more  specifically  to  the  way 
wires  break  down  in  time,  generating  a  different  implosion  dynamics  than  is  predicted  by 
0-D  Z-pinch  calculations  [32].  In  all  cases,  however,  the  different  fluid  dynamic  mechanisms 
that  move  energy  about  the  plasma  and  that  aid  in  converting  enhanced  energy  couplings 
into  X  rays  need  to  be  investigated  and  more  quantitatively  understood. 
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FIG.  Jr.  Measured  total  x-ray  output  for  three  identical  shots  having  the  wire  array  paxcimeters 
listed  in  the  figure.  The  time  aixis  has  no  fiducial. 
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FIG.  2.  Measured  K-sheU  x-ray  output  for  the  same  three  identical  shots  as  in  Fig.  1  having  the 
wire  array  parameters  listed  in  the  figure.  The  time  axis  has  no  fiducial. 
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FIG.  3.  The  running  time  integrals  for  the  three  x-ray  pulses  shown  in  Fig.  1  are  plotted  as  a 
function  of  the  integration  time. 
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FIG.  4.  A  circuit  resistance,  which  produces  a  comparable  energy  coupling  as  seen  experimen¬ 
tally,  is  plotted  as  a  function  of  time.  It  consists  of  a  part  derived  from  Eq.  (32),  using  Braginskii’s 
formula  for  electrical  resistivity,  and  another  part  that  is  prescribed.  The  dotted  points  represent 
the  circuit  resistance  that  was  calculated  using  Eq.  (39). 
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FIG.  10.  The  total,  resistive,  and  inductive  energies  coupled  to  the  Z-pinch  load  are  shown  as  a 


function  of  VEiriations  in  the  maximum  value  of  the  Fig.  4  circuit  resistance.  Spline  fits  are  drawn  to 


the  calculated  points.  The  oscillations  in  the  inductive  coupling  at  smaU  values  of  represent 


numerical  error. 
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maximum  value  of  the  Fig.  4  circuit  resistance.  A  spline  fit  is  drawn  to  the  calculated  points. 
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FIG.  13.  The  total  and  resistive  energies  that  couple  to  the  Z-pinch  are  shown  as  a  function  of 
variations  in  the  time  to  peak  resistance  of  the  Fig.  4  circuit  resistance.  Spline  fits  are  drawn  to 
the  calculated  points. 
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FIG.  15.  Four  calculated  resistivities  are  shown  as  functions  of  time  for  the  Figs.  4  and  5  cal¬ 
culation.  The  one,  labeled  rjBrag,  is  derived  from  a±rjo  {rjo  from  Eq.  (32)),  The  dashed  cmrve  is 
derived  from  the  drift-wave  formula,  Eqs.  (43).  The  one,  labeled  »7r=6,  is  derived  from  Rioad  using 
Eqs.  (32)  and  (40),  and  the  one  labeled  »7r=ro  is  linearly  extrapolated  down  from  rjr=b-  ro  is  the 
inner  surface  location  of  the  imploding  plasma  shell. 
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Confinement  and  compression  of  magnetic  flux  by  plasma  shells  is  of  interest  for  a  variety  of 
applications  associated  with  keV  x-ray  production  in  pulsed-power  driven  Z-pinch  plasma  radiation 
sources  (PRS).  Confinement  of  a  field  with  a  Bg  field  of  the  pinch  current  can  help  stabilize  the 
implosion  of  a  shell  firom  a  large  initial  radius.  Compression  of  the  azimuthal  magnetic  flux  in  the 
PRS  (secondary)  circuit  with  a  plasma  shell  driven  by  5  ^  field  of  the  primary  circuit  may  represent 
a  new  opportunity  for  using  a  low-cost,  relatively  slow  pulsed  power  to  generate  large  amounts  of 
keV  X  rays.  The  magnetic  field  has  to  be  compressed  and/or  confined  by  low-mass  plasma  shells 
emitting  soft  x-ray  radiation  that  limits  the  temperature  and  conductivity  of  the  shell  plasma.  The 
thickness  of  a  plasma  shell  is  established  self-consistently  during  the  implosion,  and  it  is  not  obvious 
whether  or  not  the  shell  becomes  thick  enough  to  confine  or  compress  the  magnetic  flux.  The  results 
of  analytical  theory  and  numerical  simulafions  demonstrate  that  the  flux-compressing  capability  of 
a  low-beta  plasma  shell  is  surprisingly  good  because  the  shell  is  shown  to  dynamically  adjust  its 
thickness  so  that  it  always  remains  of  the  order  of  its  skin  depth.  The  self-similar  profiles  of  confined 
magnetic  field  predicted  by  the  theory  are  consistently  reproduced  in  the  simulations.  ©  2003 
American  Institute  of  Physics.  [DOI:  10.1063/1.1614254] 


I.  INTRODUCTION 

Magnetic  flux  compression  (MFC)  as  a  method  for  gen¬ 
erating  high  magnetic  fields  and  currents  has  been  developed 
in  the  late  1950s,  independently  in  the  US*  and  in  the  Soviet 
Union,^  and  is  still  being  advanced  (see  Refs.  3-5  and  ref¬ 
erences  therein,  as  well  as  the  proceedings  of  the  Megagauss 
conferences).  The  concept  is  to  compress  a  closed  conduct¬ 
ing  cavity  so  rapidly  that  the  magnetic  flux  trapped  inside 
cannot  escape  via  diffusion  through  the  conducting  surface. 
Then  the  magnetic  field  inside  the  cavity  grows  in  inverse 
proportion  to  its  shrinking  cross-sectional  area,  due  to  con¬ 
servation  of  the  trapped  magnetic  flux,  and  so  does  the  in¬ 
duced  diamagnetic  current  that  generates  this  field.  In  a  cy¬ 
lindrical  geometry,  we  can  compress  either  the  axial 
magnetic  field  inside  an  annulus,  or  the  azimuthal  magnetic 
field  trapped  between  the  central  conductor  (stator)  and  the 
liner  imploding  onto  it  (armature).  The  former 
configuration — the  MC-l-type  generator,  using  the  term 
originally  introduced  by  Sakharov^ — produces  an  ultrahigh 
pulsed  magnetic  field  (up  to  ~28  MG)  in  the  annulus, 
whereas  the  latter  one,  labeled  the  MC-2-type  generator,  pro¬ 
duces  an  ultrahigh  current  pulse  (up  to  ~  100  MA)  in  a  load 
connected  to  it.  These  generators  are  explosively  driven, 
which  limits  the  characteristic  velocities  of  imploding  solid 
conductors  at  the  level  of  ~5  km/s.  With  initial  radii  about  5 
cm  or  more,  the  explosive  MFC  generators  operate  in  a  mi¬ 
crosecond  range  of  current  and  magnetic  field  rise  times.  The 
same  time  scale  is  characteristic  of  other  concepts  of  MFC 
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generators,®  e.g.,  of  those  compressing  magnetic  flux  with  a 
converging  cylindrical  shock  wave  that  converts  a  dielectric 
medium  into  a  conductor.® 

The  MFC  concept  could  be  advanced  to  generate  (mr- 
rents  and  fields  with  nanosecond-range  rise  times.  To  make  it 
work  the  velocity  of  the  conducting  liner  that  compresses  the 
magnetic  flux  should  be  of  order  of  100  km/s.  There  seems  to 
be  no  way  to  accelerate  solid  conductors  to  ~  100  km/s,  ex¬ 
plosively  or  otherwise.  The  magnetic  flux  could  be  com¬ 
pressed  fast  enough  only  with  a  low-mass  plasma  shell.  Then 
the  velocities  in  the  desired  range,  —100  km/s,  are  well 
within  the  reach:  In  the  multi-MA  Z-pinch  experiments  de¬ 
signed  to  produce  keV  x  rays^  or  soft  x  rays  to  heat  a 
hohlraum,®  imploded  annular  plasmas  are  routinely  acceler¬ 
ated  to  hundreds  of  km/s,  e.g.,  see  Refs.  9.  Hence  the  shell 
velocity  might  not  be  an  issue,  but  its  conductivity  can  be¬ 
come  one,  in  contrast  with  the  case  of  a  solid  conductor.  In 
the  latter  case,  one  can  ensure  that  the  liner  thickness  ex¬ 
ceeds  the  skin  depth,  so  that  the  magnetic  flux  cannot  escape. 
On  the  contrary,  the  in-flight  thickness  of  a  plasma  liner  is 
not  known  in  advance.  Rather,  it  is  determined  self- 
consistently  in  the  course  of  implosion.  With  the  conductiv¬ 
ity  of  imploding  plasma  determined  by  its  temperature, 
which  in  most  cases  of  interest  is  controlled  by  radiative 
losses,  confinement  of  the  compressed  magnetic  flux  is  not 
ensured  in  advance;  the  flux  can  be  lost  through  the  plasma 
shell  or  be  trapped  inside  it. 

A  similar  theoretical  problem  had  been  studied  some 
time  ago  in  connection  with  the  proposal  to  generate  ultra- 
high  magnetic  axial  field  B^  in  a  MC-l-like  configuration, 
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with  a  Z-pinch-driven  gas-puff  plasma  shell  used  as  a  con¬ 
ducting  liner.^®  It  was  noted’*  that  enhanced  flux  losses  char¬ 
acteristic  of  flux  compression  could  have  rendered  this  con¬ 
cept  unworkable.  Subsequent  theoretical  analysis  and  one¬ 
dimensional  (ID)  simulations’^  had  shown  that  efficient 
compression  of  the  axial  magnetic  flux  by  a  plasma  shell  is 
possible  if  the  shell  velocity  is  sufficiently  high.  This  conclu¬ 
sion  was  confirmed  in  experiments.’^'’^  In  particular,  there 
were  some  indications  that  in  the  MFC  experiment  at  7.5  MA 
on  PROTO-II  generator  at  Sandia,  a  seed  field  of  100  kG  was 
compressed  to  about  42  MG.’"*  In  Ref.  15,  theoretical  analy¬ 
sis  and  experimental  observations  of  improved  stability  of 
Z-pinch  implosions  with  external  axial  magnetic  field  were 
reported.  More  recently,  compression  of  axial  magnetic  field 
B,  was  studied  for  the  purpose  of  stabilizing  Z-pinch  implo¬ 
sions,  particularly  those  starting  from  large  initial  radii,  in 
order  to  increase  the  K-shell  radiation  yields  of  plasma  ra¬ 
diation  sources  (PRS).’®  ’^  In  contrast  with  the  goal  of  gen¬ 
erating  ultrahigh  magnetic  field,  the  latter  application  re¬ 
quires  R,  to  be  as  low  as  possible.  It  was  experimentally 
demonstrated”*  that  stabilization  of  a  double-shell  gas-puff 
load  is  achieved  with  a  half  the  initial  field  R^o  conse¬ 
quently,  a  quarter  of  energy  required  to  compress  it)  com¬ 
pared  to  the  case  of  a  single  shell.’*  Not  only  the  uniformity 
of  radiating  pinch  plasmas  was  much  better  in  shots  with  R^, 
but  the  K-shell  yield  of  neon”  and  argon’®  was  shown  to 
increase,  for  neon  by  about  40%. 

The  possibility  of  generating  multi-MA  currents  via  MC- 
2-like  compression  of  azimuthal  magnetic  field  by  an  annular 
plasma  liner  remained  unexplored  until  it  was  discussed  in 
Ref.  19  in  the  context  of  the  authors’  concept  of  a  staged  Z 
pinch,  where  a  plasma  shell  compresses  azimuthal  magnetic 
flux  initially  created  inside  it  by  a  small  seed  current  passing 
through  a  fiber  or  wire  on  axis.  In  late  1990s,  Center 
d’Etudes  de  Gramat  (CEG,  France)  started  an  experimental 
program  of  K-shell  radiation  production  based  on  magnetic 
flux  compression,  the  Sphinx  project^®  This  effort  is  aimed 
at  proving  MFC  to  be  a  viable  alternative  to  more  costly 
conventional  power  amplification  techniques,  pulse  forming 
water  lines  and  plasma  opening  switches.  Proof-of-the- 
principle  experiments  showing  current  amplification  by  MFC 
were  performed  on  “Z”  accelerator  at  Sandia;^’  the  goals 
and  current  status  of  the  Sphinx  project  are  reviewed  in  Ref. 
22. 

In  this  paper,  the  physics  of  flux  confinement  and  com¬ 
pression  by  a  dynamic  low-beta  plasma  is  studied  analyti¬ 
cally  and  numerically.  Our  main  result  is  that  the  diffusive 
losses  of  compressed  magnetic  flux  do  not  constitute  a  seri¬ 
ous  problem:  the  plasma  is  shown  to  adjust  in  such  a  way 
that  its  thickness  is  of  the  order  of  its  skin  depth,  and  the 
magnetic  flux  is  essentially  conserved.  In  Sec.  n  we  present 
analytical  theory.  In  Sec.  EH,  its  results  will  be  compared  to 
numerical  simulations,  and  in  Sec.  IV  we  conclude  with  a 
discussion. 

II.  THEORY 

The  problems  that  we  consider  here  involve  acceleration 
or  deceleration  of  a  thin  plasma  shell  by  a  magnetic  pressure. 


The  best-known  example  is  acceleration  of  an  imploded  an¬ 
nular  shell  to  the  axis  by  the  azimuthal  magnetic  field  Bg, 
converting  magnetic  energy  into  the  kinetic  energy  of  the 
shell.^’*  When  an  axial  magnetic  field  R^  is  present,  which  in 
most  cases  is  small  initially,  but  grows  faster  than  Rg  in  the 
course  of  implosion  (Re“/?~‘  and  B^<^R~^),  the  shell  is 
first  accelerated  by  the  pressure  of  Bg  field  and  then  decel¬ 
erated  by  the  pressure  of  Rj  field,  thus  converting  its  kinetic 
energy  into  magnetic  energy  of  the  compressed  R^  field,’^'” 
as  in  MC-1  generator.  For  a  current-driven  magnetic  flux 
compression  of  MC-2  type,^®”^  a  plasma  shell  (armature)  is 
first  accelerated  by  the  R^  field  produced  in  the  primary  cir¬ 
cuit,  and  then  decelerated  by  the  Bg  field  in  the  secondary 
circuit,  thus  using  its  kinetic  energy  to  drive  the  secondary 
current. 

Let  a  thin  annular  shell  be  initially  located  at  r=Ro  and 
have  a  radial  velocity  vq  towards  a  perfectly  conducting  sta¬ 
tor  column  whose  radius  is  /?, .  The  magnetic  field  in 
vacuum  between  the  shell  and  the  stator  is  produced  by  the 
current  flowing  through  the  stator  and  closing  through  the 
shell,  whose  initial  value  is  Iq  .  The  magnetic  field  profile  in 
vacuum,  at  R,<r<Ro,  is  given  by 


(1) 


where  Bq^ZIqIRoC-  Assuming  both  the  stator  and  the  shell 
to  be  perfectly  conducting,  we  conclude  that  magnetic  flux 
trapped  in  vacuum  between  them  is  conserved.  Its  initial 
value  per  unit  length  is  given  by 

f^O  JRn 

d)=  R(r)Jr=RoRoln^.  (2) 

Due  to  conservation  of  magnetic  flux,  the  right-hand  side  of 
Eq.  (2)  should  yield  the  same  value  of  <I>  for  Rq  replaced 
with  any  other  value  of  the  shell  radius  r=R.  Therefore,  the 
magnetic  field  at  the  surface  of  the  imploding  shell  at  the 
moment  when  the  shell  radius  equals  R,  is  given  by 


R(R)= 


BqRq  ln(Ro/Rj) 
R  ln{R/R,)  ■ 


We  find  the  maximum  values  of  magnetic  field  and  cur¬ 
rent  from  energy  conservation.  The  initial  kinetic  energy  of 
the  shell  is  equal  to  the  work  done  by  the  magnetic  pressure 
as  the  shell  converges  to  its  minimum  radius,  r=R,^„,  where 
it  stagnates: 


1  2  r 


«o  B^(R) 


lirRdR 


BlRl[HRofR,)f  r«o 


4 


R[URIR,)f 


ln(Rmi„/Rj)  ln(Ro/Rj) 


BlRlHRolRs)iB^R^,, 


(here,  fi  is  the  shell  mass  per  unit  length),  so 
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Bo  «min\  B^Bgln(Bo/B,)/’ 

^max _ j  l/jLvl  \n(RolR,) 


(5) 

(6) 


Equation  (6)  presents  the  current  amplification  as  a  ratio 
of  two  logarithms.  To  make  this  ratio  large  enough  (to  be  of 
practical  interest,^®"^^  magnetic  flux  compression  should  am¬ 
plify  the  current  in  the  above  short-circuit  regime  at  least  by 
a  factor  of  5-10),  the  logarithm  in  the  denominator  should  be 
small,  which  implies  Rmn~f^s‘^^ima-  I*  means  that  the  gap 
spacing  between  the  shell  and  the  stator  near  stagnation 
should  be  much  less  than  the  stator  radius. 

The  near-stagnation  phase  of  implosion  is  the  most  rel¬ 
evant  for  current  amplification,  as  well  as  the  most  danger¬ 
ous  with  respect  to  the  magnetic  flux  losses,  since  the  radial 
gradient  of  magnetic  field  also  peaks  at  stagnation.  There¬ 
fore,  in  what  follows  we  limit  ourselves  to  the  radii  R  that 
satisfy  the  inequality  Rq-R'^Rq.  In  this  range,  the  curva¬ 
ture  of  the  stator  and  the  shell  (whose  thickness  is  also  as¬ 
sumed  to  be  much  less  than  its  radius)  is  not  important, 
which  allows  us  to  use  the  thin-shell  approximation  and  a 
planar  geometry.  In  the  case  of  acceleration,  a  plasma  shell 
of  finite  conductivity  whose  thickness  (or  a  characteristic 
length  scale)  is  a{t),  and  initial  velocity  zero,  is  pushed  by  a 
magnetic  field  B(r)  whose  time  dependence  is  known.  In  the 
case  of  magnetic  flux  compression,  a  plasma  shell  ap¬ 
proaches  a  perfectly  conducting  stator,  compressing  the  mag¬ 
netic  field  B(f)  in  the  vacuum  gap  [whose  width  is  A(r)] 
between  the  shell  and  the  stator.  Amplification  of  the  mag¬ 
netic  field  is  due  to  compression  of  magnetic  flux  in  the  gap. 
The  magnetic  flux,  however,  is  not  fully  confined  to  the  gap 
spacing.  Some  of  it  is  trapped  in  the  plasma  shell  or  lost 
through  it  to  infinity.  We  assume  that  the  magnetic  field  is 
zero  behind  the  plasma  shell,  and  the  flux,  having  been  lost 
through  the  shell,  is  instantly  distributed  in  the  infinite  space 
behind  it  without  producing  a  finite  magnetic  field  there. 

The  planar  limit  of  Eqs.  (5)  and  (6)  corresponds  to  Rq 
— /?j= Then  ln(Bo/B,)=Ao/Bo,  the  mass  per  unit 
length  fi  is  expressed  via  the  mass  per  unit  area  M  by  fi 
=  27r/?oA/,  and  Rq/R^„=1.  Equations  (5)  and  (6)  in  this 
limit  are  reduced  to 


ited  by  the  instability  of  the  decelerated  plasma  interface,  we 
can  conclude  that  magnetic  field  compression  ratio  is  ap¬ 
proximately  equal  to  the  initial  kinetic-to-magnetic  energy 
ratio,  and  could  be  made  as  high  as  the  stabihty  limitations 
on  the  radial  compression  ratio  would  permit.  The  question 
we  seek  to  answer  in  this  study  is:  What  are  the  additional 
limitations  on  magnetic  flux  compression  due  to  finite  con¬ 
ductivity  of  the  plasma  shell?  Under  what  conditions  do  they 
become  critical  or  turn  out  to  be  inessential? 

Our  physical  model  is  given  by  the  following  ID  MHD 


equations: 

dp  d 

-+-(p.)  =  0, 

(8) 

du  du  1  dp  1  d 

dt  dx  p  dx  p  dx  BTT 

(9) 

dB  d  d  j  dB\ 

(10) 

These  equations  are  written  in  a  noninertial  reference  frame, 
in  which  the  inner  boundary  of  the  plasma  shell,  which  cor¬ 
responds  to  x=0,  is  at  rest;  the  effective  gravity  acceleration 
in  (9)  is,  of  course,  g=  -d^Aldt^.  The  plasma  is  supposed 
to  be  cold,  with  its  kinetic  pressure  much  less  than  the  mag¬ 
netic  pressure  (low  beta),^  so  that  in  what  follows  only  the 
latter  is  retained  in  the  equation  of  motion  (9).  We  assume 
the  plasma  temperature  to  be  radiation  controlled,  essentially 
constant  and  uniform  throughout  the  bulk  of  the  plasma  dur¬ 
ing  the  implosion,  as  in  most  experimental  conditions  of  in¬ 
terest  for  keV  radiation  production  in  Z-pinch  PRS.^"^^ 
(Note  that  by  neglecting  the  thermal  pressure,  we  overesti¬ 
mate  the  diffusion  of  magnetic  flux  through  the  shell,  that  is, 
underestimate  the  capability  of  the  shell  to  confine  the  flux. 
Indeed,  the  opposite  limiting  case  of  high  p^]=%TTplB^>\ 
in  the  volume  of  the  shell  accelerated  by  magnetic  pressure 
to  supersonic  velocities  corresponds  to  a  skin  depth  much 
less  than  the  shell  thickness,  that  is,  to  a  perfect  confinement 
of  magnetic  flux.  The  ID  simulations — e.g.,  see  Ref.  23 — 
show  that  1  in  the  volume  of  the  plasma  shell  at  all 
time  prior  to  stagnation.)  Therefore  the  plasma  conductivity 
<r  determined  by  its  temperature,  as  well  as  magnetic  diffu- 
sivity  defined  as 


Bnax  ^  47r3ft;g 
Bq  BqAq  Iq 


(7) 


”  direr 


=  const 


(11) 


It  is  easy  to  obtain  Eq.  (7)  directly  by  equating  the  initial 
kinetic  energy  per  unit  area  of  the  shell  Mvl/2  to  the  in¬ 
crease  in  the  energy  of  the  compressed  field  from  the  initial 
value  BqAo/Stt  to  the  peak  value  B^A^nin/Sir  at  the  instant 
when  the  shell  stagnates,  and  using  the  magnetic  flux  con¬ 
servation  equation  appropriate  for  planar  geometry, 
BmaxAmin=BoAo,  where  Amn^Rmio-Rs  IS  the  minimum  gap 
spacing  between  the  shell  and  the  stator.  Formula  (7)  says 
that  the  lower  the  initial  field  Bq  is,  the  higher  peak  field 
Bmax  would  be  produced,  which  is  not  surprising,  because  the 
compression  ratios  required  for  that  must  be  very  high: 
^o^^mo-^nax/^o-  Since  the  compression  ratio  itself  is  lim- 


are  assumed  constant  below.^"* 

Our  initial  and  boundary  conditions  are,  respectively, 

B(0)  =  Bo,  a(0)  =  ao,  A(0)  =  Ao, 

dA\ 

(the  case  of  shell  acceleration  fi-om  zero  velocity  corresponds 
to  Uo=0), 

B=B(f)  at  x=0;  B=0  at  x=a{t)  or  at  x— kx>. 

(13) 


66 


4438  Phys.  Plasmas,  Vol.  10,  No.  11,  November  2003 


Rudakov  et  al. 


The  thin-shell  approximation  assumes  the  time  variation 
of  the  boundary  conditions  to  be  much  slower  than  it  takes 
for  a  fast  magnetosonic  wave  to  travel  across  the  plasma 
shell,  thus  establishing  homogeneous  deformation  through¬ 
out  the  bulk  of  the  plasma  (see  Refs.  12  and  25  and  refer¬ 
ences  therein).  We  therefore  seek  a  solution  of  Eqs.  (8)-(10) 
with  the  initial  and  boundary  conditions  (12)  and  (13)  within 
the  class  of  self-similar  solutions  with  homogeneous  defor¬ 
mation,  as  in  Ref.  12.  Some  particular  solution  of  this  type 
was  first  investigated  in  Ref.  26.  We  choose  the  time  unit  as 

ro=-^(87rMAo)"'2,  (14) 

"0 

introduce  dimensionless  time,  self-similar  coordinate  and 
compression  ratio  as 

T=tltQ,  7}=x/a{r),  a{T)=a(t)/a(„  (15) 

and  use  the  self-similarity  ansatz 
an  da  N{ii) 

(16) 

Here,  the  normalization  constant  po  is  chosen  so  that  the 
total  areal  mass  of  the  plasma  shell  M—p^aQ.  The  dimen¬ 
sionless  density  fiinction  N{t])  should  thereby  satisfy  the 
normalization  condition 

f%(v)d7/=l,  (17) 

Jo 

where  7=  corresponds  to  the  outer  boundary  of  the 
plasma  shell  (i7»=“  if  the  shell  is  diffuse). 

Substituting  (15)  and  (16)  into  (8),  we  find  that  the  con¬ 
tinuity  equation  is  identically  satisfied  with  arbitrary  density 
function  A^(  7).  The  equation  of  motion  (9)  is  reduced  to 


Substitution  of  (19)  into  (18)  yields 


Then  the  induction  equation 
^  aj  Rm  H  drf'' 


(21) 

(10)  transforms  into 

(22) 


where  the  magnetic  Reynolds  number  is  defined  as 

Rm  =  ayvJo.  (23) 

The  left-hand  side  of  Eq.  (22)  depends  only  on  r,  the 
right-hand  side  on  7.  The  variables  in  Eq.  (22)  could  only  be 
separated  if  both  the  left-  and  right-hand  sides  of  this  equa¬ 
tion  are  equal  to  some  separation  constant.  We  denote  this 
constant  depending  on  its  sign.  Three  options  are  im¬ 
mediately  seen. 

(1)  Separation  constant  is  positive,  +  s^.  Then  the  profile 
shape  is  exponential: 

H{T})  =  txp{—  rfSyjRm) 


exp(  2  7„s  yjRm)  —  1 


sinh(7jV^ffi).  (24) 


This  choice  of  constants  satisfies  the  boundary  conditions 
(20).  If  the  plasma  shell  is  diffuse  and  extends  to  infinity 
(7„=«’),  then  only  the  decaying  exponential  remains.  We 
see  that  this  solution  exists  and  is  physically  meaningful 
[corresponds  to  a  positive  Nirt)  in  Eq.  (19)]  for  any  7„ .  In 
this  case  the  magnetic  flux  contained  in  the  plasma 


r  f  Vm 

‘J'pi(7’)  =  j^  Bdx  =  aoBoa(T)/3(r)j^  H{7i)dr) 


where  5(t)  =  A(r)/Ao  is  the  dimensionless  gap  between  the 
shell  and  the  stator,  and  the  dot  denotes  differentiation  with 
respect  to  normalized  time,  t. 

The  thin-shell  approximation  means  neglecting  the  first 
term  in  the  left-hand  side  of  (18)  compared  to  its  right-hand 
side.  It  is  justified,  at  least  for  some  time,  because  of  the 
smallness  of  the  initial  ratio  UqIAo.  For  our  solution  to  be 
applicable,  we  need  to  verify  the  smallness  of  (go/Aq) 
X(arS)  in  the  process  of  compression. 

Ihe  thin-shell  approximation  allows  us  to  separate  the 
variables  in  Eq.  (18)  if  the  7-dependent  term  is  constant.  The 
choice  of  this  constant 


aoBoa{T)^{T)  ( Vm  r;r-\ 

= - P= - tanh  —  j  -^Rm  (25) 

s^Rm  \  2  / 

grows  with  time. 

(2)  Separation  constant  is  zero.  The  magnetic  field  pro¬ 
file  is  linear, 

H{7i)  =  \-T)lTj„.  (26) 

The  shell  thickness  is  finite.  Again,  the  solution  is  physi¬ 
cally  meaningful  for  arbitrary  7„ .  The  magnetic  flux  con¬ 
tained  in  the  shell  is  constant  in  time:  as  much  flux  diffuses 
from  the  inside  as  is  lost  through  the  outside  surface. 

(3)  Separation  constant  is  negative,  —s^.  The  profile 
shape  is  trigonometric 


(19) 


satisfies  the  normalization  condition  (17),  since  the  initial 
and  boundary  conditions  (12),  (13)  for  B  imply  that 

i/(0)=l,  H(7J  =  0.  (20) 

The  solution  is  physically  meaningful  only  if  W(7)^0  for 
all  7,  so  the  derivative  of  magnetic  field  (and  thus,  the  cur¬ 
rent  density)  cannot  change  sign. 


7)  =  cos(  7]s  ^Rm)  —  cot(  rj^s  ■<jRm)sin{  rjs  ^Rm). 

(27) 

The  solution  is  only  physically  meaningful  for  ij^S'jRm 
<7r/2.  The  magnetic  flux  contained  in  the  shell  decreases 
with  time,  since  more  flux  is  lost  through  the  outside  surface 
than  added  via  diffusion  from  the  inside. 

Let  us  consider  the  first  option,  which  obviously  is  the 
most  interesting  physically.  We  derived  two  equations,  (21) 
and 
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$  a 

For  the  particular  case  of  the  exponential  solution  with 
57m we  can  reduce  the  number  of  dimensionless  param¬ 
eters  of  the  system  using  its  symmetry.  Indeed,  the  definition 
of  the  length  scale  a{t)  characteristic  of  the  exponential  pro¬ 
file  allows  multiplying  it  by  any  number  Con¬ 

sequently,  we  have  some  arbitrariness  in  the  choice  of  qq  and 
Rm:  our  exponential  solution  depends  on  xs\[RmlaQa(t) 
=xs/a{t)^lv„to,  a  combination,  which  does  not  include  Oq 
at  all.  Thus  we  can  define  oq  by  postulating 


which  simplifies  the  expressions  for  the  profiles  of  magnetic 
field  and  density, 

//(J7)  =  exp(- 7/),  N{i})  =  2expi-2ij),  (30) 

and  transforms  (28)  into 

$  a  I 

For  the  simple  case  of  a  shell  acceleration  by  a  known 
magnetic  field  B(t)  =  BoP(T)  Eqs.  (21)  and  (31)  are  imme¬ 
diately  integrated.  In  particular. 


where  oq  is  a  positive  integration  constant  When  B(t)  is 
either  constant,  or  linearly  rising  with  time,  or  decreasing  as 
1/r,  Eq.  (32)  implies 

B{t)  =  Boe{t), 

a{t)  =  '  (2v„r/3)‘'2,  B{t)  =  Botlt„ 

,  (2v„r(r-bri)//,)>'2.  B{t)  =  Bot,l{t+t,). 

(33) 

We  have  substituted  here  the  definition  (23)  of  Rm  to 
convert  the  result  into  a  dimensional  form.  Not  surprisingly, 
it  does  not  contain  the  arbitrary  normalization  constant,  oq  • 
The  square-root  time  dependence  in  (33)  is  something  one 
expects  when  dealing  with  diffusion.  Note,  however,  that  the 
magnetic  field  profiles  are  not  Gaussian,  as  typical  for  self- 
similar  diffusion  into  an  incompressible  medium  in  planar 
geometry,^“^  but  rather  exponential! 

This  simple  solution  allows  us  to  estimate  the  role  of  the 
omitted  thermal  pressure  term  in  the  equation  of  motion  (9). 
For  instance,  taking  the  self-similar  solution  (33)  for  constant 
magnetic  pressure,  we  find  that  the  ratio  of  the  omitted  ther¬ 
mal  pressure  term  to  the  retained  magnetic  pressure  terra 
equals  to 


dpidx 

(dldx)B-l%TT 


\()v(\+Z)noToaQ  ^ 

{2^., rat 


ness  of  the  thermal  pressure  compared  to  magnetic  pressure, 
and  additionally  decreases  with  time  due  to  diffusion  [the 
short  time  interval  when  the  term  {2v„t)^^  in  the  denomi¬ 
nator  is  much  smaller  than  oq  is  of  no  importance  because  it 
corresponds  to  the  transient  phase  when  the  self-similar  pro¬ 
files  are  formed,  see  Sec.  IH  B]. 

For  the  case  of  magnetic  flux  compression,  the  time  de¬ 
pendence  of  compressed  magnetic  field  jS(  r)  is  not  known  in 
advance;  rather,  it  should  be  determined  self-consistently,  to¬ 
gether  with  a(r)  and  ^(r).  To  close  the  system  (21)  and  (31), 
we  need  one  more  equation.  It  can  be  derived  from  conser¬ 
vation  of  magnetic  flux.  If,  for  instance,  we  assume  that  no 
flux  losses  take  place  through  the  outer  surface  (which  means 
a  diffuse  shell,  then  we  can  state 

<f>pl+  ‘I>vac=  ao®0«{  T")  +  t) 

=  ^>n= const  (35) 


J?oAo  ^  Aq  ^  ^0,vac\^ 


where  Tp  is  the  constant  temperature  of  the  plasma,  and  Z  is 
the  corresponding  average  ion  charge.  The  ratio  (34)  is  uni¬ 
form  over  the  shell  volume.  It  is  initially  small  due  to  small¬ 


denotes  the  ratio  of  magnetic  fluxes  initially  contained  in  the 
shell  plasma  and  in  vacuum  between  the  shell  and  the  stator, 
respectively.  In  other  words,  the  initial  ratio  of  magnetic 
fluxes  contained  in  the  gap  and  in  the  shell  is  equal  to  the 
ratio  of  their  length  scales  according  to  the  definition  of  Oq 
given  by  (29). 

Thus  we  have  defined  all  the  dimensionless  parameters 
on  which  the  solution  depends  in  the  case  of  magnetic  flux 
compression.  They  are  the  small  ratio  of  initial  magnetic 
energy  in  the  gap  to  the  initial  kinetic  energy  of  the  shell, 

e=2-^<«l,  (38) 

AttMvq 

the  small  ratio  of  initial  shell  thickness  to  the  initial  gap 
spacing, 

^0  '^O.vac 

and  the  magnetic  Reynolds  number  defined  by  (23)  with  the 
length  scale  Qq  consistent  with  (29).  The  magnetic  Reynolds 
number  could  be  either  large  or  small.  Recall  that  its  large 
value  implies  high  conductivity,  whereas  small  Rm  corre¬ 
sponds  to  a  strong  diffusion  of  magnetic  field.  Equation  (36) 
is  presented  as 

l  +  <f> 

S=  — - (pa.  (40) 

The  initial  conditions  for  the  variables  or,  jS,  ^  are 

a(0)  =  )3(0)=5(0)=l.  (41) 
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The  initial  velocity  of  the  shell  is  —  uq.  the  cor¬ 
responding  normalized  combination  being  VqIq/Aq 
=  (8'7rMi;o/®o^o)*^“(2/e)'^.  see  (14)  and  (38).  Therefore, 
^(0)=-(2/s)*^,  which  implies 

\  +  <j) 

«(0)=-^-(2/s)‘^,  (42) 

Eqs.  (21)  and  (31)  are  rewritten  in  the  form 

a-v,  (43) 


P{\  —  avRm) 
a^Rm 


(44) 


l  +  4>-a^^^Rm^ 

(1  -f  <f>—a^<p)a^Rm^ ' 


(45) 


With  initial  conditions  for  a,  p,  and  v  given  by  (41), 
(42),  our  problem  is  well  defined.  In  the  limit  of  perfect 
conductivity,  Rm^^,  we  find  fi-om  (31),  )3=l/a,  as  ex¬ 
pected  from  conservation  of  flux,  S—  a,  and  (43)  and  (45) 
reduce  to  a=a~^.  This  corresponds  to  the  maximum  com¬ 
pressed  magnetic  field  at  the  instant  when  d(r)  =  0  and  the 
plasma  shell  bounces  from  it,  ySmax"!  "•"£  *>  which  exactly 
reproduces  (7),  as  it  should  be. 

Efficiency  of  magnetic  flux  compression  is  demonstrated 
in  Fig.  1  for  Rm  ranging  between  0.1  and  100  with  s  =  0.1, 
^=0.01  and  0.1.  A  perfect  conductor  would  have  amplified 
the  initial  magnetic  field  by  a  factor  ofl  +  8“*  =  ll.  Figure  1 
shows  that  even  for  magnetic  Reynolds  number  Rm  as  low  as 
0.1,  the  peak  compression  still  exceeds  9.  The  ratio  of  the 
skin  depth  defined  as  Aj=(v„ja/dj')^^  to  the  characteristic 
thickness  of  the  plasma  shell  estimated  fi-om  its  density  gra¬ 
dient  as  agor,  is 


aga 


=  [Rm(2/e)'^-^-l]~'^X 


d(0) 


1/2 


a(r)d(r) 


(46) 


The  right-hand  side  of  (46)  remains  of  order  unity  for 
most  of  the  implosion. 


III.  NUMERICAL  SIMULATION 
A.  Formulation  of  the  problem 

The  above  analysis  was  based  on  the  self-similar  solu¬ 
tions  of  Eqs.  (8)-(10),  corresponding  to  certain  boundary 
conditions.  The  self-similar  profiles  of  density  and  magnetic 
field  correspond  to  the  case  when  the  magnetic  flux  is  suc¬ 
cessfully  confined  and  compressed.  These  profiles  are  ex¬ 
pected  to  be  rapidly  established,  which  means  that  the  shell 
acceleration  or  deceleration  distances  are  much  smaller  than 
Aq.  The  self-similar  flow,  representing  an  intermediate  as¬ 
ymptotics  of  a  wider  class  of  solutions,^  is  formed  when  the 
solution  forgets  the  form  of  initial  motion,  but  selects  from 
the  entire  set  of  information  provided  by  the  initial  condi¬ 
tions,  such  as  initial  density  profile,  a  limited  number  of 
parameters  determining  a  self-similar  solution.^  The  theory 
does  not  allow  us  to  model  the  transient  phase  when  the 
self-similar  flow  is  formed  or  to  estimate  the  duration  of  this 
phase.  It  also  does  not  describe  the  situation  near  the  turning 
point,  when  the  magnetic  fields  on  both  sides  of  the  shell  are 


Rm 


FIG.  1.  Peak  compressed  magnetic  field  (a)  and  magnetic  flux  confined  in 
vacuum  between  the  shell  and  the  stator  (b)  vs  magnetic  Reynolds  number, 
calculated  for  two  values  of  the  parameter  0  (39).  Both  and  <pnc  ^re 
normalized  to  their  respective  initial  values  in  vacuum. 


comparable  and  have  a  different  time  dependence.  To  inves¬ 
tigate  these  cases,  a  numerical  simulation  is  needed. 

As  in  Sec.  H,  we  assume  the  temperature  of  the  acceler¬ 
ated  plasma  to  be  radiation  controlled  (the  Joule  heating  is 
compensated  by  radiation  losses)  and  thereby  constant,  suf¬ 
ficiently  low  to  neglect  thermal  pressure  of  the  plasma  in 
comparison  with  the  magnetic  pressure.  Our  basic  equations 
therefore  are  those  of  continuity  and  motion,  and  the  induc¬ 
tion  equations  for  a  two-component  magnetic  field,  which 
are  presented  below  in  Lagrangian  variables  for  an  inertial 
reference  frame 


<?  /l\  du 
dt\pj  If 


(47) 


du  d  ~  . 

--b-(R^R?)=0, 


(48) 
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(49) 

where 


/: 


p{x',t)dx' 


(50) 


is  the  Lagrangian  mass  coordinate,  and  Rm  the  magnetic 
Reynolds  number.  The  system  (47)-(49)  is  further  approxi¬ 
mated  by  an  algebraic  system  of  equations  for  discrete  values 
of  unknown  functions  p,  u,  and  B,.  y  on  a  ID  spatial  grid,  f,-, 
and  solved  numerically.  We  used  a  fully  conservative  finite 
difference  approximation  scheme,  as  described  in  detail  in 
Ref.  31.  Detailed  discussion  of  this  numerical  technique  is 
beyond  the  scope  of  the  present  paper.  We  only  note  that 
implementation  of  such  schemes  was  already  shown  to  be 
efficient  even  for  solution  of  much  more  complicated,  two- 
dimensional  (2D)  radiative  magnetohydrodynamic 
problems.*^ 

We  use  dimensionless  variables  here.  Both  components 
of  magnetic  field  and  density  are  normalized  with  respect  to 
their  initial  values,  and  pg,  respectively.  The 

jt-coordinate  is  normalized  with  respect  to  the  initial  thick¬ 
ness  of  the  shell,  ag ;  the  time  unit  is 


flgVSrrpg 


B 


yO 


(51) 


This  definition  differs  from  (14)  by  a  factor  of 
(ag/Ag)‘'^<Sl.  Consequently,  this  normalization  makes  the 
magnetic  Reynolds  number  defined  by  the  same  Eq.  (23) 
larger  by  a  factor  of  (Ag/ag)‘^>  1. 

The  energy  conservation  equation  follows  from  Eqs. 
(47)-(49).  The  sum  of  kinetic  and  magnetic  energies  is  not 
conserved  because  of  the  Joule  dissipation  of  magnetic  en¬ 
ergy.  flar^{p^lRm)y~[{dByld^'^  +  idBJd^)'^]. 

The  initial  conditions  are  set  as  shown  in  Fig.  2.  The 
computational  area  extends  between  x=Xnj„<0  and  x 
=Xn^>0.  The  plasma  shell  mass  is  distributed  between  x 
=0  and  x=l  (in  Fig.  2 — uniformly,  with  constant  density 
p=  1,  but  any  other  initial  density  profile  in  the  shell  is  pos¬ 
sible).  We  assume  a  presence  of  a  uniform  low-density 
plasma  with  p=Pmio'^l  at  and  1  <x<x^.  The 

layer  between  and  x=Xsu,ot=  -  Ag  represents  a  mo¬ 

tionless,  highly  conductive  electrode,  “the  stator,”  that  con¬ 
fines  the  magnetic  flux  from  the  left  side  of  the  computa¬ 
tional  area,  and  Ag,  as  above,  is  the  initial  distance  between 
the  shell  and  the  stator,  satisfying  the  thin-shell  assumption 
of  Sec.  n,  Ag^flg.  The  corresponding  values  of  the  La¬ 
grangian  mass  coordinate  fare  0  at  x=x^a’  fsuior. 

at  the  inner  and  outer  boundary  of  the  shell,  and 
Mass  of  the  low-density  plasma  in  the  gap  is  taken  to  be 
much  smaller  than  that  of  the  plasma  shell.  Pinin'^ 
and  therefore 

As  in  Sec.  U,  the  assumption  of  constant  temperature 
implies  constant  values  of  magnetic  viscosity  v„  (11)  and 
magnetic  Reynolds  number  Rm  (23).  Ideally,  the  magnetic 
Reynolds  number  is  finite  in  the  plasma  shell,  zero  in  the  gap 
between  the  shell  and  the  stator  and  behind  the  shell,  and 
infinite  in  the  stator.  In  our  numerical  simulation,  we  assume 


Stator  Plasma 

RfTlstator 


p 

i?7WshelI 

rX.-v;-, 

0  u(t) 

pM 

Pmin 

.-.Ot 

•V-  •  ■-•..>745' 

_ W 

1  1 

•^min  -^stator  ”  •  -^max 


I — I  I  I 

0  Stator  &  il  4iax 

FIG.  2.  Schematic  of  the  geometry  of  our  numencal  problem.  Solid  lines  are 
initial  profiles  of  plasma  density  p  (the  gray  area  represents  the  plasma  shell, 
shown  for  a  flat  density  profile)  and  magnetic  Reynolds  number  Rm.  The 
shaded  area  on  the  left  side  is  the  motionless  stator.  Lagrangian  maitcer 
on  the  right  is  the  boundary  of  a  low-density  “background”  plasma  (p 
-Pma)’  which  external  magnetic  field  can  be  prescribed.  Alternatively, 
the  second  stator  could  be  added  between  and  . 


for  these  respective  magnetic  Reynolds  numbers  some  finite 
constant  values  satisfying  the  inequalities  Rm^^Rm^n 
<?Rmsator-  The  value  of  Rm^f„n  is  varied  between  the  limits 
of  weakly  conducting  (/?msheii<l)  and  strongly  conducting 
(R/Wjheii^l)  plasma. 

Now  we  define  our  boundary  conditions.  In  the  motion¬ 
less,  highly  conducting  stator,  we  assume 

u(f,/)  =  0  at  and  i5,.,(f=0,r)=0.  (52) 

The  boundary  conditions  on  the  other  side  of  the  shell 
depend  on  the  problem  addressed.  We  can  have  another  elec¬ 
trode,  like  a  stator  (the  case  of  a  magnetic  cavity  with  two 
walls  and  with  a  bouncing  shell,  see  below),  and  thus  impose 
boundary  conditions  like  (52),  where  ^staror 
placed  with  and  fniax>^taior-  Alternatively,  we  can  ap¬ 
ply  an  external  source,  “a  generator”  of  magnetic  flux  at  the 
right  boundary,  assuming  that 

«(f=?max.0  =  0,  By^=^^,t)  =  b{t),  (53) 

where  h(r)  is  a  known  function.  The  other  component  of  the 
magnetic  field  is  maintained  constant  at  this  boundary: 

equals  its  initial  value,  which  is  either  0  or  1. 


B.  Simulation  results 

Figure  3  shows  the  profiles  of  density  and  magnetic  field 
obtained  numerically  for  a  constant  one-component  magnetic 
field  instantly  applied  to  the  right  side  of  the  computational 
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FIG.  3.  Successive  density  (a)  and  magnetic  field  (b)  profiles  obtained  nu¬ 
merically  for  Rm;^cii=  10,  B^=0.  The  profiles  here  and  below  are  labeled 
by  the  respective  times  expressed  in  units  (51).  The  plasma  shell  is  acceler¬ 
ated  by  a  constant  external  magnetic  field. 

area:  in  Eq.  (53)  b{t)=6{t),  the  Heaviside  step  function. 
The  profiles  are  shown  vs  the  dimensionless  coordinate 

shifted  by  so  that  the  left  boundary  of  the  stator  corre¬ 
sponds  to  ^=0,  as  shown  in  Fig.  2.  The  magnetic  pressure 
accelerates  the  plasma  shell  towards  the  stator,  in  the  nega¬ 
tive  X  direction,  while  the  magnetic  flux  diffuses  through  the 
shell  due  to  its  finite  conductivity.  Here  we  have  chosen 
^stator^  —70  (which  means  Ao/ao—70),  and  a  high  conduc¬ 
tivity  of  the  plasma  shell:  /fwisheii’^  10-  We  have  taken  Sj 
=0.  Here  and  below,  the  choice  of  other  constants  is  arbi¬ 
trary;  we  have  checked  that  for  Pn.in^l0~^.  /?/ngap^l0“^, 
10^  it  does  not  affect  the  numerical  solution. 

The  shell,  accelerated  and  compressed  by  the  pressure  of 
the  magnetic  field  By  fi-om  the  right,  moves  to  the  left.  By 
the  time  f=2,  after  the  shell  has  traveled  about  twice  its 


-4  -3  -2  -1  0  1 


X 

FIG.  4.  Initial  density  profiles  (flat — solid  curve,  same  as  in  Fig.  3,  and 
linear — dashed  curve);  density  (for  both  cases),  velocity  and  magnetic  field 
profiles  (for  flat  initial  density  profile)  at  r=2. 


initial  thickness,  to  x=  -2,  the  theoretically  predicted  expo¬ 
nential  profiles  (which  are  linear  on  the  semilogaiithmic 
plots  of  Fig.  3)  are  seen  to  be  already  established  with  a  good 
accuracy. 

Figure  4  illustrates  that  this  agreement  holds  for  a  vari¬ 
ety  of  initial  profiles.  Here  the  density  profiles  shown  for  / 
=2  in  Fig.  2  are  compared  to  those  obtained  for  the  initial 
linear  density  profile  in  the  shell,  psbeii(-i^.0)  =  2(l -x) 
(dashed  linesX  In  the  latter  case  the  initial  shell  mass  is  un¬ 
changed  and  the  density  decreases  linearly  firom  Pshen(O) 
=  2  to  psheii(l)=0  (if  we  neglect  pmn'^l)-  By  f=2,  density 
profiles  for  both  cases  approach  the  exponential  shape  pre¬ 
dicted  by  the  theory.  The  velocity  profile  shown  in  Fig.  4  is 
close  to  linear,  also  in  agreement  with  the  theory. 

If  the  plasma  shell  were  incompressible  and  perfectly 
conducting,  then  its  zero-dimensional  (OD)  equation  of  mo¬ 
tion  would  be  p(_<f'AJdt^)=B^,  where  is  the  positive 
distance  traveled  by  the  shell.  With  p  =  1  and  S=  1,  its  solu¬ 
tion  describes  a  constant  acceleration:  Aa=f^/2.  In  Fig.  5, 
this  simple  solution  is  compared  to  the  numerically  obtained 
displacement  of  the  right  boundary  of  the  shell,  Ayy .  Figure  5 
illustrates  that,  as  seen  in  the  figure,  after  the  shell  is  formed 
and  the  analytical  profiles  are  well  established,  its  further 
acceleration  proceeds  in  agreement  with  the  OD  prediction. 

Let  us  test  the  theoretical  prediction  (33)  for  the  in-flight 
thickness  of  the  accelerated  shell.  In  our  units,  for  the  case  of 
b{t)—  9{t),  we  obtain 


(55) 


The  shell  thickness  in  the  simulation  could  be  estimated 
from  the  magnetic  field  and  density  profiles  obtained  numeri¬ 
cally  as  the  distance  corresponding  to  an  e-fold  decrease  of 
the  corresponding  variable  from  its  peak  value  at  the  bound¬ 
ary. 
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FIG.  5.  Comparison  of  analytically  and  numerically  calculated  time  varia¬ 
tion  for  the  plasma  sheU  thickness  [a(f)  given  by  (55)  vs  numerical  esti¬ 
mates  of  the  skin  depth  Sg{t)  and  the  shell  thickness  <5^(/)  (56)],  and  the 
distance  traveled  by  the  shell  [^^(r)  is  an  analytical  OD  prediction,  Ajv(r)  is 
the  numerical  displacement  of  the  right  boundary  of  the  shell]. 


The  former  estimate  corresponds  to  the  fonnula  (55)  in  the 
theory;  the  latter,  as  demonstrated  by  Eq.  (30),  to  a  half  of 
this  thickness. 

Figure  5  demonstrates  that  although  the  theoretically 
predicted  exponential  profiles  emerge  in  the  simulations  al¬ 
most  immediately,  the  effective  thickness  of  the  shell,  esti¬ 
mated  as  Sg  otlSp,  approaches  the  corresponding  theoreti¬ 
cal  prediction,  a{t),  somewhat  later,  by  r=4-5. 

Finally,  at  f>5,  after  the  shell  has  propagated  more  than 
10  times  its  initial  thickness,  both  the  profiles  of  density  and 
magnetic  field  and  the  shell  thickness  reproduce  those  found 
analytically.  This  is  demonstrated  by  Fig.  6,  where  the  pro¬ 
files  are  plotted  at  the  time  when  the  accelerated  boundary  of 
the  shell  has  propagated  12  times  its  initial  thickness  (Ajy 
=  12).  For  a  wide  range  of  magnetic  Reynolds  numbers,  the 
density  and  magnetic  field  profiles  are  seen  to  reproduce  the 
theoretical  exponential  shapes.  This  figure  shows  that  even 
for  a  relatively  low  magnetic  Reynolds  number,  Rm^^gn 
=  0.5,  with  initial  thickness  of  the  shell  smaller  than  the  skin 
depth,  the  plasma  shell  still  confines  the  driving  magnetic 
flux  quite  well.  However,  a  small  but  noticeable  part  of  the 
magnetic  flux  is  seen  to  break  through  the  shell  during  its 
formation  and  appears  in  between  the  shell  and  the  stator 
(located  at  JC5tator=  ~  70,  not  seen  in  the  figure).  The  time  t 
=t*  it  takes  for  the  shell  to  reach  x=  — 12  exhibits  a  weak 
dependence  on  the  plasma  conductance  (magnetic  diffusiv- 
ity).  The  difference  in  t*  appears  during  the  shell  initiation, 
when  the  center-of-mass  displacement  still  depends  on  the 
magnetic  field  gradient  value.  Further  on,  the  center-of-mass 
coordinate  varies  as  as  it  should  be.  The  numerical  shell 
thickness  estimated  for  the  conditions  of  Fig.  6  from  Eqs. 
(56)  is  in  all  cases  very  close  to  the  analytical  estimate  (55). 

To  test  the  theoretical  predictions  for  the  waveforms  of 
the  driving  magnetic  field,  which  are  different  from  a  step 
function,  we  choose  /?/nsheii=  10  and  make  the  calculation 
taking  in  (53)  b(t)  =  t/10  (linearly  growing  drive)  and  b{t) 
=  10/(f+ 10)  (decreasing  drive).  In  Fig.  7,  we  compare  these 
two  cases  with  the  above  case  of  a  constant  drive,  b{t) 


(a) 


X 

FIG.  6.  Density  (a)  and  magnetic  field  (b)  profiles  at  the  times  when  the 
right  boundary  of  the  shell  reaches  x-  ~  12  for  the  same  conditions  as  in 
Figs.  3-5  with  the  shell  magnetic  Reynolds  number  varied. 

=  ^f),  for  two  initial  density  profiles  (flat  and  linear,  as  in 
Fig.  4)  at  the  time  when  the  accelerated  boundary  of  the  shell 
reaches  x=  - 12.  For  all  cases,  the  theoretical  exponential 
profiles  are  reproduced.  For  the  cases  of  constant,  linearly 
increasing  and  decreasing  drive,  the  shell  reaches  this  posi¬ 
tion  at  t*  =  5,  11,  and  6,  respectively.  Substituting  these  val¬ 
ues  of  t=t*  into  the  first,  second,  and  third  lines  in  the 
right-hand  side  of  Eq.  (33),  respectively,  we  find  that  the 
theoretically  predicted  shell  thicknesses  should  relate  to  each 
other  for  these  cases  as  1 : \/ll/15:4v(5/5,  that  is,  as 
1:0.86:1.38.  The  ratios  estimated  from  the  numerical  profiles 
in  Fig.  7  are  1:0.8: 1.3,  again,  in  good  agreement  with  the 
theory. 

Figures  4  and  7  demonstrate  that  while  the  shell  adjusts 
itself  to  the  magnetic  energy  flow  through  it,  the  density  and 
magnetic  field  profiles  evolve  to  the  same  self-similar  shape, 
“forgetting”  much  of  the  initial  conditions. We  extend 
this  numerical  study  to  investigate  the  ID  stability  or 
“asymptotic  attractor”  property  of  the  self-similar  solution 
by  testing  whether  it  will  evolve  from  the  initial  density  dis¬ 
tribution  that  is  not  a  single  shell  at  all.  Let  us  start  with  a 
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HG.  7.  Density  (a)  and  magnetic  field  (b)  profiles  at  the  times  when  the 
right  boundary  of  the  shell  reaches  x  =  - 12,  as  in  Fig.  6,  for  /?msiKii=  10  and 
the  waveform  of  the  driving  magnetic  field  varied.  Two  coinciding  curves 
(squares)  for  a  constant  magnetic  field  drive  correspond  to  the  two  different 
initial  density  distributions  of  Fig.  4. 
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FIG.  8.  Density  (a)  and  magnetic  field  (b)  profiles  for  a  plasma  load,  which 
initially  consists  of  a  weakly  conducting  thin  shell  and  a  thicker  low-density 
fill  whose  effective  magnetic  Reynolds  number  is  50  times  higher  (see  pro¬ 
files  at  1=0.1).  By  t=9,  both  components  merge  into  a  single  plasma  shell. 


shell  that  consists  of  two  parts,  roughly  modeling  a  nested 
wire  array  or  a  double-shell  gas  puff  load.®’^^  We  split  the 
plasma  shell  into  two  parts,  placing  the  uniform  outer  shell 
with  Psheii^  1  as  in  Fig.  2  and  adding  a  low-density  fill  with 
uniform  density  pfii]=0.01  between  x=—25  and  x=—20. 
The  initial  density  distribution  is  very  close  to  that  shown  in 
Fig.  8  for  r=0.1.  For  both  the  shell  and  the  fill,  we  take 
This  means  low  conductivity  of  the 
shell,  whose  initial  thickness  and  density  are  unity.  For  the 
fill,  however,  the  dimensionless  parameter  Rm  that  enters 
Eq.  (49)  is  much  less  than  its  effective  magnetic  Reynolds 
number.  Indeed,  combining  the  definitions  (23)  and  (51),  we 
find  that  Rm«.aQ  /  Vp.  The  density  of  the  fill  is  100  times  less 
than  the  shell  density,  and  its  thickness  is  5  times  greater, 
implying  the  value  of  the  effective  magnetic  Reynolds  num¬ 
ber  /?/njiHe£f=0.2X50=  10.  The  shell,  therefore,  initially  is 
relatively  transparent  for  the  magnetic  flux,  whereas  the  fill 
is  not  We  drive  this  system  with  a  constant  magnetic  field. 


b{t)—  0{t),  and  place  the  stator  electrode  very  far  to  the  left, 
4t  XgJgjQi-  700. 

If  the  shell  were  a  solid  conductor,  the  driving  magnetic 
field  would  fi-eely  diffuse,  leak  through  it,  at  f>V8lr 
X/?OTshe)i=  1.  Then  the  shell  would  be  shunted  by  the  plasma 
fill,  and  the  current  switched  to  it.  The  effective  magnetic 
Reynolds  number  for  the  fill  is  large;  hence  the  magnetic 
field  would  accelerate  the  fill  towards  the  stator.  Since  the  fill 
mass  is  only  5%  of  the  shell  mass,  it  would  fly  to  the  stator, 
leaving  the  shell  behind. 

Figure  8  shows  that  it  does  not  happen  this  way.  It  pre¬ 
sents  the  simulated  profiles  of  density  and  magnetic  field  at 
times  varying  from  0.1  to  9.  At  7=0.1  the  magnetic  field 
already  penetrates  through  the  shell  and  reaches  the  fill, 
which  has  a  thickness  greater  than  the  skin  depth.  At  the 
same  time,  the  shell  also  starts  to  expand  with  the  velocity 
equal  to  a  fraction  of  its  center-of-mass  velocity.  As  a  result 
of  this  expansion,  the  heavier  shell  finally  captures  most  of 
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the  current  by  r=5,  and  the  flux  leaking  stops.  Before  t 
=  5,  the  magnetic  field  profile  is  nonmonotonic,  exhibiting  a 
peak  in  the  fill  plasma  and  a  decrease  towards  the  shell, 
where  it  increases  again.  This  indicates  a  current  reversal,  a 
presence  of  a  closed  current  loop  produced  by  the  compres¬ 
sion  of  some  magnetic  flux  diffused  through  the  shell  be¬ 
tween  the  fill  and  the  shell.  Later,  both  parts  of  the  initial 
density  distribution  merge,  forming  a  single  shell  with  a  self¬ 
similar  density  profile  (shown  in  Fig.  8  for  /=9).  We  have 
therefore  demonstrated  a  “ID  stability”  of  the  double- shell 
distribution:  after  some  transient  period,  the  two  components 
coalesce  into  a  single  shell.  This  stability,  demonstrated  for 
the  chosen  value  of  the  parameter  f?msheii=0.2,  is  not  as¬ 
sured  for  an  arbitrarily  small  value  of  this  parameter.  We 
have  found  that  for  a  much  more  diffuse  case,  1, 

there  is  no  stability,  and  the  magnetic  field  diffused  through 
the  shell  plasma  can  accelerate  its  inner  part,  forming  a  pre¬ 
cursor  flow.^  Its  description,  however,  is  beyond  the  scope 
of  the  present  paper. 

We  now  compare  confinement  and  compression  of  mag¬ 
netic  flux  for  the  cases  of  one-  and  two-component  magnetic 
field,  the  latter  referring  to  all  Z-pinch  implosions  involving 

.  The  plasma  shell  is  enclosed  in  a  cavity  bounded  by  two 
conducting  walls.  In  addition  to  the  stator  located,  as  before, 
between  and  x=x^„,  the  second  stator  is  put  on 

the  right  side  of  the  computational  area,  between 
and  x=x^  (we  take  xl^„=29).  This  configuration  does 
not  allow  driving  the  plasma  shell  with  an  external  generator, 
but  we  can  introduce  initial  magnetic  fields  on  both  sides  of 
the  shell  and  follow  their  evolution. 

We  start  with  the  case  when  the  magnetic  fields  initially 
separated  by  the  shell  are  parallel  to  each  other:  0.25, 

jjnghi_  j  •ffjgj.g  jg  ggj(j  inside  the  shell  initially.  The  pres¬ 
sure  of  the  larger  magnetic  field  from  the  right  pushes  the 
shell  to  the  left,  compressing  the  magnetic  flux  of  the  field  to 
the  left  of  the  shell  between  the  shell  and  the  left  stator.  At 
the  same  time,  the  magnetic  field  between  the  shell  and  the 
right  stator  decreases  due  to  expansion  of  magnetic  flux.  The 
increasing  magnetic  field  on  the  left  eventually  exceeds  the 
decreasing  magnetic  field  on  the  right,  starting  to  accelerate 
the  plasma  shell  back,  towards  the  right  stator.  This  gives  rise 
to  oscillations,  which  are  damped  due  to  diffusion  of  mag¬ 
netic  flux  through  the  shell. 

Figure  9  shows  the  successive  density  and  magnetic  field 
profiles  for  Rm^^=2.  At  r=6.7  the  shell  is  accelerated  to¬ 
wards  the  left  (a).  At  r=  10  the  magnetic  fields  equalize  from 
both  sides  of  the  shell  and  the  magnetic  energy  is  converted 
into  the  kinetic  energy  of  the  plasma,  see  also  Fig.  10.  At  t 
=  18  the  shell  stagnates  at  its  extreme  left  position.  The  mag¬ 
netic  flux  is  compressed  (b)  and  the  magnetic  field  and  den¬ 
sity  profiles  are  inversed,  as  was  predicted  analytically.  At  r 
=  26  the  shell  has  bounced,  it  travels  back  to  the  right  and  at 
r=37  comes  close  to  its  initial  position.  The  exponential 
profiles  are  seen  to  form  when  the  magnetic  field  is  much 
higher  on  the  one  side  of  the  shell  than  on  the  other. 

Figure  10  shows  the  time  evolution  of  magnetic  and  ki¬ 
netic  energy  and  magnetic  fluxes  on  both  sides  of  the  shell. 
We  estimate  the  relative  role  of  energy  dissipation  by  calcu¬ 
lating  the  sum  of  magnetic  and  kinetic  energy 


(b) 


FIG.  9.  Successive  density  (a)  and  magnetic  field  (b)  profiles  for  a  plasma 
shell  confined  in  a  conducting  cavity  with  a  one-component  magnetic  field, 
®)e«ll®rlsh(  • 


+£*  (which  in  the  absence  of  Joule  losses  must  be  equal  to 
1  at  all  time,  with  all  the  energies  normalized  with  respect  to 
the  initial  total  energy).  The  deviation  of  this  sum  fi-om  unity 
turns  out  to  be  small,  less  that  10%  after  3.5  full  oscillations, 
at  r=60,  which  demonstrates  that  the  Joule  energy  losses  are 
relatively  low.  The  magnetic  flux  is  seen  to  be  pretty  much 
conserved  on  each  side  of  the  shell,  and  the  damping  of  the 
oscillations  is  very  weak.  This  result  is  favorable  for  the 
one-component  magnetic  flux  compression  by  a  plasma 
shell,  like  the  compression  of  azimuthal  magnetic  flux  in  the 
secondary  circuit  with  the  azimuthal  magnetic  field  of  the 

•  •  •  TO-W 

pnmary  circuit. 

A  different  situation  is  found  for  the  second  case  with 
Bijftl Bright,  as  in  the  experiments  where  the  initial  axial 
magnetic  field  produced  by  a  Helmholtz  coil  was  com¬ 
pressed  with  the  azimuthal  magnetic  field  of  the  pinch 
current  We  keep  the  initial  and  boundary  conditions  the 
same  as  in  Figs.  9  and  10,  and  only  change  the  direction  of 
the  initially  larger  magnetic  field  on  the  right,  £"^*"=0, 
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FIG.  10.  Time  evolution  of  magnetic  fluxes  4>rt8hf  magnetic  E„  and 
kinetic  Et  characteristic  energies  in  the  conditions  of  Fig.  9.  Both  the  fluxes 
and  the  energies  are  normalized  to  the  sum  of  their  respective  initial  values 
in  the  cavity. 

gnght_  2  11  shows  the  evolution  of  density  and  mag¬ 

netic  field  profiles.  Near  the  time  of  bouncing,  when  either 
^nght<^^ieft  Qj.  ^nght^^ieft  ^^—42)^  the  profiles  are 

still  exponential,  as  they  should  be.  The  profiles  are  seen  to 
be  substantially  different,  however,  when  both  components 
of  the  magnetic  field  are  comparable  (r  =  11  and  f=27). 

Besides,  Fig.  12  shows  that  about  one-third  of  the  initial 
left  magnetic  flux  is  lost  from  the  left  side  of  the  cavity 
through  diffusion  during  the  first  deceleration,  by  /=19, 
when  the  shell  compresses  By  flux  in  the  left  semicayity, 
being  driven  by  the  pressure  of  Bj  field  fi-om  the  right,  in 
comparison  to'  the  case  of  parallel  compressed  and  driving 
fields,  where  the  losses  of  the  left  By  flux  during  the  first 
deceleration  are  almost  negligible.  Figure  12  also  demon¬ 
strates  that  a  two-component  magnetic  flux  compression  is 
not  efficient  even  if  the  magnetic  Reynolds  number  of  the 
plasma  shell  is  increased  from  5  to  10.  In  a  sense,  this  is 
consistent  with  the  structure  of  inductance  equation  (21), 
which  would  have  the  same  form  for  each  of  the  two  com¬ 
ponents  of  magnetic  field,  see  (49).  None  of  its  solutions  can 
satisfy  B=0  and  VB=0  at  the  same  time,  ensuring  a  zero 
flux  of  magnetic  energy  through  the  plasma.  Unless  the  ex¬ 
ponential  solution,  with  B,  VB—*0  at  ?;— ><»  is  valid  for  a 
sufficiently  thick  target  (high  effective  value  of  Rm),  as  in 
most  cases  described  above,  fast  losses  of  magnetic  flux  are 
inevitable. 

IV.  CONCLUSIONS 

We  have  demonstrated  that  a  magnetically  accelerated 
low-beta  plasma  sheU  behaves  as  a  self-regulating  structure, 
which  adjusts  its  thickness  to  the  acceleration  conditions  so 
that  the  magnetic  flux  is  confined  even  for  relatively  low 
values  of  magnetic  Reynolds  numbers.  The  case  of  one- 
component  magnetic  field  exactly  corresponds  to  the  use  of 
the  pressure  of  the  azimuthal  magnetic  field  in  the  primary 
circuit  in  order  to  compress  the  azimuthal  magnetic  flux,  and 
thereby  amplify  the  current  in  the  secondary  circuiL^®"^^  For 
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FIG.  11.  Successive  density  (a)  and  magnetic  field  (b)  profiles  for  a  plasma 
shell  confined  in  a  conducting  cavity  with  a  two-component  magnetic  field, 

®righl  • 

this  case,  the  profiles  of  density  and  magnetic  field  were 
shown  to  approach  the  asymptotic  self-similar  exponential 
profiles  predicted  by  the  theory.  Transition  to  the  same 
asymptotic  profiles  was  demonstrated  for  plasma  shells  start¬ 
ing  from  a  wide  variety  of  initial  conditions,  including  a 
plasma  shell,  which  is  initially  split  into  two  nested  sub- 
shells.  Confinement  of  magnetic  flux  was  found  to  be  less 
effective  for  the  case  of  a  two-component  magnetic  field, 
which  corresponds  to  the  Z-pinch  implosions  with  Bj ,  as  in 
Refs.  13-17. 

The  analytical  and  numerical  problems  addressed  in  this 
paper  are  highly  idealized.  Planar  geometry  and  a  consider¬ 
ably  simplified  physical  model  permitted  us  to  find  exact 
analytical  solutions  of  some  problems  under  study,  and  al¬ 
lowed  us  to  use  a  relatively  simple  and  robust  numerical 
approximation  technique  for  testing  the  theory.  This  is,  in  our 
opinion,  a  necessary  step  in  the  study  of  acceleration  of 
plasma  shells  for  the  purpose  of  radiation  production  with 
magnetic  field  geometry  different  fixrm  that  of  a  classical 
Z-pinch.  Even  a  realistic,  systematic  ID  radiation-hydro- 
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FIG.  12.  Time  evolution  of  mngnetic  flux  in  the  left  semicavity,  for 
the  conditions  of  Fig.  II  at  three  values  of  magnetic  Reynolds 

number.  Confinement  of  the  magnetic  flux  is  not  good  even  for  the  highest 
value  of  /?mjhdi=10. 

analysis  of  Z-pinch  implosions  with  ,  has  never,  to  our 
knowledge,  been  performed,  to  say  nothing  about  adequate 
2D  modeling.  A  similar  analysis  is  timely  for  the  future  ex¬ 
periments  on  keV  radiation  production  based  on  the  MFC 
approach.^®*^^  Basic  physics  established  in  the  present  work 
will  be  helpful  for  advancing  in  this  direction. 
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Chapter  5 

Two  and  Three  Wire  Loads  for 
Large  Current  Machines 


5.1  Introduction 

The  use  of  a  single  wire  load  on  large  machines  like  ”Z”  is  usually  avoided 
due  an  excessive  initial  inductance  and  subsequent  high  voltages  on  upstream 
components.  There  is  however  a  middle  ground  where  a  few  larger  wires 
can  start  well  off  axis  and  present  much  less  initial  inductance.  Such  loads 
will  provide  the  lower  inductance  at  the  expense  of  more  precursor  plasma 
involvement. 

In  contrast  to  closed  arrays  with  hundreds  of  fine  wires,  this  relatively 
unexplored  path  to  energetic  implosions  would  use  two  or  perhaps  three 
load  wires  of  appropriately  heavier  mass  and  aims  to  create  a  focused  axial 
collision  of  dense  wire  cores  amidst  the  assembled  precursor  plasma.  The 
precursor,  not  confined  by  a  cage  of  many  wires,  would  presumably  not  much 
soften  this  collision. 

Hence  we  must  examine  the  trade  between  the  implosion  mass  lost  to 
precursor  ablation  and  the  mass  available  to  deliver  the  wire  core’s  kinetic 
energy. 

A  clear  consequence  of  this  load  choice  is  a  quite  open  geometry  for  which 
the  transition  to  a  highly  conducting  annular  MHD  plasma  is  neither  an 
early  nor  a  necessarily  dominant  feature  of  the  electrodynamics.  Generalizing 
slightly  beyond  MHD,  we  use  a  Lorentz  gauge  direct  field  solver  to  treat 
the  TEM  to  TM  mode  set  transition  in  the  pinch  region.  The  problem 
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is  formulated  using  a  scalar  potential  $  and  an  axial  vector  potential  as 
unknowns,  with  the  radial  vector  potential  being  determined  by  the  gauge 
constraint. 

In  contrast  to  earlier  wire  dynamic  model  (WDM)  formulations  with  in¬ 
ductance  matrix  elements  good  only  in  the  thin  wire  limit^,  the  present  work 
makes  use  of  a  new  analytic  Green’s  function  for  A^  that  accounts  for  prox¬ 
imity  effects  among  the  wire  cores  and  the  return  current  structure. 

The  Green’s  function  provides  a  first  estimate  for  the  axial  vector  po¬ 
tential  solution  A2(a;,y)  at  the  pinch  midplane  and  the  net  load  inductance 
arising  among  the  heavy  wire  filaments.  That  solution  is  further  refined 
by  a  current  density  source  term  representing  the  extended  coronal  plasma 
conductivity  and  flow. 

The  use  of  empirical  mass  source  terms^,  for  estimating  the  precursor 
mass  ejection  rate,  ablative  cooling  rate,  and  ablation  velocity,  is  then  intro¬ 
duced  to  complete  the  picture. 


5.2  Mathematical  Formulation 

There  are  two  required  innovations  to  carry  out  the  desired  calculations. 
First,  we  must  go  beyond  the  thin  wire  inductance  limit,  which  is  fine  for 
wire  core  dynamics  and  current  sharing,  but  starts  to  break  down  when  the 
tenuous  corona  plasma  jackets  the  wires  and  entrains  current.  Moreover, 
in  the  larger  wires  used  for  the  loads  we  consider,  proximity  effects  become 
noticeable  as  the  wires  approach  one  another.  Hence  a  means  of  describing 
the  extended  conduction  medium  of  the  corona  while  maintaining  contact 
with  the  lowest  order  inductive  picture  enforced  by  the  wire  cores  is  key  to 
a  versatile  picture. 

Second,  we  require  a  fluid  description  that  admits  a  large  adaptability 
in  scale  lengths  in  order  to  track  the  history  of  coronal  plasma  elements 
that  drift  radially  to  form  the  axial  precursor.  A  smooth  generalization  of 
the  discrete  wire  filaments  into  fluid  particles  allows  a  large  dynamic  range  in 
mesh  size  for  the  field  solver.  The  new  fluid  particle  methods  used  here  adapt 
easily  to  such  grids  because  the  projection  of  the  fluid  particle  variables  like 
density,  velocity,  and  energy  onto  the  Eulerian  grid  is  virtually  exact.  The 
resulting  errors  in  working  with  the  tabular  collisionaJ  radiative  equilibrium 
(TORE)  algorithms  will  also  be  minimized. 
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5.2.1  Energy  Transfer  in  a  Circuit  Model 

For  the  original  WDM  the  direct  interaction  among  an  ensemble  of  current 
elements  is  equivalent  to  a  potential  energy  defined  by  the  array  inductance. 
Within  the  ensemble  are  symmetry  groups  of  multiplicity  S,  e.g. 

^  =  =  (5.1) 

where  a,-  is  any  ^-wise  invariant  current  fraction  contained  in  any  particular 
wire  path  represented  by  a  series  resistance  and  inductance. 

By  construction  then  the  sum  of  all  a{  is  one.  Since  the  wire  paths  are  in 
parallel,  the  voltage  seen  by  each  path  a,-  is  equal  in  the  absence  of  (small) 
wave  transit  time  effects.  If  we  examine  the  voltage  Ve{t)  impressed  at  the 
entrance  or  feed  of  a  wire  array  cage,  then  from  the  following  figure,  the 
resolution  of  the  set  a,-  and  the  total  current  I  admitted  to  the  array  are 
seen  to  be  two  distinct  problems.  From  the  view  of  the  external  generator 
Vg{t),  the  array  is  a  two  terminal  device  and  the  energy  into  it  must  be  the 
same  for  all  sets  a,-  that  show  the  same  net  impedance. 


Figure  1.  WDM  circuit  elements  can  decompose  into  symmetry  groups, 
with  each  element  representing  a  resistance  and 

As  illustrated  in  Fig.  1,  maintaining  equal  voltage  at  each  circuit  element, 
the  equations  fixing  the  current  fractions  a,-  for  any  given  Ve{t)  =  Vg{t)  — 
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Zgsnl  and  I  can  be  written  as, 


Ve(0  =  Ri{t)ail  + 


(5.2) 


and  then  integrated  over  a  time  interval  5  =  eliminate  the  time 

derivative.  Denote  by  A’^e  the  change  in  flux  over  this  interval,  and  the 
matrix  relation 


=  SRjlij  •  aj  +  Lij  ■  Qj  (5.3) 

emerges,  with  I,j  the  diagonal  identity  matrix,  u,-  a  unit  vector  over  the  local 
group. 

In  the  limit  of  vanishing  inductance,  or  slow  timescales  (5  ':$>  Li/Ri),  it 
is  easy  to  see  that  the  a,-  =  Z^\/Ri(i)  and  thus  they  clearly  add  up  to  unity. 


5.2.2  Field  Methods  for  Open  Geometries 

Loads  with  fewer  wires  using  larger  initial  wire  diameters  and  load  radii  will 
present  a  lower  initial  inductance  at  the  expense  of  more  precursor  plasma 
involvement.  In  contrast  to  closed  arrays  with  hundreds  of  fine  wires,  this 
very  open  field  geometry  is  positioned  on  the  edge  of  validity  for  a  conven¬ 
tional  inductance  and  MHD  picture  due  to  the  large  dynamic  range  in  scale 
lengths  and  magnetic  Reynolds  number.  The  limit  of  spatially  constant  volt¬ 
age  throughout  the  corona  is  probably  pretty  good  but  by  no  means  certain. 


Mode  TVansitions  from  TEM  to  TM 

With  Zo  =  yjtolyio  =  367.70,  an  inlet  TEM  boundary  condition  on  the  scalar 
potential  $(r,  z,t) 


Kn  =  Zl,h  = 


I^M-) 

27r 


4, 


(5.4) 


determines  all  the  fields  near  the  inlet,  where  A  ->  0. 

Elsewhere,  the  time  integral  of  the  scalar  potential  forms  a  useful  gen¬ 
eralization  of  the  familiar  inductive  “flux  function”,  viz.  let  ^(r,  2:,t)  = 
P  and 
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'f(r’>)  -  ^(r<)  =  /  Ar{zc,ri,t)dri  ,  (5.5) 

Jr< 

^(z>)  -  ^(2:<)  =  /  A^{zi,rc,t)dzi  ,  (5.6) 

Jz< 

on  radial  Zc  and  axial  conductor  boundaries  to  ensure  that  the  tangential 
component  of  electric  field  vanishes. 


y 


—  1.971  nH/cm 


Figure  2.  Three  1000  fim  Wire  Vector  Potential  Az(x,y) 

Utility  of  the  Lorentz  Gauge 

The  Lorentz  gauge  condition  choice 

r~^dr{rAr)  +  dzA^  +  c~^dt^  =  0  ,  (5.7) 

allows  the  dynamics  to  concentrate  on  the  solution  of  Az{z,r,t)  only.  One 
may  in  fact  examine  only  two  coupled  wave  equations  for  the  needed  fields. 

=  0,  (5.8) 
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with  the  axial  current  density  near  the  wires  given  by 

J.  =  ^tjl[E,  +  V,  xB,(A)]  . 


(5.10) 


For  the  inhomogeneous  wave  equation  involving  Ai(x,  y),  a  Green’s  function 
can  be  developed  for  the  problem  at  the  diode  midplane  which  captures  the 
proper  boundary  conditions  on  the  outer  return  current  surface  and  on  the 
finite  size  conducting  wire  cores.  The  function  can  be  built  up  piecewise  from 
image  currents  generated  by  a  source  cylinder  (a),  any  number  of  floating 
cylinders  b,b’,. . .  and  the  return  cylinder  (c), 

0{x,  y\s  =  G,  6, 6' . . .)  =  ga{x,  y)  +  pi,(x,  y)-|- 

9b'{x,y)-\- gxh{x,y)-\- gxh'{x,y)-\- ...  ,  (5.11) 

and  over  the  symmetry  groups  Az{x,y)  =  Es^(x,t/|s).  This  Aj  “diagonal¬ 
izes”  the  inductance  matrix. 

An  example  is  shown  in  Fig.  2  above  with  quite  large  wires  close  to  the 
return  current  can  where  the  thin  wire  inductance  would  be  inaccurate. 


Variational  Form  of  the  Field  Problem 

The  dimensionless  spatial  variables  are  based  on  a  characteristic  time  tg  of 
the  input  pulse  line,  Xg  —  ctg,  and  a  skin  depth  Sg  =  \r^^. 


With  Ez  =  toE',  action  integrals  advance  all  fields, 


£($)  =  J  drdA  [(V$)"  -  [dr^f]  , 


C{Az)  =  /  drdA  {VAzf  -  {dr 


-{drAzf-Az{-ffez\  . 


(5.12) 


(5.13) 


The  finite  element  (FE)  representation  employs  value  and  time  derivative 
degrees  of  freedom,  e.g.  $  =  [$,5,-#].  The  variation 


.+1  W.+i  .’+*  $  +.+1  w.- #  +.+1  w._i  $  =  0  , 


(5.14) 
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is  equivalent  to  solving  the  PDE  problem  (Eqn.  12)  .  The  symmetric  positive 
definite  weight  matrices  have  the  general  form,  jW„i  =  Tim  ®  {Gx  +  Gy)\l  — 
Tim  <8>  V  with  integrals  of  squared  time  values  T,  time  derivatives  T',  spatial 
gradients  G,  and  spatial  values  V  taken  over  each  cell.  A  similar  formulation 
resolves  Eqn.  13,  here  A  = 

,+iW.-+i  A  +.+1  W.-  i  +H1  W._x  A  =  S-  £{A,  $)  ,  (5.15) 

with  <S  =  {^)Tim  ®  y  and  5  =  to{-{^^drA,  +  +  V,  x  (V  x  A)). 

With  the  equation  for  $  homogeneous,  the  potential  variable  can  actually  be 
replaced  by  is  axial  derivative,  the  static  electric  field  component,  —dz^. 


5.2.3  Particle  Methods 


An  early  formulation^  of  the  fluid  particles  we  want  used  freely  drifting  par¬ 
ticles  only  and  focused  on  particle  annihilation  and  creation  rules  to  track 
emerging  features. 

The  novel  particle  solutions  used  here  are  drift  kinetic  fluid  particles 
(DKFP)  and  track  precisely  the  distribution  function  evolving  under  the 
action  of  the  drift  velocity  C  and  acceleration  A  which  are  functions  of  local 
position  X,  C  =  V  +  X  •  S/h  ,a  =  A  +  x  •  a/h  . 

These  DKFP  contain  three  common  factors  which  represent  the  dilation 
of  the  initial  size  (due  to  the  shear  in  the  velocity  and  acceleration),  and  the 
(asymmetric)  movement  of  the  initial  domain  boundaries.  For  a  ID  particle 
class  these  are: 

=  +  T  + 

hi(t)  =  ±h  +  (V+±5)t+l(A  +  ±a)e  .  (5.17) 

The  expected  profile  for  number  density  (per  unit  length,  area,  or  volume) 
n{X,i)  =<  N  >  IS  then: 


n{X,t)  = 


N 


2hD(t) 


with  similar  expressions  for  momentum  and  enthalpy.  The  fluid  proper¬ 
ties,  needed  by  the  electrodynamics  to  support  the  evaluation  of  a  magnetic 
Reynolds  number  {IZm  =  ^^BioVrl^o)  I'he  source  term  above,  are  then 
easily  projected  into  each  node  or  cell  of  the  field  solution  grid  with  virtually 
no  error. 
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Source  Boundary  Conditions  at  the  Wire  Core 

The  detailed  accounting  of  material  phase  changes  as  the  wires  vaporize  and 
ionize  is  left  to  future  refinements.  A  reasonably  simple  and  apparently 
accurate  approximation  is  to  balance  the  force  of  the  wire  core  attraction 
with  a  “rocket”  force  due  to  the  inward  flow  of  coronal  plasma^.  Here  a 
strict  balance  is  too  specialized,  for  what  is  needed  is  a  path  to  account 
for  wire  loads  that  do  not  hang  at  the  initial  radius  until  they  virtually 
disintegrate.  From  earlier  work^  we  know  that  A1  foil  switches  exhibit  an 
areal  mass  loss  rate  that  scales  directly  with  the  surface  magnetic  energy 
density,  viz.  dtm  —  —  ,  and  a  flow  velocity  determined  by  a  sound 

speed^  or,  as  seen  for  switch  foils,  an  Alfven  speed  characterized  by  the  local 
field.  Here  the  “private  flux”  near  the  wire  plays  the  same  role  and  this  leads 
to  a  mass  loss  estimate,  for  a  wire  of  length  i  and  radius  r^^re,  that  scales 
like 

m  «  \^K{t|r^ire)IlMA][^J■9|ns]  .  (5.19) 

Typical  values  for  «  are  [5-25]  iJ,gcm/ergs  and  we  find  that,  in  keeping 
with  the  other  models  cited  here,  only  a  small  fraction  of  the  wire  current 
(«  10"^)  can  be  viewed  as  directly  participating  in  any  core  erosion  process. 

IcM  Pressure  EntTalpy  and  Resistivity  for  Titanium 
linear  r|(H)  implies  progressive  wave  solutions 
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The  larger  current  fraction  must  be  concentrated  outside  the  core  in  the 
vapor  and  plasma  emitted  by  the  wire.  If  in  fact  the  material  outside  the 
wire  tends  to  a  state  such  the  magnetic  Reynolds  number  is  about  one,  then 
a  similar  analysis  reveals  that 

rh  ^  {llA\MZ\nKas.{Z .  (5.20) 

The  expected  resistivity  of  the  wire  core  can  be  modeled  as  a  function  of 
absorbed  enthalpy  H[x)  as  shown  above  in  Fig.  3.  At  the  larger  temperature 
values  on  this  plot  we  join  the  TORE  model  table  for  radiation  and  transport 
functions,  using  a  magnetoplasma  conductivity  rather  than  •q{H{x)). 

Magnetoplasma  Conductivity  and  Effective  Resistance 

The  Epperlein-Haines®  formulation  for  the  required  skin  depth  So  in  the  wire 
corona  can  be  written 


l/(r±  =  f?x  =  »?o  ax(Z,coT) 
with,  ax(Z,ojT)  a  rational  function,  and 


s]. 


(5.21) 


(5.22) 


Using  the  resistivity  of  Fig.  3  for  the  core  and  this  formulation  for  the  ablative 
corona,  the  parallel  combination  of  resistances  becomes 


'  VA.(xtcRb)  ' 

and  deteremines  the  ohmic  heating  rate  for  the  core  material  consistent  with 
the  mass  loss  model. 

As  needed  in  the  source  term  for  the  coronal  plasma. 


=  1.218472  •  10"^ 

Oo 


\2  rpZl2 

2folnA  Z  a±(Z,Li)T) 


is  the  (dimensionless)  equivalent,  for  a  fiducial  time  tg. 


(5.24) 
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5.3  Early  Inductance  And  Stagnation  Energy 

First  we  have  examined  the  question  of  initial  inductances  for  loads  that 
might  be  contemplated  on  larger  machines  like  “Z”.  For  1,2,  and  3  wire  loads 
set  at  0.9  of  the  return  current  radius,  one  sees  quite  favorable  inductances 
for  diameters  near  or  above  100/rm. 

In  the  following  Table  1,  the  subscripts  on  the  inductance  denote  the 
number  of  wires  and  the  larger  masses  implied  for  diameters  above  100pm 
clearly  would  not  be  viable  implosion  loads,  but  here  hollow  loads  of  large 
radius  and  lower  mass  should  be  examined.  The  advantages  of  larger  radius 
will  also  be  reflected  in  the  mass  ablation  and  core  cooling  rates. 

Clearly  once  two  or  three  wires  are  used  the  initial  inductance  is  quite 
modest  and,  as  expected,  this  presents  no  noticeable  voltage  load  for  the  run 
down  phase  of  the  implosion. 


Table  1.  Initial  Inductances  and  Masses  for  Ti  Loads 


Dia. 

pm 

Li 

nK 

cm 

mPBrn 

Lz 

hK 

cm 

Mass 

cm. 

1000 

5.88 

2.95 

1.97 

500 

7.27 

3.64 

2.43 

.  .  . 

200 

9.11 

4.56 

3.05 

1,426 

100 

10.49 

5.25 

3.51 

357 

50 

11.98 

5.95 

3.97 

89.2 

20 

13.71 

6.86 

4.58 

14.26 

10 

15.10 

7.55 

5.04 

3.57 

With  a  load  mass  m,  in  pgm  per  wire,  inductance  L  in  nH,  driver  voltage 
V  in  MV,  and  pinch  dimensions  of  length  and  radius  tq  in  cm,  the  available 
wire  core  energy  at  stagnation  scales  as, 

ATimp, 2  =  522.4  (ro  (m,T)'/2  v/l)  [fcj],  (5.25) 

One  expects,  Kimp,3  =  \/4/3Kimp,2  for  the  three  wires. 

As  shown  in  Table  2,  in  an  initial  survey  of  the  expected  behavior  for  these 
loads,  the  model  was  used  in  the  limiting  case  of  (small)  fixed  wire  resistance 
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and  a  fixed  mass  erosion  constant  so  that  for  given  current  histories  smaller 
diameter  wires  would  lose  a  higher  mass  fraction. 


Table  2.  Energy  Transfers  for  Ti  Loads 


Dia.  &  Case 
^m  —  Z,DE 

H2 

kJ 

Hs 

kJ 

Mass 

1000  —  Z 

428.2 

591.6 

1,125 

552.6 

802.2 

1,125 

hbh 

631.8 

1,017 

1,125 

HnQBanil 

662.4 

1,004 

1,125 

100  — DE 

78.8 

94.47 

180 

The  2  and  3  wire  loads  were  set  at  0.8  of  the  return  current  radius,  and  run 
for  a  fixed  mass  erosion  coefficient.  On  Z  the  optimizing  3  wire  case  picked 
the  best  final  mass  at  implosion  for  the  200  case.  The  lower  masses  lost 
available  stagnation  kinetic  energy  due  to  mass  ablation,  the  larger  ones  lost 
it  due  to  poor  machine  coupling.  On  DE  the  wires  were  placed  at  0.6  of  the 
return  current  radius. 

With  a  good  load  on  “Z”,  Fig.  4(a,b)  shows  the  radiation  yield  and  power 
accessible  in  a  high  density  stagnation,  using  the  further  features  of  TORE 
radiation  loses^  and  an  explicit  calculation  of  the  stagnation  heating  with 
viscous  drag  forces  slowing  the  model  wire  cores.  These  very  prompt  and 
large  heating  rates  are  seen  to  extract  over  half  the  available  energy  (0.57) 
and  show  net  cooling  rates  clearly  comparable  to  expected  stagnation  heating 
rates.  The  achieved  temperatures  are  also  quite  high  enough  to  extract  some 
K-shell  radiation  from  the  Ti  load. 


5.4  Conclusion 

We  have  (i)  modified  the  wire  dynamic  model  (WDM)  to  the  calculation 
of  kinetic  energy  transfer  to  such  low  wire  number  loads  on  both  Double 
Eagle  and  ”Z”  class  drivers,  and  (ii)  assessed  the  load  energy  and  upstream 
voltages. 

For  DE  any  initial  mounting  of  two  or  three  50  to  200  /im  Ti  wire  loads 
at  radii  in  excess  of  0.6  of  the  return  current  radius  will  easily  preclude  early 


voltages  near  the  insulator  stack  from  rising  above  one  half  the  open  circuit 
voltage  cis  the  load  initiates  and  starts  to  run  in.  The  best  available  kinetic 
energies  are  about  80  kJ  for  two  wire  loads,  and  about  100  kJ  for  three  wire 
loads.  The  (60-80  /zg/cm)  masses  required  are  typical  for  DE,  the  inductive 
current  “bite”  is  clear. 


1.70E*2  1.80E»2  1.90E->2  2.00E*2  2.10E-»2  2.20E*2  2.30E<2 
Time  (nsl 

%oad  Poynting  Kinetic  &iergy  Radiation 


Vport 


E>ithalpy(H)  Mag.  Field  Energy 


Tise  (ns] 


I 

3 


I 


For  “Z”  again  initial  mounting  of  similar  wire  loads  at  radii  in  excess 
of  0.6  of  the  return  current  radius  will  also  keep  early  voltages  in  bounds. 
The  best  available  kinetic  energies  are  about  600  kJ  for  two  wire  loads,  and 
about  100  kJ  for  three  wire  loads.  The  (250-375  /ug/cm)  masses  used  are  an 
excellent  match  in  “Z”.  In  both  cases  the  result  of  too  much  mass  erosion  is 
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to  degrade  the  available  kinetic  energy  at  stagnation. 
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Chapter  6 


Flashover  and  Energy  Coupling 
in  the  DQ  Machine 

6.1  Observations  and  Analysis 

The  last  experimental  campaign  on  Decade  Quad  comprised  a  set  of  7  gas 
puff  shots,  preceded  by  some  wire  calibration  shots.  For  the  gas  puff  series  it 
would  appear  as  if  the  insulator  stack  flashed  over  on  all  of  them,  although 
two  of  the  shots  had  no  “tube  current”  monitor  working  to  confirm  the 
event.  When  the  tube  current  monitor  was  operating,  all  these  events  show 
a  marked  surge  in  the  tube  current  and  usually  show  an  associated  weakening 
of  the  downstream  MITL  and  load  currents.  On  one  shot,  548,  where  both 
front  and  rear  insulators  appeared  to  flash  and  the  tube  current  surge  was 
the  highest,  the  MITL  and  load  currents  appear  to  rise  in  late  time.  Such 
behavior  could  arise  with  a  load  plasma  rebound  against  a  fully  closed  short 
at  the  insulator  stack  that  has  electrically  isolated  the  front  end. 

Within  experimental  uncertainty,  the  K-shell  x-ray  yields  on  these  shots 
are  quite  comparable  with  earlier  gas  puff  experience.  In  fact  one  of  the  best 
yields  is  seen  on  shot  548  with  the  strongest  flashover  event.  So  while  it 
might  be  tempting  to  hope  that  the  correction  of  such  insulator  breakdown 
will  lead  to  a  better  performance,  the  studies  summarized  here  show  that, 
energetically  at  least,  the  flashover  is  too  little,  too  late  to  be  having  any 
serious  effect  on  the  K-shell  yield  in  these  events. 

The  experimental  implosions  were  examined  with  three  distinct  computa¬ 
tional  methods.  First  a  transmission  line  code  was  set  to  model  the  insulator 
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Volt«  and  Mpa  i  I  I  i  I  I  Volta  and  Aa^a 


breakdown  as  a  time  dependent  shunt  resistance  which  initiates  at  a  fixed 

8.00E+6  I—  ■■  I  — I.  I  ^  MM.-.  7,ooE^.s 

r  ^  Fig.  1 

7.00E+6 - f - >  '  _  _  ^ 


FoiLand  current  closer  to  shot 
conditions  .  Shuntto  37  Sm  at 
240  ns  ,w  ellbefore  peak  current. 
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Fig.  2 
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time  and  drops  it  resistance  to  a  specified  asymptote.  The  upstream  circuit 
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model  includes  the  monitors  for  tube  current  and  voltage,  while  downstream 
it  is  equipped  with  a  MITL  model  and  a  Z-flow  shunt  which  joins  a  post  hole 
convolute  model  and  a  PRS  load  model. 

In  the  figures  above  note  the  tube  current,  foil  current  and  voltage,  the 
flashover  current,  the  energy  into  the  flashover,  and  the  foil  kinetic  energy. 
An  early  flashover  is  shown  to  rob  foil  energy  while  a  late  one  does  not. 
Tube  voltage  suppression  occurs  to  the  same  value  for  equal  shunt  impedance 
at  either  timing,  but  this  is  not  a  strong  indicator.  In  Fig.  1,  are  the 
characteristics  of  an  early  breakdown  with  a  large  tube  current  surge;  in  Fig. 
2,  a  breakdown  in  keeping  with  the  experimental  timing  shows  virtually  no 
energy  loss  to  the  implosion. 

The  shunt  37.5mfl  impedance  used  here  was  found  to  be  about  the  best 
mimic  of  the  experimental  behavior  allowing  peak  shunted  currents  on  the 
order  of  2  to  3  MA,  and  the  decay  time  of  the  shunt  resistor  was  kept  short 
compared  to  the  other  timescales  in  the  problem  at  about  2ns.  Hence  within 
the  limitations  of  this  single  line  circuit  model,  the  flashover  event  is  captured 
reasonably  well. 

Fig.  3  Decade  Shot544  : 

0  bserved  Loss  C  unentFittD  a  C  ubi:  att=  159 OE -6  s 

I(t)  =  -534263294E-6*t3  +  7  55674  806E -4 +  -2  .28795080E -3  *t 


tine  (after  Qa^rx^)  [ne) 


loss  currait 


Alternatively,  using  the  observed  current  loss  from  shot  544  (c.f.  Fig. 
3)  instead  of  an  arbitrary  shunt  resistance,  the  free  variables  are  reduced  to 
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only  the  breakdown  time.  The  only  modification  of  the  equivalent  circuit  is 
to  lower  the  “output  voltage”  by  the  machine  impedance  times  the  observed 
shunt  current,  viz. 


Kc(^)  Voc{t)  -  Z,ghbs{i)  ,  (6-1) 

with  Iobs{t)  a  known  function  of  time  as  modeled  in  Fig.  3.  In  order  to 
examine  K-shell  yield  sensitivity  rather  than  simple  energy  transfer  the  ID 
code  DZAPP  was  also  applied  to  these  cases  and  equipped  with  a  standard 
equivalent  circuit  driver.  Again,  using  the  experimentally  observed  current 
loss,  the  shunt  timing  was  the  primary  factor  in  setting  the  expected  K-shell 
yield.  An  early  shunt  at  140  ns  quenched  the  yield,  while  a  late  one  (at  210 
ns  like  in  the  experiment)  showed  virtually  no  effect. 

Table  1.  Shunt  Effects  with  Shot  544  Loss  Current 


Shunt  Time 
[ns] 

Yt 

kJfcm 

kj/cm 

Current 

[MA] 

none 

52.0 

10.0 

5.45 

210 

50.0 

10.0 

5.45 

140 

38.0 

0.059 

4.75 

When  the  same  model  is  applied  in  2D,  the  Mach2  code  shows  the  same 
result.  These  two  cases  differ  only  in  that  Fig.  4  has  no  shunt,  while  Fig.  5 
uses  the  experimental  trace  and  timing.  The  change  in  K-shell  yield  is  almost 
imperceptible  even  though  the  late  time  currents  and  magnetic  energy  in  the 
front  end  differ  markedly. 


6.2  Conclusions 

From  the  foregoing  analysis  it  would  appear  that  DQ  was  not  hampered 
significantly  by  the  insulator  flashover  events,  even  though  the  yields  were 
not  much  improved  over  earlier  lower  diameter  nozzle  shots.  Insofar  as  the  2D 
calculation  showed  some  Rayleigh  Taylor  modulation  and  some  zippering  (in 
keeping  with  the  experiment)  these  effects  alone  were  not  enough  to  quench 
the  K-shell  yield.  This  is  not  to  say  that  azimuthal  perturbations  and  more 


pronounced  r-z  perturbations,  perhaps  brought  on  by  an  asymmetric  power 
flow  from  the  flashover  or  a  tube  misalignment,  could  not  account  for  the 
observed  performance.  Indeed  since  the  one  flashover  event  (548)  that  seems 
to  have  involved  both  front  and  back  insulators  showed  a  slightly  enhanced 
yield,  the  largest  tube  current  surge,  and  currents  rising  after  compression 
in  the  load  region,  there  is  a  strong  anticorrelation  of  flashover  with  poor 
yield  in  the  dataset!  The  fact  that  a  clean  strong  flashover  isolating  the  front 
end  might  have  helped  the  pinch  may  indicate  that  tube  misalignment  is  a 
good  candidate  problem  area  to  clear.  Perhaps  the  poor  yields  are  arising  in 
part  from  asymmetric  power  flow  causing  more  zippering  or  other  incoherent 
assemblies. 

However,  for  the  Decade  Quad  in  the  shot  series  examined  here,  what 
seems  clear  is  that  (i)  the  observed  yield  variations  do  not  correlate  well 
with  “strong”  or  “weak”  flashover  events,  and  (ii)  the  general  energetics 
and  observed  timings  of  the  flashover  events  do  not  imply  significantly  lower 
yield  expectations.  Moreover,  these  conclusions  are  largely  independent  of 
the  analysis  scheme  (OD,  ID,  or  2D)  used  to  examine  the  issue. 
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